ADA051067 


ELECTROMAGNETICS  LABORATORY 
TECHNICAL  REPORT  NO.  77-24 

December  1977 


MUTUAL  ADMITTANCE  BETWEEN  SLOTS  ON  A 
CYLINDER  OR  CONE 


S.  tf.  Lee 


R.  Mlttre 


ELECTROMAGNETICS  LABORATORY 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ENGINEERING  EXPERIMENT  STATION 
UNIVERSITY  OF  ILLINOIS  AT  URBANA-CHAMPAIGN 
URBANA.  ILLINOIS  61801 


Supported  by 

Contract  No.  N00019-77-C-0127 
Department  of  the  Navy 
Navy  Air  Systems  Command 
Washington,  D.C.  20361 


report  DOCUMEWTATIOH  PAGE  j BEF^E^c^gggfgg^oRM 

1.  neroRT  NUMten  li.  govt  accbuion  NoIJvNeciPiENT's  catalog  number 


title  (mt  SiMitI*) 


i^irUAL  ^MTTANCE  BETWEEN  gLOTS  ON  A CTLINDEiToR] 
ft  CONE  » " — " 


*UTHOWr<0~ 

fs.  wT/Lee  ”' 

IR./Mlttra 

•■  PERFONMINO  ONOANIZATI'SN  NAME  AND  ADDRESS 

Electromagnetics  Lahoratoryv/ 

Department  of  Engineering,  University  of  Illinois 
Urbana,  IL  61801  ^ 

n CONTNOLLINO  OFFICE  NAME  ANO  AOOSESS  (/L 

Department  of  the  Navy 
Naval  Air  Systems  Command 
Washington,  D.  C.  20361 

U.  MONlTodlNd  ASENCY  NAME  * AOORElVIf  Inm  CwiIMfllnl  6wc») 


Final  <ei^t« 

16  Nov  76  — 15  NOV  77 

t.  PEWr0WMINr.f.«O.  ■E>0«T  MUUI»ri. 

rEM»77-24«  UILU-ENG-77-2267 
TTfl^nfXPrETWSirrNuiJSeSflj^ 


19-77-C-9127 


ID.  enOORAM  element.  RROJECT,  TASK 
AREA  • WORK  UNIT  NUMRERS 


PR-SJ002 

l£££2S5J£ILJ__: 

^1.  NUMBER  OE  RAOEI 

m J 

TiriECURITVCLAMrToTlAA 


if.  distribution  statement  (*i  tAi«  rarwu 

Distribution  unlimited.  (Reduction  In  whole  or  in 
purpose  of  the  United  States  Government.) 


unclassified 

Is*.  declassification''oownoradino 
schedule 


part  Is  permitted  for  any 


Mr.  distribution  statement  (tl  «•  «S«lr«e«  Mlartf  In  Bf««S  SO,  II  tlihrmi  Inm  RtR«M> 


MS.  surrlementarv  notes 


Itt.  Kgy  WOWOI  fCiKIww  Oft  rtycM  f/ite  U m0t999Wf  i4milfr  >r  >!♦<* 


GTD 

ray  techniques 
conformal  array 
mutual  coupling 


slot 

high-frequency  diffraction 


fts^iSTW  ACT  (CmtIniM  ftw  wvff  914$  II  wftcftftfQr  iWwtliy  *y  Moct  numkf) 

In  $he  present  contractr-w^ studied  the  mutual  admittance  between  two  thin  slots 
on  a conducting  cylinder  or  a cone  by  ray  techniques.  Our  main  contribution  In 
the  development  of  an  approximate  asymptotic  Green's  function  for  the  surface 
function,  expressed  In  terms  of  Fock  functions.  Is  simple  to  evaluate,  and  gives 
excellent  numerical  results  when  compared  with  known  exact  solutions  and/or 
experimental  data.  This  report  contains  a brief  administrative  summary  plus 
three  attachments  vdilch  give  the  technical  details. 


DO  I JM*?*  1473  edition  ok  I NOV  f I IS  OBSOLETE 


\ unclassified 

SECURITY  CLASSiniCATION  OF  THIS  PACE  D*l»  BnnrmO 


^/(i  ^ i.() 


y 


UILU-ENG  77-2267 


Electromagnetics  Laboratory  Report  Mo.  77-24 


MUTUAL  ADMITTANCE  BETWEEN  SLOTS  ON  A 


CYLINDER  OR  CONE 


S.  W.  Lee 
R.  Mlttra 


Technical  Report 


December  1977 


Supported  by 

Contract  No.  N00019-77-C-0127 
Department  of  the  Navy 
Naval  Air  Systems  Command 
Washington,  D.C.  20361 
Mr.  James  Willis,  AIR-310B 
Contract  Manager 


UCCFSSION  tor  

NTIS 

poc  B fl  Sc., .on 

UNANNOUNCED 

JUSTI  ICC'i:  1 


OlSlRIByilOH/WAI'ABI'ITl: 
nT'* ' CiAL 


Electromagnetics  Laboratory 
Department  of  Electrical  Engineering 
Engineering  Experiment  Station 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 


Distribution  is  unlimited.  (Reproduction  in  whole  or  in  part  .Is  permitted 
for  any  purpose  of  the  United  States  Government.) 


□ n 


HIILWII  luywiJ  II , IMJIJU  UII 


■n**'<PPIiliPPiPPPP 


iii 


TABLE  OF  CONTENTS 


Page 

I.  INTRODUCTION 1 

II.  PERSONNEL 2 

III.  TECHNICAL  RESULTS  3 

IV.  PUBLICATIONS  AND  PRESENTATIONS 4 

ATTACHMENT  A:  MUTUAL  ADMITTANCE  BETWEEN  SLOTS  ON  A CYLINDER  5 

ATTACHMENT  B:  SIMPLE  APPROXIMATE  FORMULA  FOR  MUTUAL  ADMITTANCE  BETWEEN 

SLOTS  ON  A CYLINDER 107 

ATTACHMENT  C:  MUTUAL  ADMITTANCE  OF  SLOTS  ON  A CONE:  SOLUTION  BY  RAY 

TECHNIQUE 147 


mmmmm 


3. 

,1 


1 

I 

! 


\ 


9 

i 

\ 


I.  INTRODUCTION 


The  contract  N00019-77-C-0127  was  awarded  to  the  University  of  Illinois 
by  the  Naval  Air  Systems  Command  for  "Mutual  Admittance  Between  Slots  on  a 
Cylinder  or  Cone"  for  a one-year  period,  16  November  197C  to  15  November  1977. 
Mr.  J.  Willis  of  AIR-310B  is  the  contract  monitor. 

This  is  the  final  report  for  the  contract,  covering  personnel  (section  II), 
technical  results  (section  III  and  attachments) , publications  and  presentations 
(section  IV). 


The  attachments  have  their  own  pagination.  The  report  is  page  numbered  in  lower 
left-hand  and  right-hand  corners. 
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II.  PERSONNEL 


S.  W.  Lee,  Professor  of  Electrical  Engineering,  University  of  Illinois, 
principal  Investigator 

R.  Mittra,  Professor  of  Electrical  Engineering,  University  of  Illinois, 
principal  investigator 

J.  Boersma,  Professor  of  Mathematics,  Technological  University  of  Eindhoven, 
The  Netherlands,  and  visiting  professor.  University  of  Illinois 

V.  Krlchevsky,  Associate  Professor  of  Electrical  Engineering,  University  of 
Illinois 

S.  Safavl-Nalnl,  Research  Assistant,  University  of  Illinois 
L.  Grun,  Research  Assistant,  University  of  Illinois 

P.  Chang,  Computer  Programmer,  University  of  Illinois 
C.  L.  Law,  Computer  Programmer,  University  of  Illinois 


III.  TECHNICAL  RESULTS 


In  the  design  of  a slot  array  on  a conformed  surface,  a most  important 
parameter  is  the  mutual  admittance  between  two  slots.  In  the  present 
contract,  we  have  studied  the  following  problems  about  Yj^2* 

(a)  When  the  conformal  surface  is  a conducting  cylinder,  a GTD  solution 
of  Yj^2  bas  been  successfully  developed.  It  applies  to  cylinders  with  radius 
greater  than  one  wavelength,  and  gives  excellent  numerical  results  (error  is 
within  0.25  dB  in  magnitude  and  several  degrees  in  phase).  Details  are  given 
in  Attachment  A. 


(b)  For  the  slot  array  on  a cylinder,  a simple  approximate  solution  of 
Yj^2  derived.  It  is  generally  valid  when  the  separation  between  the  slots 
is  greater  than  two  wavelengths  (Attachment  B). 


(c)  The  GTD  solution  described  in  (a)  has  been  generalized,  so  that  it 
now  can  be  used  to  calculate  Yj^2  between  slots  on  a general  convex  conducting 
surface.  In  particular,  it  was  applied  to  slots  on  a cone,  and  the  numerical 
results  of  Yj^2  good  agreement  with  the  experimental  results  (Attachment 


IV.  PUBLICATIONS  AND  PRESENTATIONS 


(1)  S.  W.  Lee,  S.  Safavl-Nainl,  and  R.  Mittra,  "Mutual  Admittance  Between 
Slots  on  a Cylinder,"  University  of  Illinois,  Electromagnetics  Lab.  Tech. 
Rept.  77-8,  Urbana,  Illinois,  March  1977. 

(2)  S.  W.  Lee  and  S.  Safavl-Nalnl,  "Simple  Approximate  Formula  for  Mutual 
Admittance  Between  Slots  on  a Cylinder,"  University  of  Illinois, 
Electromagnetics  Lab.  Tech.  Rept.  77-13,  Urbana,  Illinois,  July  1977. 

(3)  S.  W.  Lee  and  R.  Mittra,  "Mutual  Admittance  Between  Slots  on  a Cylinder 
or  a Cone,"  Final  Report,  Contract  N00019-77-C-0127  (Including  the  above 
two  reports  as  attachments),  December  1977. 

(4)  S.  W.  Lee  and  S.  Safavl-Nalnl,  "Approximate  Asymptotic  Solution  of 
Surface  Field  Due  to  a Magnetic  Dipole  on  a Cylinder,"  to  appear  In 
IEEE  Trans.  Antenna  Propagat..  1978. 

(5)  S.  W.  Lee  and  S.  Safavl-Nalnl,  "An  Unexpected  Result  About  Surface  Rays," 
Digest  of  International  Symposium  on  Antennas  and  Propagat..  pp.  52-55, 
Stanford,  CA,  20-22  June  1977. 

(6)  S.  W.  Lee,  S.  Safavl-Nalnl,  and  R.  Mittra,  "Mutual  Admittance  Between 
Slots  on  a Cylinder,"  paper  presented  at  Program  Review  of  Naval  Air 
Systems  Command.  Hughes  Aircraft  Company,  Culver  City,  CA,  29  March  1977. 
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1.  INTRODUCTION 

In  the  design  of  a conformal  slot  array  on  the  surface  of  a 
conducting  cylinder,  the  calculation  of  the  mutual  admittance 
is  a crucial  step,  which  has  been  studied  extensively  in  recent  years. 
In  this  paper,  we  summarize,  in  a handbook  format,  all  of  the  final 
formulas  of  Y^^*  present  some  typical  numerical  data. 


2.  STATEMENT  OF  PROBLEM 

Referring  to  Figure  1,  two  Identical  slots,  circumferential  or 
axial,  are  located  on  the  surface  of  an  infinitely  long  cylinder.  The 
geometrical  parameters  are 

R « radius  of  the  cylinder  (2.1) 


I 


(a,b)  ■ dimensions  of  the  slot  along  (^,z)  directions  (a  is 


the  arc  length  along  the  cylinder)  (2.2) 

(Zq,R^q)  - center-to-center  distances  between  slots  (2.3) 

Sq  •Jzl  + (R^q)^  (2.4) 

0Q  » tan"^(zQ/R(}>Q)  (2.5) 


m 


The  problem  is  to  determine  the  mutual  admittance  between  these  two 
slots  when  kR  is  large. 
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First  let  us  define  mutual  admittance.  Throughout  this  work 
we  always  assume  that 

(i)  the  slots  are  thin,  and  (2.6a) 

(ii)  their  length  is  roughly  a half-wavelength.  (2.6b) 


Then  the  aperture  field  in  each  slot  can  be  adequately  approximated  by 
a simple  cosine  distribution,  which  is  the  so-called  "one-mode"  approximation. 
For  example,  if  slot  1 is  a circumferential  (lower  slot  in  Figure  la), 
its  aperture  field  under  the  "one-mode"  approximation  is  given  by 
(exp  + jdit  time  convention) 

E - Vj^e^  , H - (2.7a) 

where 


X X e. 


(2.7b) 


R(|i  . 


(2.7c) 


(Vj^.Ij^)  are  respectively  the  modal  (voltage,  current)  of  slot  1.  The 

mutual  admittance  Y^2  defined  by 

I. 


Y ■ Y 
12  *21 


‘21 


'l2  ■ vjvj  II, 


E2  X . dSj 


(2.8) 


where  l2^  is  the  Induced  current  in  slot  2 when  slot  1 is  excited  by 
a voltage  and  slot  2 is  short-circuited.  An  alternative  expression 
for  Y22  is 


(2.9) 


where 

A2  " aperture  of  slot  2 

•»  magnetic  field  when  slot  1 is  excited  with  voltage  and 

slot  2 is  covered  by  a perfect  conductor 

E2  ° electric  field  wl>cn  slot  2 is  excited  with  voltage  V2,  and 

slot  1 is  covered  by  a perfect  conductor. 
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Because  Hj^  = 121^2  z ~ ^^2^2’  ^ simple  matter  to  verify  that 

(2.8)  and- (2.9)  are  equivalent  [1]. 

There  Is  an  alternative  definition  of  mutual  admittance.  Instead 
of  (2.7),  a modal  voltage  (with  a bar)  may  be  defined  through  the 
expression  for  the  aperture  field  of  slot  1 as  follows: 

E = z f V,  cos  — y 
b 1 a ■’ 


(2.10a) 


or  equivalently 
fb/2 


fD/Z 

1 “ J 


(2.10b) 


Then  a different  mutual  admittance  is  defined  by  (2.9)  after  replacing 
(Vj^,V2)  by  (Vj^.V^).  It  can  be  easily  shown  that 

"l2'i^2  * ^2.11) 

Two  remarks  are  in  order:  (1)  In  the  limiting  case  that  b 0,  goes 
to  zero  as  b , whereas  approaches  a constant  independent  of  b. 

(li)  For  the  special  case  a « A/2  and  R •>  «,  it  is  that  is 
identical  to  the  mutual  Impedance  between  two  corresponding  dipoles 
calculated  by  the  classical  Carter's  method  (2],  [3],  [4],  (ill)  When 
the  slots  are  excited  by  waveguides  (transmission  lines),  one  often  uses 
^12  ^^12^'  ^tom  here  on,  we  will  concentrate  on  Y^^^  instead  of  Yj^^* 

The  mutual  admittance  defined  in  (2.8)  and  (2.9)  includes  the  self 
admittance  Y^^^  as  a special  case  which  occurs  when  two  slots  coincide. 

(All  the  formulas  of  Yj^2  Slven  in  this  paper,  except  for  the  one  in 
Section  4,  can  be  used  for  calculating  Y^^^  by  setting  (jip  ->  0 and  ■+  0.) 


This  solution  is  suitable  for  numerical  calculation  if  (i)  z^  < b for 
circumferential  slots,  and  z^  < a for  axial  slots,  (11)  kR  is  less  than 
20,  and  (ill)  the  medium  is  slightly  lossy  so  that  k has  a small 
(negative)  imaginary  part.  Based  on  this  solution,  extensive  numerical 
results  have  been  reported  by  Hughes  Aircraft  Company  at  Culver  City 
[7],  [8],  [9]. 


4.  EXACT  UI  MODAL  SOLUTION 


Under  the  one-mode  approximation,  another  exact  modal  solution  is 
given  in  (10].  This  solution  is  derived  from  the  Hughes  (SGP)  solution 
in  Section  3 by  a deformation  of  integration  contour  and  an  application 
of  the  Duncan  transform  (11).  THe  final  result  reads 


Circumferential  slots 


Y^2  " + .iii 


f ~ CCS  11,1^1 

C;  r.  ) eos  1/  vHn.k  )lt(u,!<  ) (!!•• 

7 0 ir-0  ' 

cos  ml,  ^ I 

® “ I 7 {"  R(ni,k  )i;'(m,k  ) sin  k z . dk 

npO  ‘•m  's.Jo  z 0 7 


+ j R(m,Jn)iJ'(m,Jn)e  ' dn| 


(4.1a) 


(4.1b) 


(4.1c) 


L m l ' / in  l 


2 , m - 0 
1 , m 0 


t|((in,k  ) is  defined  in  (3.2)  and  k in  (3.4) 
2 • 


Axial  slots 


8Y„  ® cos 


Y a — - 

12  TTkR 


« cos  m?-  .,  _ OK 

Or  0 , X -ik  z 

i —r «(m,k  )e  z 0 -5 

dpO  m ^0  N (k.R) 


f ♦(m, 

•'0 


jn)e  *^*0 


where  $(m,k^)  is  defined  in  (3.6) 


This  solution  is  valid  only  if  Zq  > b for  circumferential  slots  and 
Zq>  a for  axial  slots.  It  is  suitable  for  numerical  calculation  if 
kR  is  less  than  20. 


5.  ASYMPTOTIC  SOLUTION 


The  two  modal  solutions  given  in  Sections  3 and  4 are  based  on 
fields  in  the  Fourier  transform  domain.  An  alternative  calculation  of 
Yj^2  involves  the  field  in  the  spatial  domain,  namely, 

Circumferential  slots 

^12  “ ib  L ‘^^I'^'l  a ^1^^"=°® 

• A • A ' 


Axial  slots 


12  ab 


[ dyj^dZj^  I dy2dz2  [cos -g  z^^) {cos  ^(Z2  - ZQ)]gj,(8,e)  (5.2) 

'A  * A 


where  (y  , z ) - a typical  point  in  the  aperture  of  slot  n (n  - 1 or  2) . 


A^  « aperture  of  slot  n 


s = i^y2  - + (^2  ' 


e = tan'  [(Z2  - Zi)/(y2  " 


Several  versions  of  the  Green’s  functions  g.  and  g have  been  approximately 

(f  z 

determined  under  the  condition  that  kR  » 1.  They  are  listed  as  follows: 


OSU  Asymptotic  solution  [12]  [13] 
g^  'V  G(s)  [v(5)  sin^e  + (j^)u(5)  cos^e] 

-v  G(s)  (v(0  cos^e  + sin^e] 


(5.7) 


(5.8) 


PINY  Asymptotic  solution  [9]  [14] 

G(s)  v(O(sln^0  + (1  - 3 sin^0)] 

V > 

+ sec^0[u(C)  - Vj^(5)  sln^0] 

> 

g^  'V.  G(s)  V(5)  [cos^0  + (2  - 3 CO8^0)] 


(5.9) 


(5.10) 


UI  Asymptotic  solution  [15] 

g^  -v-  G(8)  ^v(O[sin^0  + CO820]  + (^)  u(5)(cos^0{l  " |^)  + ('^)  8in^0] 

+ j (1/2* kR/cos^0)  ^^^[v’(5)  8ln^9  + (tan^0  '^■^)'^'(^)  cos^0]J  (5.11) 
g^  • G(8)  v(C)(co8^0  - ^ CO820]  + (^)u(0[sin^9(l  " + (^)  cos^0] 

+ J(i^  kR/co8^0)“^^^[v'(C)cos^0  + (1  + J-)u’(C)8ln^0]j  (5.12) 


where 


k Yq  -jks 

^ V-  ’ ^0  - ^120.) 

? ■ (k  co8^0/2R^)^^^8 


(5.13) 


(5.14) 


The  Fock  functions,  u,  v,  etc.,  can  be  calculated  from  the  following 
two  representations: 


For  0 < C 0.7 


20 


8 


,(5)  .1-4  J^/\y2  + ii  53  * ^ ,-J./4^9«  . , ^ ^„-3^6 

u(5)  .1-4  + |i  5^  ^ 54^-J,M^9/2  . 3_,„^  ^ 33-2^6 

»1«)  -1+4  - S «'  - ^.-J’'«C»«  + 4.555  X lo-^s* 

v*(5)  '1'  ~ e"^3Tr/4  1/2  ^7^-2  „-1ir/4  7/2  - 2.485  x lo”^5^ 

8 -i  ^ 20  ^ ^ 1024  ® 5 


u’  (5)  ''-  f ^ ^2  ^ ^-jir/4  7/2  _ 

4 ^ 4 ^ 128  ^ 

For  0.7  < i < <« 

v(5)  = Y (t'  )“^e  " 


221  X 10"^5^ 


(5.15) 

(5.16) 

(5.17) 

(5.18) 

(5.19) 


n*l 

10  -^F^■ 

u(5)  = I e ” 

n«l 

v,(C)  = {1/1  Y e'^“'" 

n*l 

v'(e)  = i .-1’/'*^  {-1/1  “ a - j2{f.)(t'.)V"'- 


n«l 


10 


n n 
-Kt 


u'(5)  = .I'/S/:  {1/1  J (l  - J f {t  )e  ” . 


n»l 


(5.20) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 


where  - |tjj|exp(-jiT/3)  , t'^  - It'j,|exp(-jiT/3),  and 


n 

l'„l 

|t'l 

' n' 

1 

2.33811 

1.01879 

2 

4.08795 

3.24820 

3 

5.52056 

4.82010 

4 

6.78671 

6.16331 

5 

7.99413 

7.37218 

n 

1^1 

ii;i 

6 

9.02265 

8.48849 

7 

10.04017 

9.53545 

8 

11.00852 

10.52766 

9 

11.93602 

11.47506 

10 

12.82878 

12.38479 

It  has  been  verified  through  several  hundred  numerical  examples  that  the 


UI  asymptotic  solution  given  above  is  in  excellent  agreement  (within  a 


quarter  db  in  magnitude  and  a few  degrees  in  phase)  with  the  exact  model 


solution  for  all  slot  separations  (^q,Zq)  provided  that  kR  > 5. 


In  using  the  asymptotic  solutions  for  calculating  the  self  admittance 


care  must  be  exercised  in  avoiding  the  singularity  in  the  Green's 


function  which  occurs  at  s ■ 0.  A most  convenient  way  to  avoid  this 


apparent  difficulty  is  to  (1)  use  a large  nuiid>er  of  points  for  the  two 


surface  integrals  in  (5.1)  and  (5.2),  and  (li)  shift  slightly  the 


integration  nets  for  this  two  surface  integrals. 


6.  EXACT  PLANAR  SOLUTION 


In  the  limit  kR  the  Green's  function  of  the  UI  solution  in  (5.11) 


and  (5.12)  is  reduced  to 


g - G(a)(8ln^0  + ^ (2  - 3 8in^e)(l  - ^)] 


(6.1) 


gjj  ■ G(s)(cos^0  + ^ (2  - 3 CO8^0)(1  - ^)1  , 


(6.2) 


When  (6,1)  and  (6.2)  are  used  in  (5.1)  and  (5.2),  we  obtain  the  exact 


solution  (under  the  "one-mode"  approximation  of  course)  for  two  slots 


on  an  infinitely  large,  conducting  plane. 


7.  APPROXIMATE  SOLUTION 


Based  on  the  UI  asymptotic  solution,  a simple  approximate  solution, 


is  reported  in  [10],  l.e.. 


Circumferential  slots 


(S(b  sin  6)  C(a  sin  0)]^  g 

J * 


(7.1) 


Axial  slots 


Yi2  = - ^ tS(a  cos  6)  C(b  sin  6)]^ 


where 


• c(»)  ■ SM.  (Mm  . 
1 - (kx/ir) 


The  (simplified)  Green's  functions  i,  and  g are  given  by 

9 2 


- G(s)[v(0  (sln^  e + ^ cos  20]  + ^ U(0  coa^  0 
+ ju'(C)(»'2  kR  cos  ain'  0] 

ijj  • G(s)[v(C)  {cos^  0 - ^ cos  2o|  + ^ u(C)  sin^  o]. 


This  solution  gives  an  accurate  numerical  result  (within  several 

.0 


percent  in  magnitude  and  less  than  5 in  phase)  provided  that 
kR  > 10  and  the  slot  separation  la  greater  than  two  wavelengths. 


8.  CONCLUDING  REMARKS 


Based  on  extensive  numerical  data,  we  conclude  that  '*^12  (Including 


as  a special  case)  can  be  best  calculated  by 


(1)  Hughes  modal  solution  If  kR  < 5 and  Is  less  than  the  axial 
dimension  of  the  slot. 


(11)  UI  modal  solution  If  kR  < 5 and  Zq  is  greater  than  the  axial 
dimension  of  the  slot,  and 


(7.2) 


’7.3) 


(7.4) 

(7.5) 


(ill)  UI  asymptotic  solution  if  kR  > 5 for  all  slot  separations. 

If  several  percents  of  error  are  acceptable,  the  approximate  solution 
can  be  used  if  kR  > 10  and  the  slot  separation  is  greater  than  two 
wavelengths. 
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I Slot  Type  Dimension  Suggested  by 


A 

Circumf, 

0.9"  X 0.4" 

(f  - 9 GHz) 

Aerospace 

Hughes 

B 

Circumf, 

0.5X  X O.OIX 
(R  in  inch) 

Hansen 

C 

Axial 

0.4"  X 0.9" 

(f  * 9 GHz) 

Aerospace 

Hughes 

D 

Circumf. 

0.5X  X O.OIX 
(R  in  X) 

Hansen 

E 

Circumf, 

0.5X  X 0.2X 



Hansen 

I F ■ Axial  , 0.5X  X 0.2\ 

I 

I I 

: 1 


In  all  tables,  listed  in  (db,  phase  in  degree)  format  where 
db  = 20  (Yj2  in  mho).  In  all  figures,  the  normalized  phase 
of  Yj^2  is  equal  to  Arg(Y j2expjksQ) . 
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DATA  SET  A OF  MUTUAL  ADMITTANCE 

(1)  The  mutual  admittance  ¥^^2  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

* (Exact)  Hughes  modal  solution 

* (Exact)  UI  modal  solution 

* Ul  asymptotic  solution 

* OSU  asymptotic  solution 

* PINY  asymptotic  solution. 

The  parameters  are 

* Frequency:  f = 9 GHa,  k =»  4.7878  (inch)“^,  X = 1.3123" 

* Cylinder:  R 1.991"  unless  specified  otherwise 

* Slot  A:  Circumferential 

a =>  0.9"  =»  0.6858X 
b = 0.4"  * 0.3048X 


‘11' 


1.70747  X 10"^  mho 

.-3 


■55.35  db 


Y = 1.8155  X 10  mho 
g 


* Center-to-center  distance  between  two  slots  is  (R(||q,2q). 

(2)  Yj^2  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  logj^^  ^1^12^  mho). 

(3)  Data  are  presented  in 


TABLE  A-1: 


= 0 and  various  Zq 


A-2:  Zq  = 2"  and  various 
A- 3:  Zp  = 0 and  various  <(|q 
A-4: 


*'0 


0 and  various  z^. 


Figure  A-1:  Mutual  admittance  Yj^2  between  two  circumferential  slots 
as  a function  of  ^q* 

A-2  Mutual  admittance  Yj^2  between  two  circumferential  slots 
as  a function  of  Zq. 

A-3:  ^ cylinder  (UI  modal  solution)  and  that 

on  a plane  as  a function  of  Zq. 

A-A:  Yj^2  ^ cylinder  as  a function  of  the  radius  R of  the 


cylinder. 
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TABLE  A-4 


UI  SOLUTIONS  OF  Y,-  OF  SLOT  A FOR  ())  = 0 
li  o 


^0 

Modal 

Asymptotic 

z 

o 

Modal 

Asymptotic 

0.5" 

-62.62  db 

-62.54 

11" 

-84.06 

-84.06 

-72° 

-72° 

-70° 

-68° 

1 " 

-66.82 

-66.71 

12" 

-84.61 

-84.65 

155° 

155° 

15° 

18° 

Olt 

-71.78 

-71.66 

13" 

-85.12 

-85.20 

z 

l-» 

l-» 

o 

-116° 

o 

o 

o 

103° 

3" 

-7U.18 

-74.67 

14" 

-85.63 

-85.70 

-31° 

-30° 

-175° 

-172° 

4" 

-76,89 

-76.89 

15" 

-86.09 

-86.17 

54° 

54° 

-90° 

-86° 

5" 

-78.51 

-78.44 

16" 

-86.48 

-86.60 

139° 

141° 

-4° 

-1° 

6" 

-79.85 

-79.77 

17" 

-86.85 

-87.01 

-136 

-134 

81 

84 

■7** 

-80.94 

-80.88 

18" 

-87.24 

-87.38 

/ 

-51° 

-49° 

166° 

170° 

Q** 

-81.84 

-81.83 

20" 

-87.91 

-88.08 

O 

34° 

37° 

-24° 

-19° 

Q" 

-82.65 

-82.66 

30" 

-90,33 

-90.68 

y 

119° 

122° 

110° 

115° 

II 

10 

-83.40 

-83.40 

40" 

-91.95 

-92.46 

-156° 

-153° 

-115° 

-110° 

UI  ASYM 


. Pi-yj 


DATA  SET  B OF  MUTUAL  ADMITTANCE 


(1)  The  mutual  admittance  Yj^2  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

* (Exact)  UI  modal  solution 

* UI  asymptotic  solution. 

The  parameters  are 

* Frequency:  f = 9 GHz,  k » 4.787787  (inch)  X = 1.3123" 

* Cylinder:  R « 1.991"  , 3.777",  6" 

* Slot  B:  Circumferential 

a = 0.656168"  - 0.50X 
b « 0.013123"  * O.OIX 

* Center-to-center  distance  between  two  slots  is  (R(j)Q,ZQ). 

(2)  Y^2  listed  in  (db  value,  phase  in  degree),  where 


db  value  ® 20  logj^g  mho), 


(3)  Data  are  presented  in 


TABLE  B-1:  ({»q  = 0 and  various  Zq 

B-2:  Zq  = 2"  and  various  (J>q 
B-3:  Zq  “ 8"  and  various  (t>Q. 

B-A:  Comparison  of  Hughes  and  UI  solutions 


TABLE  B-2 


UI  SOLUTIONS  OF 


OF  SLOT  B FOR  z 

o 


2" 


3: 


25 


4»  ?C'Vt*Vf'/  - 


TABLE  B-4 


COMPARISOK  OF  HUGHES  AND  UI  SOLUTIONS 


45"  1 2" 


1 

1 - 1.991' 

1 

Hughes 

UI 

Modal 

Modal 

Aaymp 

-92.3  db 

-92 

-92.03 

rt 

-79' 

-79° 

-78° 

-96.5 

-96.31 

-96.28 

153° 

152° 

153° 

-112.02 

-111.5 

-111.56 

33° 

36° 

37° 

-117.08 

-116.13 

-116.35 

-6° 

-4° 

-1° 

-112.73 

-112.51 

-112.98 

91° 

92° 

93° 

Hughes 

Modal 


R - 3.777" 


UI 


Modal  Asynp 


8 -92.52 

•.rtO 


R = 6" 

Hughes 

UI 

Modal 

Modal 

Asymp 

-92.87 

-92.70 

-92.74 

o 

1 

-76° 

-76° 

-97.34 

-97.24 

-97.19 

156° 

157° 

157° 

-114.42 

-113.85 

-113.81 

44° 

46° 

46° 

-118.38 

,0 


-125.40 


-119.36  -119.33 
10°  10° 


DATA  SET  C OF  MUTUAL  ADMITTANCE 


(1)  The  mutual  admittance  ¥^^2  between  two  axial  slots  on  an  infinitely 


long  cylinder  is  calculated  from  the 

* (Exact)  UI  modal  solution 

* UI  asymptotic  solution 


The  parameters  are 

* Frequency;  f = 9 GHz,  k = 4.7877  (inch)  X = 1.3123" 

* Cylinder:  R = 1.991",  and  other  values 

* Slot  C:  Axial 

a = 0.4"  » 0.3048X 
b = 0.9"  = 0.6858X 

* Center-to-center  distance  between  two  slots  is  (R(j)^,?^). 

(2)  Yj^2  listed  in  (db  valun,  phase  in  degree),  where  db  value  = 20  logj^Q 

(|y^2I 

(3)  Data  are  presented  in 

TABLE  C-1:  <{|q  = 0,  F = ].991",  and  various  Zq. 

C-2:  z»  = 1.5",  R * 1.991",  and  various  ij>„ 

0 u. 

C-3:  = 0,  Zq  = 8",  and  various  R. 

Figure  C-1;  ^ cylinder  (UI  nodal  solution)  and  that  on 

a plane  as  a function  of  z^. 


TABLE  C-2 
_ OF  SLOT  C FOR  z 


^0 

Modal 

As3nnptotic 

0° 

-86.58  db 

-86.31 

151° 

149° 

-86.41 

-85.15 

30° 

-26° 

-38° 

-87.43 

-85.77 

60° 

o 

00 

72° 

-93.02 

-91.04 

90° 

169° 

156° 

infinitely  long  cylinder  from  the 


* (Exact)  UI  modal  solution 

* UI  asymtotlc  solution 
The  parameters  are 

* Cylinder;  R = IX,  2X,  4X,  lOX,  “ (planar) 

* Slot  D;  Circumferential 


a « 0.5X 
b - O.OIX 

* Center-to-center  distance  between  two  slots  is  (Rci  ,z  ) 

^o  o 

(2)  Y^2  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  ■ 20  logj^Q  0^12! 

(3)  Data  are  presented  in 


TABLE  D-1: 

<^0 

D-2: 

<^0 

D“3: 

^0 

D-4: 

z 

0 

D-5; 

z 

0 

D-6: 

z 

0 

0,  R " 2X  and  various  z 

I 

0 and  various  R and  z 

o 

0 and  various  R and  z 

0 

0 and  various  R and  (b 

^o 

IX  and  various  R and  d> 

o 

5X  and  various  R and  (b 

o 


TABLE  D-1 

UI  SOLUTIONS  OF  Y,,  OF  SLOT  D FOR  (b  - 0 and  R 

LZ  O 


9 


-115.18 

73° 


TABLE  D-2 

UI  ASYIIPTOTIC  SOLUTIONS  OF  Y,,  OF  SLOT  D FOR  (|) 


66° 


8A  ^50 


65° 


R = 2X 

R - 4X 

R - lOX 

81.51 

81.51 

81.51 

90° 

90° 

90° 

-98.56 

-99.15 

-99.51 

71° 

74° 

76° 

-103.84 

-104.63 

-105.13 

75° 

79° 

81° 

-106.96 

-107.92 

-108.52 

75° 

80° 

83° 

-109.16 

-110.25 

-110.94 

75° 

80° 

84° 

-110.85 

-112.05 

-112.81 

75° 

80° 

84° 

-112.21 

-113.51 

-114.34 

74° 

80° 

84° 

-113.35 

-114.74 

-115.63 

74° 

80° 

84° 

-114.33 

-115.80 

-116.75 

73° 

79° 

84° 

-115.18 

-116.72 

-117.73 

73° 

79° 

84° 

-115.94 

-117.55 

-118.61 

72° 

79° 

84° 

Planar 
(R  = ®) 


-99.76 

-,-,0 


-105.47 


-108.93 


-111.40 

c-,o 


-113.33 

87° 


-114.91 

f>r>0 


-116.25 

88° 


-117.40 

88° 


118.43 


119.34 


TABLE  D-6 

UI  ASYIIPTOTIC  SOLUTIONS  OF  OF  SLOT  D FOR  z 


\ 

R = IX 

R = 2X 

tl 

R = lOX 

^ \ 
o \ 

-108.90  db 

-111.13 

-112.73 

-115.49 

O 

o 

66° 

70° 

61° 

-25° 

-109.31 

-111.98 

-114.67 

-121.98 

ro 

o 

o 

61° 

54° 

6° 

39° 

-109.98 

-113.34 

-117.66 

-129.87 

30° 

53° 

29° 

-83° 

1 

O 

-111.47 

-116.23 

-123.42 

-141.86 

45° 

37° 

-26° 

91° 

-72° 

-113.50 

-119.88 

-130.02 

-153.21 

60° 

14° 

-99° 

o 

1 

159° 

DATA  SET  E OF  MUTUAL  ADMITTANCE 


(1)  The  mutual  admittance  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

* UI  asymptotic  solution 
The  parameters  are 

*Cy Under:  R = IX,  2X,  4X,  lOX 
*Slot  E:  Circumferential 
a = 0.5X 
b = 0.2X 

*Center-to-center  distance  between  two  slots  is  (Ri  ,z  ) 

o o 

(2)  Y^2  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  logj^g  ^1^12^ 


(3)  Data  are  presented  in 
TABLE  E-1: 

E-2: 

E-3: 

E-4: 

E-5: 

E-6: 


z =0,  various  ‘t>  and  R 
o o 

z = 0.5X,  various  •I'  and  R 
0 o 

z = IX,  various  <{'  and  R 
o o 

z ■ 2X,  various  <t>  and  R 
o o 

z = 4X,  various  and  R 


z = 8X,  various  4’  and  R 


E-7:  Comparison  of  UI  asymptotic  and  UI  modal  solutions 

E-8:  Comparison  of  UI  asymptotic  and  UI  modal  solutions 

Figure  E-1:  Mutual  admittance  Y 2 between  two  circumferential  slots 
as  a function  of  i{>g. 

E-2:  Mutual  admittance  Y2^2  between  two  circumferential  slots 
as  a function  of  (^q. 

E-3:  Mutual  admittance  Yj^2  between  two  circumferential  slots 
as  a function  of  ())g. 

E-4:  Mutual  admittance  Yj^2  between  two  circumferential  slots 
as  a function  of  z^. 

E-5:  Mutual  admittance  Yj2  between  two  circumferential  slots 
as  a function  of  Zg. 

40 


52 


. ah  a 


TABLE  E-1 


TABLE  E-3 

UI  ASYMPTOTIC  SOLUTIONS  OF  Y 


R = IX 


-72.28  db 
68° 


-74.71 

49° 


-77.44 

26° 


-82.65 

-24^ 


-88.42 

-90° 


II 

>- 

It 

o 

II 

-73.34 

1 

-73.92 

-74.28 

73^ 

76° 

78° 

-75.33 

-80.15 

55° 

9° 

105° 

1 

00 

o 

ou> 

-91.02 

-161° 

-110.75 

-116° 

-86.49 

-100.56 

-122.14 

-76° 

-20° 

-18° 

-95.69 

-112.38 

-136.13 

132° 

-22 

111*^ 

-104.13 

-122.57 

-148.13 

-42° 

-41° 

-124° 

TABLE  E-4 


I 


UI  ASYMPTOTIC  SOLUTIONS  OF  Y,.  FOR  z = 2X 

12  o 


R = IX 

R = 2X 

11 

R = lOX 

0° 

-77.24  db 

-78.67 

-79.46 

-79.96 

70° 

76° 

80° 

82° 

O 

O 

-77.52 

-79.44 

-81.77 

-89.55 

1 

67° 

65° 

37° 

-148° 

o 

o 

CM 

-78.37 

-81.60 

-87.38 

-103.31 

57° 

31° 

-79° 

76° 

30° 

-79.73 

-84.81 

-94.18 

-114.68 

42° 

-21° 

112° 

-144° 

0 

45 

-82.59 

-90.76 

-104.37 

-128.88 

9^ 

-130° 

162° 

15 

60° 

-86.17 

-97.24 

-113.88 

-141.23 

-35° 

94 

177° 

153° 

Mutual  admittance 


1 


< 

Q 

U1 

O 

<n 

s 

< 

I 

CL 

2 

* 

h- 

o 

U 

>• 

(f) 

X 

o 

< 
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UJ 

< 

o: 

UI 

N 

CL 

_j 

H 

a. 

H 

Z) 
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< 

3 

oc 

o 

i 

0 

z 

1 

Mutual  admittance  between  two  circumferential  slots 

as  a function  of  z_- 


0.5  X 


igure  E-5.  Mutual  admittance  between  two  circumferential  slots 

as  a function  of  z^. 


DATA  SET  F OF  MUTUAL  ADMITTANCE 


(1)  The  mutual  admittance  between  two  axial  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

*UI  asymptotic  solution 
The  parameters  are 

^Cylinder:  R = lA,  2A,  4A,  lOA 
*Slot  F:  Axial 

a = 0.2A 
b = 0.5A 

*Center-to-center  distance  between  two  slots  is  (Rij)  , z ) 

o 0 

(2)  Y^2  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  0^2.2^  mho). 


(3)  Data  are  presented  in 
TABLE  F-l 


F-2: 

2 “ 
0 

0.5A 

, various  * and  R 
0 

F-3: 

z = 
0 

lA, 

various  di  and  R 
0 

F-4 : 

2 = 
o 

2A, 

various  A and  R 
0 

F-5: 

2 = 
O 

4A, 

various  A and  R 
o 

F-6: 

2 = 
O 

<< 

CO 

various  i and  R 
o 

z - 0»  various  (j)  and  R 
0 o 


66 


F-7:  Comparison  of  UI  asymptotic  and  U1  modal  solutions 
F-8:  Comparison  of  UI  asymptotic  and  UI  modal  solutions 
F-9:  Comparison  of  asvmntotic  solutions 

Figure  F-l:  Mutual  admittance  Yj^2  between  two  axial  slots 
as  a function  of 

F-2 : Mutual  admittance  Y^2  between  two  axial  slots 
as  a function  of 

F-3:  Mutual  admittance  Yj^2  between  two  axial  slots 
as  a function  of  4»q. 
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TABLE  F-5 


TABLE  F-7 


COMPARISON  OF  UI  ASYMPTOTIC  AND  UI  MODAL  SOLUTIONS 


z 

0 

R = IX 

R = 2X 

^0 

Modal 

Asym. 

Modal 

Asym. 

IX 

0° 

-87.06  db 
-171° 

-86.65 

-173° 

-86.83 

-172° 

-86.63 

-172° 

10° 

-87.69 

176° 

-87.35 

172° 

-88.23 

139° 

-87.92 

134° 

20° 

-88.91 

139° 

-88.37 

128° 

-87.64 

35° 

-86.51 

26° 

30“ 

-89.40 

85° 

-88.07 

70^ 

-87.01 

-72° 

-85.77 

-81° 

A 5° 

-89.19 

2° 

-87.32 

-15° 

-88.67 

119° 

-87.30 

109° 

60° 

-89.84 

-83° 

-87.72 

-99° 

-91.86 

-61° 

-90.36 

-72° 

61 


73 


R = IX 


R = 2X 


Plan. 


! Modal 


0.5X 


IX 


-87.06 

-171° 


Asym. 


-70. 1^'  db 
25° 


-86.65 

-173° 


Modal 


-86.83 

-172° 


Asym. 


-70.11 

25° 


-86.63 

-172° 


(Exact) 


-70.08 

26° 


2X 


-99.97 

-17A° 


-99.37 

-177° 


-99.61 

-176° 


-99.34 

-176° 


-99.32 

-176° 


4X 


-112.43 

-175° 


-111.56 

-178° 


-111.93 

-177° 


-111.54 

-178° 


-111.52 

-178° 


8X 


-124.33 

-174° 


-123.63 

-179° 


-124.12 

-177° 


-123.61 

-179° 


-123.60 

-179° 


COMPARISON  OF 


^0 

R = 2A 

^0 

UI  Asym. 

PINY 

-99.34  db 

-99.42 

0° 

-176° 

-172° 

-100.00 

-99.93 

10° 

157° 

164° 

-99.39 

-99.71 

2\ 

Ni 

O 

o 

85° 

98° 

-97.05 

-97.85 

o 

O 

4° 

17° 

-95.40 

-96.16 

45° 

-126° 

-115° 

-96.10 

-96.68 

O 

O 

89° 

97° 

.E  F-9 

ASYMPTOTIC  SOLUTIONS 


R 

= lOA 

OSU 

UI  Asym. 

PINY 

OSU 

-105.44 

-99.33 

-99.41 

-105.42 

-172° 

-176° 

-172° 

-172° 

-105.37 

-92.2 

-92.51 

-92.53 

152° 

1 

o 

-142° 

-143° 

-101.89 

-92.24 

-92.46 

-92.45 

00 

O 

62° 

64° 

65° 

-98.23 

-96.^9 

-96.30 

-96.30 

4° 

-163° 

-161° 

-160° 

-96.09 

-103.04 

-103.25 

-103.25 

-120° 

-7° 

-5° 

-4° 

-96.6 

-110.19 

-110.41 

-110.41 

96° 

129° 

131° 

131° 

63 


75 


0.5X  R=2T0I0X 


NORMALIZED  PHASE 


This  appendix  contains  the  program  listing  of  all  solutions,  f 


pvnpnf  t-hp  pvpr‘^  Hnehps  mndal  soliUinn.  HI  rpiisspH  In  ^hp  rpvf. 


PK'OGK'AM  C YORPL  < INPUT  f t API: 6 o OUTPUT^  f APE7=QUTPUT » TAPES^'INPUt ) 


/j\  ^ ^ ^ ^ ^ ^ tV^  ^ ^ ^ * * *t'  ^ ^ ^ ^ ^ ^ ^ ^ •t'  ^ ^ ^ ^ ^ ^ <n  * ^ ^ ^ ^ ^ ^ ^ ^ * *,^ ‘ 


MUTUAL  ADMITTANeE  OF  SLOTS  ON  A CYL.tNI)L-R  OR  PLANE 


S.  W.  LEE 
S,  SAFAVI-NAINI 
P,  CHANG 
C . 1. . LAW 


IiATE : 


0/:l0/7/ 


UNIUERGITY  OF  ILLINOIS 


^ ^ ^ ^ ^ ^ <*p  ^ ^ ^ ♦*  ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 


C THIS  IS  A MOD I Fi ED  UERSION  OF  THE  PROGRAMS  BOLOIHO  MUTUAL 
C ADHITTANCfc  OF  SLOTS  ON  A CYLINDER  OR  ON  A PLANE. 

C THS  PROGRAM  IS  MODIFIED  TO  RUN  UNDER  CDC  CYDER  170  SERIES 
C COMPUTER  SYSTEMS. 

G THE  ORIGINAL  VERSION  WAS  REPORTED  IN  “ELECTROMAGNETICS 
C LADORATORY  (ECHNICAL  REPORI  N(!.7Y--e"  DAIE:  MARCH  I97Y. 

C. 

C THIS  PROGRAM  CONTAINS  7 HE  FOLLOWING  SODPROGRAMS: 

C 


(1)  PLANAR 

lA)  EXACT 

IB)  APPROXIMATE 


(PLANAR) 

(APROXl  X aPR0X;7) 


(2)Cylindri(;al 

2A)  UI  ASYMPIOTIC 
2B)  UI  EXACT  MODAL 
20  OSU  ASYMPTOTIC 
2D)  PINY  ASYMPinilC 
20  APPROXIMATE 


(CYL.INII) 

(PROGl  S PR002) 
(CYLtND) 

(CYLtND) 

(APROXl  f.  APR0X2; 


2F)  UI(2)  ASYMPTOriC  (CYLIND) 


C EACH  SUBPROGRAM  IS  A SUBROUTINE  TO  THE  MAIN  I'ROORAM. 

C 

ccccccccccccccccccccccf;cc(;ccc(u;ccccccccc(:ccct.ccccccci.'i'(;;ccc(:cccr 

c 

80 


JT’? 


nrinnoon  ooorjoonnonrjoorsorjn 


C THE  INPUT  FORMAf  IS  AS  FOLLOW  J 

C ♦.♦FCiR'MAT**  ' 

C 

C — (l)/(2)  . ♦ AlO' 

■C  -~~lA/lB/2A.L=:t<,"2C/2n/2E/2l^  ♦ 2Ai6 

C — AXIAL/CIRCUNi-ERENflAL  ♦ A0 

C ■— FfeoUEi^dY  ♦ ois.o 

C —SLOT  DIMENSION  <A»B)  ♦ &ei5.0 

C 

C IF  USING  APPROXIMATE  OR  UI  EXACT  MODAL  JUMP  TO 
C THE  NEXT  SECTION 

C —INTEGRATION  GRID  <IPA»IPB)  ♦ 215 

C 


C IF  USING  OTHER  THAN  UI  EXACt  MODAL  JUMP  tO  THE  NEXt  SECTION 
C (UI  EXACT  MODAL  IS  USING  TRAPEZOIDAL  RULE  FOR  iNtEGRATION. 

C NCYCLE=N0.  OF  SUBSECTIONS  BETWEEN  ANY  TWO  SUCCESSIUE  ZEROS 
C OF  INTEGRAND 

C MMAX  IS  THE  MAXIMUM  NO.  OF  TERMS  WHICH  HAS  BEEN  USED  IN 
C CALCULATION  OF  INFINITE  SERIES) 

C — NCYCLE»MMAX  ♦ 215 

C 

C THE  NORMALIZATION  FACTOR  (YU) 

C —NORMALIZATION  FACTOR  # GIS.O 


INPUT  ZO 

— TOTAL  NO.  OF  ELEMENTS  IN  ZO  (MMAX. -20) 
—20(1) 

“-Z0(2) 


* 12 

♦ Gt5.0 
t GtO.O 

# (U5.() 


IF  USING  PLANAR  “(D*  INPUT  THE  FOLLOWING  »IF  NOT 
JUMP  TO  THE  NEXT  SECTION 

—TOTAL  NO.  FO  ELEMENTS  IN  YO  (MAX. =20) 

— YO(l) 

— Y0(2) 


* 12 

♦ G15.0 

♦ G15.0 

♦ G15.0 


IF  USING  CYLINDER  ‘(Z)"  INPUT  THE  FOLl.OWINGf 
—TOTAL  NO.  OF  ELEMENTS  IN  PHI  (MAX.  =20) 

— PHKl) 

— PHI(2) 

— — ^ 

— TOTAL  NO.  OF  ELEMENTS  IN  RADIUS  (MAX.  =20) 
— RADIUS(l) 

— RADIUS(2) 


♦ 12 

♦ G15.0 

♦ G15.0 

♦ GtS.O 

♦ 12 

# G15.() 

♦ G15.0 

* GIS.O 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
IMPLICIT  COMPLEX  (C.H.Z) .REAL(A-B»D-Gf P-Y) 

REAL  PHI ( 20 ) » RADIUS ( 20 ) > Z ( 20 ) » YP ( 20 > 

REAL  TN(10)rTNPI(10) 

INTEGER  PINY»OSUrUI»TEbT 

81 


69 


9099 


8083 


8002 


LOGICAL;  eUHy AXIAL 
real  ^HAG;y:>0 

real  gone  y CdNFl y GdNF2  y KO  y ZO 
Cd'MMdN  PI  y tZi  y f Z2 y f Yi  y TY^  LRy  THETHA 
COMMON/DAI AI/TN  y f NPI y RHO  y Cl y C2  y F2  y lOP  y CC  yRAfiNj DEG 
CdHHdN/£iATA4/AGd2  y SN2  y TN2 y R2  y ACONi  y AddN2  y.dO.i  y diO 
GOMMdf'j/DAtAS/AO  y gOy  ZO  y YO 

CdHHON/DATAS/Kd  y NDYCLE  y PHIO  y ZO  y Yl  1 y MHAX  y A ftf 
DATA  TN/2 ♦ 33811  y 4 ♦ 00795 y 5 * 5215vi) y<S ♦ 7867:1  y 7 4 99417 r 
<1  9 4 02265  y 10  ♦ 6^<017 y 11 . 00852  y 1 3 ♦ 93602 y 1 2 , 82878/ 

DA  lA  TNP I / 1 4 0 1 8 79  y 3 4 24820  y 4 4 820 1 0 y 6.4 1 6331  y 7 4,372 i0:y 
i|.  8 4 48849  i-  9 ♦ 53545 y 10.4  52766  y il , 47506  y .12  ♦ :38479/ 

DATA  UI/O/yOiiiU/O/yPINY/O/yTYPEl/lOHEXACt  /,y 
$ PR02/10HCYUlNDRICA./yPR01/10HPLANAR  Ay  ' 
i TYPE2/10HAPPR(iXIMAT/  y TVPE3/10I  IU  I.  EXACT  B/ r 
$ TYPE4/.10H0I  A8YMPTd/yTYF’E5/10H08U  ABYMPf/y 
$ rYPE6/10HPINY  A8YMP/  y TYPE -VlOHUil  ( 2 ) A8 YM/  y 
OONFl  / 1 OHAXI AL  / y C0NF2/10HCIRCUMFERE/ 

DATA  IPAA/l/yIPBD/1/ 

PI==4)KATAN(J  .EO) 

RADN“IM/180,  $ DEG=;ie0  4/PI 

AC(:)N1~"187  4/64.  $ AC0N2.'=>24/34 

C01"<0.yl.)  3!  C10==(l4y04) 


DATA 


DATA 


ACQNl  “187  4 /64 . $ AC0N2i 
C01"<0.yl.)  3!  C10==(l4y' 

WRITER  6 y 555) 

READ ( 5 * 9099 ) PRO 
FORMAT (AlO) 

.I.F < PIO.)  4 EO 4 PROl ) GOTO  8082 
IF  ( PRO  4 E(3 4 PROP ) 001  0 0083 
WRITE<7y8084) 

8T0P 
IPLAN™2 
WRITE(6y666) 

GOT08093 
I PLAN  ^!. 

URlTE(6v  777) 

READ<5y80.-n  ) TYPEyTYPEE 
READ ( 5 y 909? ) CONF 
READ<5y8086)XK 
READ  t 5 y 8087 ) AA  y BD 
FORMAT (201540) 

IF  ( TYT'i:  4 i.:0 . 1 YPEt ) 0010  809/ 

IF  ( TYPE  4 EO  4 7 rP!L2 ) GOTO  0099 
I F ( TYPE  4 EO  4 T YPE3  ) SO  ? 0 ‘.•'000 
i F ( TYPE  4 EO  4 T YPE 4 ) 00 1 (1  809  / 
TF ( TYPE 4 EO 4 T YPE5 ) GOTO  809 / 
IF ( TYPE 4 EO 4 TYPE6 ) GOTO  809/ 

I F ( TYPE  4 EO 4 T YPE / ) GO  TO  (-09 7 
WR.i;TE(7yy084) 


809/ 

809/ 

H0<7/ 


9000 


8097 


READ ( i y 8098 1 NCYCLE  y MMAX 
(vOrO  8099 

K£Ai:M'5yB09>!;  'PAAv  ll'DB 


82 


70 


809*? 

8092 

8095 

8094 
' 2230 

8089 

8096 

9092 

9093 
771 


7 72 
773 

880 


999 

113 
1 12 


111 


333 

I 14 
414 
I 15 


‘8L0T  DIMENSION 
$ • t B< ALONG  Z)=  'ffO.S) 

WRITE( 6^773; 

URITE(6»772) 

IF  ( CONF . EO . CONI-  I ) GOTO  1 1 4 
URITE<6»333) 

FORMA  r < 1 X F * * • > 3X » ■ SLO  T ORIENTATION 

GOTO  115 

URITE(6f444) 

FORMAT ( IX  F ■ X! " F 3X » • SLOT  ORIENT  AI  lUN 
WRl rE(6F773) 


71 


READ(5f8086)Y,11 

ZYl=GMPLX,(YllF6iE0)  ■ - 

READ  (Si!  80914  NbZ 
DO  8092  I^lvNIiZ 
read (5f 8086) Zti) 

IF(IPLAN-..Lt.2)G8tO  2230 
READ  (S.F  8091  )NDPHt 
DO  O0$5  I^iFNIiRHI 
READ^5f8006)PHI<I) 
ftEAD(5F869i)NDR 
DO  80^4  i=lFNDR 
' READ.<5f8086) RADIUS! i) 

Goto  8096 
READ (5f 8091 )NDY 
DO  8089  I-1fNDY 
READ<5F8086)YP(r) 

IF<  CONFi . HO . CONF ) GOTO  9092 
AXIAL=. FALSE. 

CUM=.TRUE. 

GOTO  9093 
AXIAL=.TRIJE. 

Cl)M=.  FALSE. 

URITE<6f555) 

URnE<6F771)TYPEvTYPEE 

FQRMAT(///1Xf20<  "XO/IXf  •♦VIXf  *!»:“  f3Xf  "METHOD  OF  SOLUI  tON 
$ f4Xf2A10/1Xf  "JK'/IXfZO!  ■*" ) ) 

WRITE! 6 F 772) 

F0RMAT(/1Xf20("X!")/1Xf*)|c") 

FORMAT, (IXf  ’**/1Xf20<  “#* ) ) 

WRITE! 6 F 888) XK 
FORMAT ! I X r • ,t " » 3X  F * FREQUENCY 
WRITE !6f 773) 

URITE!6f772) 

IF!IPLAN.EQ.l)GOTa  113 
WRTTE!6f999)PF!02 
FORMAT  < 1 X F • X:  • F 3X  F " GOMETRY  .‘ 

GOTO  1.12 

WRITE!6F999)PR0t 
WRITE!6f773) 

URITE!6f772) 

WRITE!6f111)AA.PB 
F0RMAT!1Xf  "X!"  fSXf  ' 


'fE14.6) 


fAIO) 


A!ALDNG  PH.1) 


CIRCIJMFERENItAI  " ) 


AXIAL" ) 


83 


I 

I 


7.75 

CtSSl 

^ < f'  X*  - 

; ; :f 82 


I&3 

.4384 

885 


9094 


9095 


9097 

9096 

8081 

8084 

8086 

8091 


IF  of  YPE  i EQ  i'fyplS ) GOf 6 774 
UR’itE(;6»772) 

IF  < f YFiE  ♦ eO  ♦ f VPE3  XGOf 6 775 
(JRif  E (6 » 99 ) IPAA  f'iPBB 
FORHAT ( i X f • * S 3X i * I NTEGRAT ION  GRI D 
GOTO  776 

WRif  E ( 6 r.88i  ) NCYGEE rMMAX 

FORMAT  ( IX » " i|t ' » 3X»  •NOYOLE^  * »T5» " 


WRI-TE(6»773) 


WRrrE(6f772X 


■;»  i3f  • X*  V:I^;); 

? MHAX=^  * 


WRBTE(6»8S2)Y11 
F0RflAt(lXy  “>l:%3X»  "NORMALIZATION  FACTOR  t Yll~  %F7i2> 
WRlfE(6»773) 


WRITE (6>a83) 

FORMAT(/7/lOXr "$$$$$  DATA  OUTPUT  $$$$$•//) 

FORMAT ( /2X y "PHi-  •i»F7,2y'  <DEG>  t Z0~  'yF7.3y 
$ " i RADIUS-  "yE14.5) 


FORMAT (/2Xy"  Y=  *»F9.4»"  i Z-  "fFZ^S) 

IF ( TYPE . EO . TYPE3 ) GOTO  9095 
IF  ( TYPE ♦ EO ♦ TYPE4 ) IOP= 1 
I F ( TYPE . EO ♦ T YPE5 ) I 0P=2 


IF ( TYPE . EO ♦ TYPE6 ) I0P=3 
IF < TYPE ♦ EO » I YPE7 ) IOpu-4 
IF (CUM) GOTO  9094 
A--DB 
B=AA 


IPA-IPBB 

1PB--IPAA 


GOTO  9091 

A--AA 

B-=DB 

IPA=IPAA 
IPB^IPBB 
GOTO  9091 

A=AA  $ D"HD  $ KO--XK 
DO  9096  I=lyNDR 
RH0~RAD1US(I) 

DO  9096  Il~-tvNDZ 
ZO=Z(II) 

DO  9096  .1 1 1-=t.NDPHI 


15) 


PHIO-PHI(III))KRADN 
WR IT E ( 6 y 88  4 ) PH I( 1 1 1 ) , ZO  V RHO 
IF ( CUNFl . EO • CONF ) GOTO  9097 
CALL  PR0G2  ( RHO  r MAG » PHAST,  y DB » PHN ) 

GOTO  9096 

CALL  PROG 1 ( RHO » MAG  y PHASE  y DB » PHN ) 

WR I I E ( 6 y 78 ) MAG  y PHASE  y DB  y r’HN 
STOP 

FORMAT(AIOyAlO) 

FORMAKOXy  ERROR  y PLEASE  CHECK  YOUR  INPUT  AGAIN") 

FORMAT (61 5.0) 

1-ORMA  l < 12) 


84 


72 


■ 


8098 

C 

9091 


FdRMAT(2I5) 


13 


IGUMSI 
|A5(lAii  =2 
iR^yOTiGUM)  IGUMi^2 

nri  ^XIAU  i AJ<I AL=1 

IJP  30  IJ~ICUM>j:AXrAL 

IFCI4,EQ,2)Gdfd  fi 

A0=A!|(XK 

B6=B*Xk 

GO  TO  6 

SAVE=A 

A=iB 

B=SAVE 

A0=A)KXK 

bo^bjkxk 

ITEMP=IPB 

IPB=:IPA 

IPA=,TTEMP 

WIDTH1=A0/J;PA 

WIDTH2=B0/IPB 

Cl =CEXP ( CMPLX  < 0 , EO » -PI/3 , 

F2=SnRT<PI) 

Yl=~A0/2,~UiriTHl/2. 
21=-B0/2 , -UIDTH2/2 . 
IF(IPLAN.EQ«1)G0T0  7 
GOTO  8 
NDR=1 

RABIUSCI )=! ,E20 
CONTINUE 

DO  50  IRAD=1,NBR 
RHO=XK)(!RADIlJS(lRAn) 

DO  60  J2--l,Nn2 
ZO=Z<JZ)*XK 

IF<  IPLAN . EQ . 1 ) NriPHI=NnY 
DO  70  IY=1,NDPHI 
IF ( IPLAN, EQ.l)  00  TO  13 
YO=RHO)t:PHI(IY))»cRADN 

Y=YO/Jk^^^’'^‘^*®‘^^  Y0=0.00d 

GO  TO  lA 
Y=YP(IY) 

Y0=Y*XK 

CONTINUE 

IF(TYPE,EfJ.TYPE2)G0T0  17 
Y2=Y0-A0/2.-WiriTHl/2. 
Z2-Z0-B0/2 , - UIUTH:V2  . 

ZSUM=0 , 

DO  80  K=^1,IPA 
TYl=Yl+UIDrH.I  *K 


DO  90  L=lr.TPB 
fzi-:iZl  iW;ibTH25!CU 
DO  100  H~  U,;CPA 
I YS^YSf  W ) riTl  il,*M 
DO  11.0 

rz2“7;;'}wiDri-i2*N 

i?-gOR r ( ^ TV 2~‘TYt ) ( 7 Z2“TZ1 ) t*2 ) 
fHETllA-ATANSiX  (172-12:1. ) i>  ( fY2- TYl  > ) 

AC02-C0S  ( THETHA ) 
k2~Rtrt 

XF(iPLAN.F:n.;l,)  GO  TO  AOO 
GN2«3:i.fl  (THETHA.  **2 
rN2  -rAH(Tl-iETHA):t!!i<2 
tAi  L (;:VL'i;Mli(.tJi.ZGUM) 

GO  TO  !.  10 

OAL.L  P--.f-iNAR(.CJ»ZSUM)  • 

CONT.Tf  L!E 
CONTTHUE 
OON  r.(  NUb. 

CO.Jr.bNlJE 

ZY2- ( W i:  (I  rn  I >^0  IDl  H2 ) **2*  ( -2 . ) / ( AO*DO ) 
fH  (.T PLAN. F.0. 1)00 1'O  2/. 

WRJTE(  Av08*1)PH:1  ( T Y ) » 2 ( J2 ) f RADIUB(  TRAD) 
(SOTO  10 

WRJ.rE(6.000)YP(.i:Y)  >2(JZ) 

Gorn  10 

ITdPLAN.EC^.DOOTO  15 

WR [ TE  < A . 00 1) PHI ( 1 Y ) y 2 ( JZ ) » RADIUS ( IRAD ) 


WRI I E < 6 . 005 ) YP ( I Y ) y 2 ( JZ ) 

IF ( (:;ONF . EO . (JONFl ) CiOTO  16 
CALL  APROXl < AO y DO y 20 y YO y PI y 2Y2 ) 

GOTO  10 

CALL  APR0X2<A0yD0y20y Y0yPly2v2) 
riAG-C;ABS(2Y2) 

PHASE-ATAN2(ATMA0(2Y2) yREAL(2Y2)  dcDEO 
ZEXPON  = CEXP  ( CMPL  X ( 0 . EO  y SORT  < Y0)K)|t2  bZ0*!t:2 ) ) ) 

ZPR0D=2Y2!K2EXP0N 

PHN~ATAN2 ( AIMAG  < 2PR0D ) y REAL ( ZPROD ) ) ♦DEG 
I1B=20 . ♦ALOGIO < CABS  ( ZY2/ZY1 ) ) 

WRI PE ( 6 y 70 ) MAG  y PHASE  y DB y PHN 

CONTINUE 

CONTINUE 

IF ( IPL AN, KO.l)  GO  TO  GO 

CONTINUE 

CONTINUE 

FORMAT  (IX  y "Yl  2-  "«E13.-1y*  <MHO:-’ yf- 2.2y  * ■(DEOy " ySXy  “DB-  dE12.: 
I “y  NORM  PHASt=  'yFZ.Z.  “ <UEG:“) 

FORMAT (/J OX y 

FORMA!  < lOXy  " MUTUAL  ADMITTANCE  OF  SLOTS  ON  A CYl.  J f-ht R'  *) 
FORMArclOXy “ MU  I UAL  AUMTTTANCE  OF  SLOTS  ON  A PLANE  *) 


'''  /»,>•'•.•-«  . , , -rr-'  .•'  " yf^ 


SUBROUT I NE  APROXl < AO » BO » ZO  r YO » PI » ZY2 ) 

IMPLICIT  COMPLEX  <C»H»Z) rREALlA^BiB^Os P-Y) 
real  TrJrio)»TNPi<10) 

REAL  KAfZO 

COMMON  /riATAl/ TN> TNPI » RHO r Cl  »C2 » F2 
COMMON  /CF/CMF  > CtlF  y COl F y GOPF  y GUPF 
THElHA'=ATAN2(ZOyY(U 

IF  ( niETHA.0E*89y<;^9)«PI/i80. ) tHETHA~89y 99*PI/ld0y 
f("=SORT  < Z0*i!2+Y6**2  > 

M-IOO=:F;UO/COS  < THETHA  ) **2 

ZOR=  ( 0 ♦ y -1  i ) JkCEXP  < CMPLX  ( 0 ♦ EO  r ~R ) > / ( 240  ♦ !KP  l!KiK2*R ) 
XA==R5i<ABJU<l./<2y)|:RHaG)K*2))**(ly/3. ) ) 

IF (KAyLT. 0.7)  GO  TO  20 
CALL  FOCK<KA) 

GO  rU  30 
(.ALL  FOCKKKA) 

Z r';l  ■=>!;UF- ^ < S in < THETH a > **2+  < O ♦ y l . ) *COS  < ZJIcTHETHA  ) /R ) 

ZT2-  ( 0 . y .1 . ) *CUF:((C0S  < THETHA ) »*2/R 

Zr3~<0.  yl.  )*<l./(SaRT<2.  )*RHO)|cCOS<THETHA)  ) )**(2./3.  > 

8'  )KSIN(  THETHA  ))«*43((nUPF 
HFHI"ZGR)C  ^ ZTl  (ZI 3 ) 

IF  (IIIETHA.EO.O.  ) GO  TO  40 

I M l-S ] N < BOJKSIN  < ■(  HETHA ) /2 ♦ ) / ( B0*SIN  ( THETHA ) /2 . ) 

(.iO  TO  SO 
CMl'-l. 

rM2  - COS  ( AOjKCOS  ( THETHA ) /2 , ) / < 1 . - ( AO»COS  < THETHA ) /PI ) **2 ) 

ZY2--8 . !FAO*BO*  ( TM,1*TM2/PI  )**2*HPHI 

RETURN 

EMfi 

SCBROUTI Ni;  APROXP ( AO y BO  y ZO y YO  y PI  y 2Y2 ) 

IMPLICI 7 COMPl.EX  (CyHyZ)  yREAL^A-ByD-GyP-Y) 

REAL  7N<  U))yTNPI(10) 

KA /ZO 

CO MMOM  /DA  7 A .1  /TN  y TNP 1 1 RHO  y C 1 y C2  y F2 
COMMON  /Cr/COF  yCUFyCVlFyCVPFyCUPF 
THETHA- ATnN2(Z0yY0) 

Tr( IHETMA.nE.89,99*Pl/180. ) THETHA=e9,99*PI/180. 
R==SQRT(Z0**2»YC*;|!2) 

(oHO(i:--^RHO/COS  < THETHA ) 

ZGR=-<0.  y"1  . )*C:EXP( CMPLX (O.EOy-R)  >/(240.!KPl!X*2*R) 

KA-=  RJCADS  ( ( I . / ( 2 / *RH0G**2 ))»)!:( 1 . /3 .) ) 

’i- (KA.LT.O./)  GO  TO  20 
:'AI.L  FOCKCKA) 

CO  10  30 
(.'Alt  M)CI\.l(KA> 

/ I I -/CVE  * ( COS  I THETHA ) **2-  < 0 . y 1 . ) »COS  ( 2 . *THETHA  > /R  > 

Zr2"<0.  yl,  )*CUF*SIN(THETHA))(oK2/R 
il/^-Z0R)(t(ZTi-FZT2) 

7(11  • (3IN(AO)CCOS(THE7HA)/2.  )/(AO!t!COS( THETHA )/2.  ) )**2 

TM2-=  ( COG  ( BO»S IN < 7 HL THA ) / 2 ) / < 1 . - ( BO*S I N ( THETHA ) /P 1 ) )(:*2  > ) )7:!»c2 

ZY2-=  -8 , »A0;(:B01:7  hl«  rM2*HZ/PlJjc*2 

R(.7URN 


i 


FWii 

■=•••  •'■■■  • ■ 3i/S!f  !/■:;■>: 


Ic^l 


iO 


20 

2r5 

80 

JCO 


200 

:>M0 


END  " 

SUBROUTINE  FHFN(X»N»FM^FN) 
real  DUM1(400)»DUM2(460,X 

REAL  FJ(400)»XBrBSSy’([4o6)»FH(400)»FN<.406)  : 

PI=3, 14159265 
XB^=X 

CF(X-Oa)  10»l0»20-- 
0AML0G=AL05 < X/2 .) +0 . 5772156649 
X2=X*X 
■ X3=X2)KX 
X4=X)KX3 
X5=X!KX4 

BSSY 1 =2  ♦ » ( GAMLOO*  ( 1 , -X2/4 « +X4/64  ♦ ) +X2/4 . ~3  ♦ *X4/1 28 . ) /PI 
BSSY2=-2 . / (pI*X ) +2 . * < GAMLOG* ( X/2 . -X3/16 , +X5/384 . ) ~X/4 . +1  ♦ 25iKX3/16 , 
&-3433333JKX5/768,  )/PI 
GO  TO  25 

CALL  BESY<X»0»BSSY1»IER) 

CALL  BESY<X»1>BSSY2»IER) 

CONTINUE 

BSSY(l)-BSSn 

BSSY(2)«BSSY2 

riBSSYl~--BSSY<2) 

I“1 
I-~I  + 1 

BSSY(I  + l)=2,*(I-l))»!BSSY<I)/X-BSSY(I-l) 

BSSYI1~BSSY(I+1) 

IF<ABS(BSSYI1),GE.1.E10).  GO  TO  100 

GO  TO  80 

NMAX=I+1 

IF(NMAX.GE.N)  NMAX==N 
Nl-N-1 

CALL  BSL J2 < XB , F J r NMAX+ 1 » 0 . DOO 1 7 » lERR » BUM 1 r DUM2 ) 

IiFJl*-FJ(2) 

FN  ( 1 ) ^-=1 , / ( BSSY  ( 1 ) »*2+F J ( I ) **2 ) 

I-  N ( 1)  =■■  I . / < BBSS Y1  !t:*2+DF J1  **2  > 

BO  200  I = 1>N.1 

IF(T.GE.NMAX)  GO  TO  250 

BBSSY=BSSY<  I)-I)KBSSY<  I + l )/X 

BF  J=F  J < I ) - I*F  J ( IT- 1 ) /XB 

FM  < I f 1 ) - 1 . / < BSSY  (111)  )(c*2+F  J < I + t ) !C!(t2  > 

FN(m)  = l./<BBSSY**2+BFJ!(t*2) 

CONTINUE 
CONTINUE 
N---NMAX 
KC [URN 
END 

SUBROUTINE  'BESY' 

PURPOSE 

COMPUTE  THE  Y BESSEL  FUNCTION  FOR  A GIVEN  ARGUMENT  AND  ORDER 


USAGE 

CALL  BESYCX.NfBYf lER) 


89 


77 


m 


3 

'1 


■M 


'I 

I 


* M 


riES(3RIPf.fbN  OF  PARi^MEtERb 

X -THE  ARGUMENT  OF  THE  Y BESSEL.  FUNflTiOH  BESiRED 
H -THE  ORDER  OF  THE  Y BESSEL  FUNC'ti&N  DESiRED 
BY  -THE  RESUUTANt  V DESSEI.  FUNC'noiil 
IER^REfUJI..'!‘AHf  ERROR  CODE  WHERE 
1ER==0  NO  ERRdR 

TER=1  N .tS  NEGAT-tOE  ’ 

1.ER=2  X CS  NEGAtiVE  OR  ZERO 

TER^ri  BY  HAS  EXCEEDED  MAGNITUDE  OF  10)l(*70 


VERY  SMALL  VALUES  OF  X MAY  CAUSE  THE  RANGE  OF  THE  LIBRARY. 
FUNCT.fON  ALOG  TO  BE  EXCEEDED 

X Ml'sr  De  greater  than  zero 

N MUST  BE  GREATER  THAN  OR  EOUAL  TO  ZERO 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  RFOIIIRED 
NONE 

METIiOn 

RECURRENCE  RELATION  AND  POLYNOMIAL  APPROXIMATION  TECHNlOUE 
AS  DESCRIBED  BY  A.  J.M.HnCHCOCK* 'POLYNOMIAL  APPROXIMATIONS 
TO  BESSEL  FUNCnONS  OF  ORDER  ZERO  AND  ONE  AND  TO  RELATED 
rUN'inONS'f  M.T.A.C,.  V.  1.1, 1957 »KP.£)(f.~OBy  AND  G.N.  WATSON > 
A TREATISE  ON  THE  THE'DRY  OF  BESSEL  FUNCTIONS CAMBRIDGE 
UNU'ERSH  T PRESS » ,t9SS.  P,  6'J 


" 


SUBROUT INL  DE  SY  < X » N ? BY  f lER ) 

CUf.'CK  F OR  ERRORS  .IN  N AND  X 

lF(N)I80!'l<)f  10 
to  lli:R“<) 

^|•(X)190^190»20 

BRANCH  IF  X LESS  THAN  OR  EUUAI  4 

20  tr (X -4. 0)40. 40. TO 

COMPUTE  YO  AND  Yl  FOR  X GR'EATER  THAN  4 

30  Il~4.0/X 
r2=r.t:S!ri 

P0-(  ( ( (-•  ,0000037043KT;T{  . 0000 1 73SA5 ) K r2  .00004S7613))FT;! 
I t . 0001  7S‘)Z ) 'M  ('2- . oO  l ) t TS  l . )90y4:;3 

NO'.  < ( U .0  >O00.i2.i?i2!KT?  ■.  00001  4207(0 *121 . 0000342460 ) *1 2 
I .OOOOti.iV/yi  .0004564.J24))t:i;?-  .01246694 

PJ  - ( ( ( I . 0000()4,.>4.1  4*  T2-  . 0000200920)'KT2f.  0000500759)*  12 
I • .00022320..T)1'!21.002921«26)*('21 .3909423 


til's*  HWW>.i«iii»iWPi«riii|«>^^ll>'y«<^|l'^^»W>»^ 


t 


C 

t; 

t: 

c 

r: 

c 


c 

c 

c 


c 

c. 

c 


Ql=( ( < < -.0000036594*72+. 00001 A22)»T2-. 0000398708 )*T2 
1 t + .6001064741 ) ¥72= i 0006396460 )*f 2+ i6374d684' 
A=2.6/SaRT(X) 

B=A*7l 

C=X-i 7853982 

Y0=A*P0*SIN ( C ) +B*Q6*G0S ( C ) 
Yl=-A*PljicC6S(C)+B*Ql*SIN<C) 

GO  70  90 

COMPUtE  Vo  AND  Y1  FOR  X LESS  7HAN  OR  EQUAL  70  4 

40  XX=X/2. 

X2=XX*XX 

7=AL0G(XX)+. 5772157 
SUM=0. 

7ERM=7 

Y0=7 

DO  70  L=l»15 
1F<L-I)50r60f50 
SO  SUM=SUM+1. /FLOAT (L-1) 

60  FL--=L 
7S=T-SUM 

TERM= ( TERM* ( -X2 ) /FL**2 ) * ( 1 . -1 . / ( FL*TS ) ) 

70  YO=YO+TERM 

TERM  = XX*<T-.5) 

SUM=0. 

Y1==TERM 
no  80  L=2»16 
SUM=SUM+1. /FLOAT (L-1) 

FL-L 

FL1=FL-1. 

TS=T-SUM 

TERM=(TERM*(-X2)/<FL1*FL))*( (TS-.5/FL)/(TS+.5/FLl)) 
80  Yl=YtfTERM 
PI  2--.  6366 198 
Y0=PI2*Y0 
Y1---PT2/X+PI2*Y1 

CHECK  IF  ONLY  YO  OR  Y1  IS  DESIRED 

90  IF(N-l)100rl00a30 

RETURN  EITHER  YO  OR  Y1  AS  REQUIRED 

too  IF<N)110»120»110 
110  BY-Yl 

60  TO  170 
120  BY=YO 

GO  TO  170 

PERFORM  RECURRENCE  OPERATIONS  TO  FIND  YN(X) 


-■'4  •; 


."1 


•■I 


79 


91 


-IJl  l^l, . 


«S-^-A»!lr»»'»l«li 


PPBgP 


.i.;30  YA=YO 
YD-Yl 
I<<.=  1 

140  ■I=^FL0AT<2*K)/X 
YC=T*YU'-YA 

,1:F  ( ABS  ( YC ) ~1 . 0E70 ) 145 1. 145  f 14  J 

141  IER=3 
RETURN 

145  K--KM 

IF<K-N)150.  L/jO  1-150 
150  YA=YB 
YB==YC 
CiO  TO  140 
1<S0  BY=Yn 
170  RETURN 
180  TER-1 
RETURN 
190  TER-S 
RETURN 
END 

C OUUROUTINE  riSLJZ(X  , FJ  i-  NMAX  » A r ND  » .TERR  r FJARRX  ? RR> 

0 

C:  THI8  IS  ONE  OF  THREE  ROUTINES^  "BSLJZ'f  "BBLIZ"^  ANP  "BSCJZ'r 
f BABFB  on  algorithm  23A  FROM  “COMMUNICATIONS  OF  THE  A.C.M.“» 

C AUGlIi-^r  I9<f.4,  THIS  ONE  EUALUAIES  THE  BESSEL  FUNCTIONS  OF  THE 
C'FIRSI  KIND  FOR  REAL.  ORIiFRS  ANT  NON-NEGAI  CUE  REAL  ARBUMENIS* 

C 

C THE  F'ARAMETERS  ARE  DIJSCIvMBlH  AS  FOl.LOUSf  WITH  “(I)“»  “(U)"r  AND 
G “(I/O)"  INdlCAVtNOf  RESPECrU'ELY^  THAT  A PARAMETER  IS  TO  BE  SET  ON 
C L'NCRYy  UII.L  BE  SET  BY  THE  ROUIINEr  OR  BOTH  : 


C 

C Y:s;1'  ALL  PARAMETERS  EXCEPT  "NB*  f '1ERR“  > “NMAX*  ARE  *** 
C SINCUE  PRECISION  REAL  NUMHEPS  OR  ARRAYS.  *** 


C 

i.;  ( f ) X 

C ( 0 ) F J 

C 

<•* 

t/ 

c 

c 

C (I ) NMAX 

C ( I ) A 

C 

f' 

C ( 1 ) NB 

C 

C (0)  I ERR 

C 

f. 


C iU)  FJAPRX 


THt.  (NON-NEGATIOE)  ARGUMENT  TO  THE  BESSEL  FUNCTIONS. 
AN  ARRAY  IN  WHICH  THE  UALULS  OF  THE  BESSEL  FUNCTIONS 
ARE  STOREB»  AS  FOLLOWS:  LET  J(X;H>  DENOTE  THE  UALUE 
OF  rUE  BESSEL  FUNCTION  OF  ORDER  B WITH  ARGUMENT  X. 
(HENy  FOR  I - 1 TO  ADS ( NMAX )ily 

FJd)  = J(X:a  •»  d-1  ))t:GlGN(NMAX) ) . 

Rl'Fl.R  TO  "FJ", 

T-;i:i-l:ls  TO  "f  J“,  NORMALLY y 0 <=  A < ly  BUT  THE  ALGOR- 
ITHM WORKS » WITH  SOME  LOSS  OF  ACCURACY r FOR  A >-  1. 
S'Eir  THE  PROGRAM  MOTES  BELOW. 

THJS  GIVES  THF  NUMBER  OF  SI-GNIFICANT  FIGURES  OF 
accuracy  desired  [N  the  function  VALUES, 
rues  IS  AN  ERROR  FLAG  WHJCH  TO  SL'T  TO  0 IF  THE 
INPUT  PARAMETERS  ARE  OKAYy  AND  TO  SOME  POSITIVE 
VALUE  II-  JNE  OF  THE  PARAMETERS  IS  INVALID.  REFER 
TU  HIE  ERROR  EXITS  AT  THE  END  UF  THE  CODE  FOR  A 
Bl  lAXLLD  LIST  OT  HIE  VALUES  OF  H-RR. 

A SCRAILH  ARRAY  USED  BY  THE  EcUUTTNE.  IT  MUST  HAVE 
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80 


I 

i 

i 


c 

c 

c 

c 

v, 

c 


(0) 


RR 


AT  LEAST  ABS<NMAX)+1  ENTRIES. 
ANOTHER  SCRATCH  ARRAY.  IT  TOO  MUST 
ABS(NMAX)41  ENTRIES 


HAVE  AT  LEAST 


OrilER  ROUTINES  CALLED,*  ( * INDICATES  A LOCAL  ROUTINE  ) 


)K 


C 

(.; 

c 

c notes: 
c, 
c 
c 

i: 

c 
c 
c 
c 
c 
c 
c 
c 

(; 
c 
c 
c 

i; 

c: 


c 

c 

c 

c: 

c 


NBSOIZ  INVERSE  FUNCTION  OF  X*LOG(X> 

UNDERZ  — • ROUTINE  TO  CONTROL  UNDERFLOW  INTERRUPTS  ON  THE  IBM  360, 
MOAMMA  --  GAMMA  FUNCTION  FROM  THE  IMSL  LIBRARY. 

AL0(3  “ LOGARITHM 

ADS  — ABSOLUTE  VALUE 

MOD  --  REMAINDER 

AMAX.I  — MAXIMUM  OF  2 REALS 


THE  METHOD  OF  COMPUTATION  IS  A VARIANT  OF  THE  BACKWARD 
RECURRENCE  ALGORITHM  OF  J.C.P, MILLER  (REFERENCE  1),  THE 
PURPORTED  ACCURACY  IS  OBTAINED  BY  A JUDICIOUS  SELECTION 
OF  THE  INITIAL  VALUE  “NU"  OF  THE  RECURSION  INDEX  (REP- 
RESENTED IN  THE  CODE  BY  THE  VARIABLE  "XNU*)»  TOGETHER 
WITH  AT  LEAST  ONE  REPETITION  OF  THE  RECURSION  WITH  "NU* 


REPLACED  BY  "NU‘T5.  NEAR  A ZERO 
FUNCTIONS,  THE  ACCURACY  OF  THAT 


OF  ONE  OF  THE  BESSEL 
PARTICULAR  BESSEL  FUNCTION 


MAY  DETERIORATE  TO  LESS  THAN  “ND"  SIGNIFICANT  DIGITS.  THE 


ALGORITHM 

LARGE. 


IS  MOST  EFFICIENT  WHEN  X IS  SMALL  OR  MODERATELY 


rue,  ABOVE  PAICAORAPH  IS  TAKEN  FROM  GAUTSCHI'S  PRESENTATION 
OF  ALGORITHM  236  IN  C.A.C.M.  THE  SELECTION  OF  THE  INITIAL 
•NU"  IS  DONE  WITH  THE  AID  OF  THE  FMJNCTIDN  NBSOlZ,  ALSO 
BY  GAUTSCHI  (AND  CALLED  *T"  BY  HIM).  IN  THIS  CODE,  THE 
FOLLOWING  SPECIAL  CASES  HAVE  BEEN  ADDED: 

A.  X~0  WHEN  NMAX  > 0 OR  A«(> 

B.  A-O  AND  NMAX  C 0 

C.  A ;:•«  1 : THE  ALGORITHM  WORKS  IN  THIS  CASE,  BUT  THE 

INITIAL  CHOICE  OF  ’Nir  IS  NO  LONGER 
OPTIMAL,  AND  SOME  ACCURACY  IS  LOST,  SIMPLE 
TESTS  INDICATE  THAT  ONLY  A FEW  DECIMAL 
PLACES  ARE  SACRIFICED  AT  WORST.  A LIMIT  OF 
•ABIG*  IS  PLACED  ON  A TO  AVOID  OVERFLOW  IN 
THE  GAMMA  FUNCTION.  TO  AVOID  COMPLICATIONS, 
NMAX  IS  REQUIRED  TO  BE  NON-NEGATIVE  IF  A > 1. 


C 


references: 

1.  GAUTSCHI,  W,  'RECURSIVE  COMPUTATION  OF  SPECIAL 
UNIVFRSIIY  OF  MICHIGAN  ENGINEERlNf  SUMMER 
ENCES,  NUMERfCAt  ANALYSIS,  1963, 


FUNCTIONS' , 
CONFER- 


C ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

SUBROUTINE  BSLJZ(X  , FJ  , NMAX  , A , ND  , lERR  , FJAPRX  , RR) 
REAL  NBSOlZ 

DIMENSIUN  FJ(1)  , FJAPRX(l)  , RR(1) 

LOGICAL  NEVEN  , AFLAG 
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F 

‘I 

i-.- 

I 

t; 

iv 


^ I 


f HALF/.5i:iO/ 
r Cl/. 73576110/ 
» C4/1.3f.)63IiO/ 
f TW0F’5/2.5D0/ 


nOTA 

■ * 

;C 
:f: 

'l: 

. ^ 

C INITFALIZE:  ri-IF  error  PARAMETEF?  » turn  underflow  off  » AND  CHECK 
"C  THE  PARAMETER?)  FOR  UAL  ID  [TY  AND  FOR  THESE  SPECIAL  CAfiES; 

C A.  X'^O  Wni-I  NMAX  > 0 UR  0:^=0 

V.  B.  (1"0  AND  UMAX  < 0 


ONE/IDO/  > TW0/2D0/  ' 

TEN/IOD'O/  » SMALL/lD-15/ 

C2/1.3591D0/  » C3/2.302ADO/ 

ZERO/ODO/  r ADXG/55D0/ 

ALEPH/3777  0000  0000  0000  OOOOB/)'  F0UR/4D0/ 
C5/2000  4000  0000  0000  OOOOB/ 


c 


nil  HLLIBERATELY  AUOIDS  TESTING  MORE  THEN  ONE  THING  IN  EACH 

i.OCiCAL  "IF"  BELOW  BECAUSE  OF  t.B.M.  FORTRAN  INEFFICIENCY  IN  THIS 
REGARD . 


c :i:r-  a>:u 

:;<>.<  :L’i: 

C 


NMAX  MUST  NOT  BE  NEGATIUE. 


I ERR  - 0 

CALL  UND 

[-RZ(  'OFF'  r 

SAUE ) 

IF  (A 

,Lr.  ZERO) 

()0T0 

999 

IF  (A 

.GT,  ABIG) 

GOTO 

990 

1F(X 

• I.T,  ZERO) 

GOTO 

997 

IF' (NMAX 

,0E.  0 ) 

GOTO 

10 

■iF(A 

.E«.  ZERO) 

GOTO 

H) 

fl  (A 

,LE.  FiMALL)  GOTCJ 

i 996 

IF  (A 

.GE.  ONE 

) GOTO 

1 994 

FF  (X 

.OT.  ze:ru) 

GOTO 

40 

IF (NMAX 

, GE  . 0 ) 

GOTO 

20 

IF  (A 

.OT,  ZERO) 

GOTO 

995 

I',' 


(.:  IF"  NMAX  . Or  NMAXT  IS  SET  HERE  SO  TFIAT  ONLY  J(XrA>  IS  CALCULATED, 
..  ?HE  LOOP  FOLLOWING  STATEMP.Nr  BOO  THEN  CALCULATES  IflE  REMAINING 
FUNCTIONS  BY  A SIMPLE  RECURRENCE, 

C,  ir  A-’O*  NMAXT  IS  SET  SO  n-lAI'  J(X»AFN)»  N»0» . . , r -NMAX  ARE 
C CAl.CULATEDr  THE  CODE  AFTER  800  TFIEN  REUERGES  THE  SIGN  OF  EUERY 
C 01  HER  ONE, 

C 

C WE  FIRST  HANDLE  THE  CAFJE  X=0, 

,iC  NTEMP  = lABSi(NMAX)  •(  1 

BO  30  f .!rNTi;MP 

30  F JU)  - ZfJNJ 

IF  (A  .i:0.  ZERO)  FJd)  - ONE 
GOrO  1000 

CL  ;!<>!! 

40  aFLAG  - (A  .F:U.  ZERO)  .AND.  (NMAX  ,LT,  0) 

NMAXT  =■-  NMAX 

ir(NMAX  *Lr.  0)  NMAXr  - 1 
NTEMP  - MAXO<NMAXil»l) 

JK.NOr.  AFLAG)  GOTO  60 
NMAX  I --  • NMAX 
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itj. 


NTEMP  = NMAXT  + 1 
60  EPSLON  = TEN*)|t(-Nn)/2 
DO  80  I = 1»NTEMP 
80  PJAPRXd)  = ZERO 

CALL  MGAMMA(ONE+A  r RESULT  » lER) 

SUM  = (X/TWO))«)«A/RESULT 
rn  = C3!(cND  + CA 
R = ZERO 

1P(MMAXT  .GT,  0)  R = NHAXT  * NBSOlZCHALF^ni/NMAXT) 

S ” C2  * X * NCS01Z(C1*D1/X) 

C THE  RECURSION  INDEX  "Nil"  IS  DELIBERATELY  CALCULATED  AS  A FLOATING 
C POINT  NUMBER  RATHER  THAN  AN  INTEGER^  AND  ALL  COMPARISONS  WITH  IT 
C ARE  DONE  AS  FLOATING  POINT  COMPARISONS. 

)KH(^ 

XNU  --  ONE  (■  AMAXKRrS) 

XLIMIT  =!  XNU/2 
TWOA  “ A f-  A 
XN  - ZERO 
FL  ••■=  ONE 

i; 

L !!li:  OUrER  ITERATION  LOOP  STARTS  HERE. 

C 


C THE  FOLLOWING  l.OOP  IS  DONE  ENTIRELY  IN  FLOATING  POINT  FOR 
C EFFICIENCY, 

200  XN  = XN  ONE 

FL  - FL  * (XN  T A)/(XN  F ONE) 

IF<XN  .1.  r.  XLIMIT)  GOTO  200 
OLDFL  ~ FL 
OLDXN  = XN 

c; 

N - 2*XN 
XN  = N 

NEVEN  ^ . IRIJE. 

R - ZERO 
!!  ZERO 

r ; - TWO/X 

C IN  THE  FOLLOWING  L0(1P»  THE  SUCCESSIVE  VALUES  OF  ’R*  ARE  PARTIAL 
C FRACTIONS  OF  A CONTINUED  FRACTION. 

MO  DENOM  = TEMPI  *(A  + XN)-R 

IF(ABG(DENOM)  .LE.  SMALL)  DENOM  = DENOM  i SMALL 

R = ONE/DENOM 

FLMBDA  = ZERO 

IF (.NOT,  NEVEN)  GOTO  400 

FL  = FL  * (XN  F TUO)/(XN  + TWOA) 

FLMBDA  = FL  K (XN  F A) 

400  S = R * (FLMBDA  F S) 

IF(N  .LE.  NMAXI)  RR(N)  = R 


N -“  N - 1 
XN  •==  XN  - ONE 
NEVEN  = .NOT.  NEVEN 
.1;F(N  .GE.  1)  GOTO  300 

l-JCI. ) =•  GUM/ (ONE  -1-  S) 

i;f(nmaxt  .eq.  o)  goto  aoo 

DO  GOO  N ==  1»NMAXT 
300  FJ(N{-.l)  - RR(N)  * FJ(N) 

(I-IE  LOTEGT  AFFROXIMATIONS  ARE  CHECKED  FOR  IMPROOEMENT: 

<fiOO  lU)  BOO  N - lyNTEMP 

1:F(APB(I  J(N)  --  FJAPRX(N))  .LE.  ABB(FJ(N>  ))HEPBl.ON)  GOTO  BOO 
DC)  700  M lyNTEMI-' 

700  FJAPRX<M)  = FJ(M) 

Xfi  ™ OLDXN 

l-l..  !-  01.I.IFL 

xij:m3  ( “ XL:iM.ir  ■»•  iwopg 

GOTO  :»<)() 

CilNTlNUE 

ircNMAX  .ni:;.  o>  goto  h)oo 

V.  ,![•  NMAX-.;()y  WE  HAOE  FrNTBHEI.i  OBTAINING  J(X?A>  y AND  NOW 
C r»ERAT!:.  TO  F INU  AI.L  THE  BEFSIRED  FUNCTIONS, 
f! 

C FIRSl  WE  CHECK  FOR  THE  SPFCTAL  CASE  A“0. 

IF(,NOT.  AFI-AG)  GOTO  FJGO 
NMAXr  - -NMAX  f 1 
liO  BPO  N 2vNi'lAXTy2 
S20  FJ<M)  - !“.J(N) 
iiOro  1000 

030  I JCi)  TWO  >l<  A * FJ(.L)/X  - FJ(1>) 

. [I•<NMAX  .EO.  - J, ) GOTO  1000 

C ,'!•:£  FOLI..OWING  COBF:  IS  A RENIUTION  OF  THE  LOOP 
0 DO  900  N -■=  2yNMAXT 

C 900  l-.KNll)  <2/X)J|!(A~-N)*FJ(N)  “ FJ(N“1) 

C 

C W,l.  FH  OULRFLOW  UF  FI-CriON,  AS  SOON  AS  IHF  NUMBERS  GET  TOO  BIGy  THEY 
0 .'ikL  SCAlEn  DOWN  (!!'  A POWER  OF  THE  MACHINE  BASEy  SO  .AS  TO  AUOID 
C loss  OF  PRECISION)  AND  THF  CALCULAIION  CONI INUES.  A SEPARATE  LOOP 
L IRANiiFORMO  I'HE  SCALF:D  OALUES  TO  I'llE  CORRECI  OUTPUF  OALUESy  SETTING 
C T(TO-LAR()E  ONES  TO  Pl.DS  OR  MINUS  IN'' IN  IT  Y. 

C K*  IK  :{<■{:  .'l: 

NMAX  I = -NMAX  ■»  I 
r ,INil/  ■■  I .)i  ! ) 

, ;Ni'1  I Ji  •) 
ilMI'1  ■■ 
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XNMl  = TWO 

0 

DO  880  N = 3fNMAXT 

r-JN  = TEMPI  * (A  - XNMl)  * FJNMl  - FJNM2 

1-JNM2  = PJNMl 

FJNMl  = FJN 

FJ<N)  = FJN 

XNMl  = XNMl  + ONE 

Fm(N>  = OVER 

IF(AEIS(FJN)  ,L.T.  C5)  GOTO  880 
OVER  = OVER  + ONE 
FJNMl  “ FJNM1/C;;=I 
FJNM2  = FJNM2/C;S 
880  CONTINUE 
0 

Il'  CNMAXT  .LE.  3)  GOTO  1000 
OVER  = ZERO 
SCALE  = ONE 
C 

DO  900  N - 4.NMAXT 

IF (OVER  .LT.  FOUR)  GOTO  890 

FJ(N)  = 8IGN(ALEPHfFJ(N) ) 

GOTO  900 

090  IF(RR(N)  .GT.  OVER)  SCALE  = SCALE  * C5 
FJ(N)  = FJ(N)  * SCALE 
OVER  :=  RR  < N ) 

900  CONTINUE 
GOIO  1000 

(.« ^ 

C ERROR  EXtlfi  FOLLOW.  MEANINGS  OF  THE  EXIT  VALUES  OF  "lERR"  AREJ 
C 0 i NO  ERROR 


1' 

.1  J A < 0 

r 

2 } A > ABIG 

1 

c 

3 .*  X < 0 

4 ! NMAX  < 0 AND  0 < A < 

SMALL 

C 

S J X«0;  NMAX  < 0»  AND  A 

> 0 

C 

6 S NMAX  < 0 AND  A >=  1 

' 

99^ 

lERR  .=  TERR  i 1 

998 

lERR  ==  lERR  T 1 

{ 

vy6 

[ERR  ■=  lERR  + 1 

1 

99  Z 

lERR  ■ lERR  + 1 

! 

998 

lERR  = lERR  + 1 

! 

999 

JERR  ==  lERR  + 1 

f 

1000 

CONTINUE 

1 

C 

CALL  UNDERZ( 'S' »SAVE) 

1 

\ 

RE  TURN 
END 


REAL  FUNCTION  NBS01Z<X) 

CjK***** 

C IHIS  IS  A NUCLEUS  FOR  FHE  THREE  BESSEL  FUNCTION  ROUTINES 
C "BSLJZ"  f 'BSLIZ*  » "BSCJZ*  BASED  ON  ALGORITHM  236  FROM 
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C ■ COMMUNICATIONS  OF  THE  A.C.M.*, 

C IT  l-:OALllAIL=:S  THE  INVERSE  FUNCTION  OF  X*LOG(X)  FOR  X >-  1 TO  AN 
•C  ACCURACY  OF  ABOUT  ONE  PER  CENT. 

C t OR  THE  INTERVAL  0 <=  X <=  10  A FIFTH  DEGREE  APPROXIMAtlON  IS 
C USED.  OBIAINED  BY  TRUNCATING  AN  EXPANSION  IN  CHEBYCHEV  POLYNOMIALS, 
i;  FOR  X J 10 f A DIFFERENT  APPROXIMATION  IS  GIVEN ^ AS  CAN  BE  SEEN. 

DATA  Cl/. 00005794100/  » C2/-.0017A148D0/  r 


Cl/. 00005794 IDO/ 
i:T/.0200.445DO/ 
C5/.85777D0/ 
ALF-HA/.  775110/ 


C2/-.0017A148D0/ 
C4/-  ..12901 3D0/ 
C6/1.10125D0/ 
IEN/.10B0/ 


'FIX  .GT.  TEN)  GO  TO  10 

NoSOIX  = ((((CIJUX  C2>!)(X  1 C3)*X  + C4))l!X  T C5)!l:X  + C6 
RETURN 

TEMPI  = ALOG<X)-ALPHA 

71:MP2  = (ALPHA-ALOGITEMPl)  )/(HTEMPl) 

NDSOIX  = X/(  (liTEMP2)!KTEMPl) 

RETURN 

END 

SUBROUTINE  PROOl ( RHO r AMPY » PHASEY » AMPYDB » PHASNM ) 

PROGRAM  TO  COMPUTE  THE  MUTUAL  ADMITTANCE  BETWEEN  TWO  IDENTICAL 
AXIAL  SLOTS  ON  A CYLINDER  < UI  MODAL  SOLUTION) 

REAL  KO.NZ.KTvISrKZKIRO 
COMPLEX  IlfYi;>fPSIEXP»YN12 

REAL  PI  (400).  F'M  < 400  ).FN(  400  ).AIMAG.  REAL.  ATAN2 
COMMON/OA'I  A5/K0 . NC  YCLE . PHIO . ZO . Yl  1 . MMAX . A . B 
INPUT  PARAMETERS} 

KO.--WAVF  NUMBER  IN  FREE  SPACE  IN  TERMS  OF  1/INCH 

NCYCLE-NO.  OF  SUBSECTIONS  BETWEEN  ANY  TWO  SUCCESSIVE  ZEROS  OF  INTEGRAND 
IN  TRAPEZOIDAL  RULE.  FOR  NUMERICAL  INTEGRATION 
A!«B*-  SLOT  DIMENSION  B>A  <1NCH> 

riiu-radius  of  cylinder  <INCH> 

PI-nO^ANGULAR  SKl’ARAdON  OF  THE  SLOTS  (CENTER  TO  CENTER)  <RADIAN> 

Z0='  SEPARATION  OF  THE  SLOTS  IN  Z-DIRECTION  ^INCH.> 

Yll"^  NORMALIZATION  FACTOR 

MMAXm  MAXIMUM  NO.  01  TERMS  WHICH  HAS  BEEN  USED  IN  CALCULATION  OF 
INFINITE  SERIES 
PI=3. 14159265 
Y0=-1./(120.*PI) 

I RLU-.'i . E 1 0)KK0/  ( 2 . !t!PI  . 54 ) 

ARA=NO)i<A 

BKB-KOiKB 

RK-KO^RHO 

i'H3A'=HALI'  ANGULAR  WIDTH  OF  THE  SLOT 
PH  r •',••'2 . *AS1  N ( A/  ( 2 . *RHO ) ) 

COf-'-UTATION  OF  INPlNTIt  SERIES 

MMA.X12-MMAX11 

BO  .'00  M=1.MMAX12 

,'!1--M-  I 

IT(M.LO.t)  GO  TO  99 

1-1  (M>-(  0K(Ml*PHI0))|f(SIN(Ml)tcPHIA/2.  )/(MlH<PHIA/2. ) ))KJ|c2 
GO  TO  100 
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"iTTiiTv/ri  ifiii'rr  T'1'  fjryifeirirff 


99  F1<M)=0.5 
100  CONTINUE 

INTEGRATION  OF  PSI(KZ)#R1  <M»KZ))|!EXP(-J*KZ#ZO)  BETWEEN  0 AND  KO 
DELTA=  NEIGHBOURHOOD  OF  THE  SINGULAR  POINT  KZ=KO  IN  WHICH  THE  INtEGRAL 
HAS  BEEN  CALCULATED  ANALYTICALLY 
liELTA=t,E-7)t!K0 

NSECT1=N0.  OF  SAMPLES  IN  THE  INTERVAL  (0. f KO-DELTA) ♦ 

NSECT 1= ( i F I X ( ( B+20+RH0 ) JKKO/PI ) +2  ) iKNCYCLE 
DKZ ^ ( KO-DELTA ) /NSECTl 
NSEC r=NSECIl+l 
ll  = (0,»().) 

I1-F[R§T  INTEGRAL  (BETWEEN  0.  AND  KO) 

DO  200  1=1 F NSECT 
1\Z=<I-1)*DKZ 

IF<KZ.ER,0)  KZ=0.00001*KO 
CIN=1. 

IF<(I.E«.l)  .OR,  (I.EQ. NSECT))  CIN=0.!5 
K T=SQRT  ( KOjKKO-KZjKKZ  ) 

TF ( ABS  ( KZ!KB/2 . -PI/2 . ) , LE . 1 , E-8 ) KZ=1 , 000001*KZ 
PSIEXP=<C0S(KZ*B/2.  )/<  (KZ*B/2.  )»*2-(PI/2.  )))t*2)  )!(t*2*CIN)KDKZ 
X*CEXP(<0. F-J . )*KZ*ZO) 

MMAX1=MMAX 

ROKT=RHO*KT 

COMPUTATION  OF  FM<N)  = 1 ./( JN(X))tt*2+YN(X)*)(t2)  AND  FN(N)=1  ./(DJN<X)!|c»2+ 
DYN(X))t!)tt2)  FOR  X=RUKT  AND  N=0  TO  MHAXl  i WHERE  MMAXl  IS  A NUMBER  AFTER 
WHICH  THE  CONTRIBUTIONS  OF  FM(N)  AND  FN(N>  TO  THE  INFINITE  SUM 
BECOME  NEGLIGIBLE.  MMAXl  IS  A FUNCTION  OF  THE  ARGUMENT  X AND  IS  ALWAYS 
LESS  THAN  OR  EQUAL  TO  MMAX.  MMAXl f FM(N>  AND  FN(N)  ARE  CALCULATED 
BY  SUBROUTINE  FMFN ( X , MMAXl fFMfFN) . 

CALL  FMFN  < ROKT  f MMAXl f FM  f FN ) 

DO  200  M=1fMMAX1 
M1=M-1 

200  11  = 11  f FN ( M ) )KPSI EXPJdFl  ( M ) 

COMPUTAriON  OF  12  (BETWEEN  ZERO  AND  ETAMAX  ) WHERE  ETAMAX  IS  A NUMBER 
AFTER  WHICH  THE  INTEGRAND  BECOMES  VERY  SMALL) 

12=0. 

ETAMAX=14./(20~B) 

THE  INTEGRATION  IS  CARRIED  OUT  BY  TRAPEZOIDAL  RULE.  AT  FIRST  THE  WHOLE 
RANGE  OF  INTEGRATION  (O.fETAMAX)  IS  DEVIDEU  INTO  TWO  SUBINTERVALS  *. 
(O.fETAD  and  (ETAIfETAMAX)  f WHERE  ETAl=ETAMAX/2. . THEN  THE  NUMERICAL 
COMPUTATION  OF  THE  INTEGRAL  IS  PERFORMED  IN  THESE  SUBINTERVALS  WITH  THE 
NO.  OF  SAMPLES  IN  THE  FIRST  SUBINTERVAL  TWO  TIMES  THAT  IN  THE  SECOND  ONE. 
ETAl=7./(20-B) 

NSECTl  = (IFIX(SQRT(K0*K0+ETAl)((ite2))HRH0/PI)+2)*NCYCLE 

BETA 1=ETA1 /NSECT 1 

DETA2=2.!)(DETA1 

NSErT2=lF IX( (ETAMAX-ETAl )/DETA2)+l 

NSECT=NSECTl+NSECT2+2 

DO  300  1=1 F NSECT 

IFd.LE.NSECIl  + l)  GO  TO  220 

El  A=ETA1  f ( I-NSECTl -2 ) JUDE  TA2 

DETA=I1ETA2 
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l30  ro  1’40 

220  r£TA=--<I~J  ):t:DETAl 

11' (ETA.ER.O. ) ETA=0.000t/A 
nETA.=riETAl 
240  CtN-l. 

Cr  ( (T.EQ.l)  .(JR.  (I.EQ»NSECTH1)  .OR.  (I.EC).NSECTlf2)  .OR,  <I.En.NSECT) 

%)  CtN=04  5 

PSEX=  ( COSH  ( ETA:KB/2  . ) / ( ( ET A>CEi/2 . ) *i»:23  ( PI/2 . ) !K^:2  > ) )(!)K2)(:nETA!Kf;iN 
!i,»EXP(-ETA*ZO) 

KT=SQRr(KOJ<KO-(  ErA)K)K2) 

MMAX1=M?^.AX 

CALL  I- MEN ( RHOiKK  r r MMAXl  t FM  r F N ) 

1)0  300  M=l>' MMAXl 
ht-M-1 

300  r2-irMFN(M))KPSEX:)!Fl(M) 

Y12=(iH  (0,  .1  . ).tI2))KA!KB:«YO/(PI)KKO*RHO)K*2> 

C NORMALIZATION  OF  THE  PHASE  OF  Y12 

VNl2= Y 1 2*CEXP  < ( 0 . > 1 . ) * < KO*SORT  < ZO*ZO+  ( RHO!t:PHI 0 > **2 ) ) ) 

C COMPUIATION  OF  HIE  ACTUAL  PHASE  'PHASEY'  AND  NORMALIZED  PHASE  'PHASNM' 

C OF  Y12. 

PHASEY"ATAN2<AIMAG<Y12)  rREAL(Y12)  ))(tl80./PI 
PHASNM'^A I AN2  ( A T MAO  ( YN  J 2 ) » REAL  ( YN12 ) ) *100 . /PI 
r COMPUTATION  OF  THE  MAGNITUDE  OF  THE  Y12  IN  TERMS  OF  <MHO>  AND  <DB>. 
AMPY=CABS(Y12) 

AMPYDEi=AI..OG  J 0 ( AMPY/ABSi  ( Y1 1 ) ) *20 . 

RPHIK=RO*RHO*PHIl) 

;:ni;-=Ko*?o 

T'H  [OD«PHiO*lHO./P.I 

T-'E  TURN 

LND 

OUI'ROUI  INL*  l••R0G2  ( RHO » AMP Y y PHASEY  y AMPYDD  r PHASNM  > 

C PROGRAM  FOR  COMPUTATION  OF  THE  MUTUAL  ADMIVTANCE  DETWEEN  TWO 

C IDENTICAL  CIRCUMFERENTIAL  SLOTS  ON  A CYLINDERdJl  MODAL  SOLUTION) 

REAL  KUyKZ.KTyl2yKZKlR0 
COMPl.  r X 1 1 y Y 1 2 y PS  1 EXP  y YN 1 2 

REAL  1-1(400)7  T M(400)  »f'N(400)»AlMAGyREAl  fATAN2 
COMMON/DA  TA3/N0  y NC YCLE  y PHI U y 20  y Y 1 1 y MMAX  y A y D 
C INPUT  PARAMETERS  i 

C KO-LIAOE  NUMDE.R  IN  FREE  SPACE  IN  TERMS  OF  l/TNCH 
C RHO'-RADIUS  or  CYLINDER  ' TNCH..' 

C PHIO=ANGULAR  SEPARATION  OF  THE  SLUTS  (CENTER  TO  CENTER)  -^^RADIAN.)- 

C 70-  SEPARATION  OF  THE  SLOTS  IN  Z-DIRF.CTION  ■nNCH.'; 

C Yll=^  NORMALIZATION  FACTOR 

C MMAX=  MAXIMUM  NO.  OF  ILRM!)  WHICH  HAS  DEEN  USED  IN  CALCULATION  OF 
C INFINITE  SERIFS 

C rCYCLIj^NU.  OF  SUBSEC  I CONS'  BETWEEN  ANY  TWO  SUCCESSIOE  ZEROS  OF  INTEGRAND 

C IN  TRAPEZOIDAL  ROLF  FOR  NUMERICAL  INTEGRATION 

F'I--3. 1415927,5 
YO'H  ./(120.*P1  ) 

FRLQ  =■3 . T J ()*K0/  < 2 . *P  J *2 . 54  ) 

ALA-  KO*A 
T'LP'^MJ*D 


100 


88 


i 


I 


RK=KO)t!RHO  I 

C PH1R=HALF  ANGULAR  WIDTH  OF  THE  SLOT  I 

.■■'HIB=ASIN(A/(2.)KRH0)  ) I 

C COMPUTATION  OF  INFINITE  SERIES 
MMAX12=MMAX-»1 
DO  100  M--lrMMAX12 
Ml=M-t 
EPM"1. 

;tF(M.li:0.1)  EPM=2. 
r•HTD^=PHIB 

[F(ABS(PH[B!KMl-PI/2, ) .LE.l»E-7)  PHIBl=PHIBJ(tl . 001 
too  FJ  (M)=COS(Ml)ltPHIO))(((-PI*C;OS<Ml)|tPHIBl>/((  (M1*RHIB1)!|c*2-(PI/2.  )*)H2 
X)))**2)K<1,/EPM) 

C INTEGRATION  OF  PSI  (KZ)*R1  (M»KZ>*EXP(-J)KKZ)K20>  BETWEEN  0 AND  KO  • 

C DELTA=  NEIGHBOURHOOD  OF  THE  SINGULAR  POINT  KZ=KO  IN  WHICH  THE  INTEGRAL 
C HAS  BEEN  CALCULATED  ANALYTICALLY 
DELTA=0.0001*K0 

C OELFAls  NEIGHBOURHOOD  OF  THE  SINGULAR  POINT  KZ=KO  WHERE  THE  INTEGRAND 
C VARIES  RAPIDLY  AND  'NDELTA'  SAMPLES  HAVE  BEEN  USED. 

DELTAl=0.01))!Kn 

NDELI'A^lOO 

DKZ2---  ( DELTA  1 -DELTA ) /NDELTA 

C NSECT1=  NO.  OF  SUBSECTIONS  BETWEEN  0 AND  KO-DELTAl 
NSECTl  = < ] FI  X ( ( B»  ZOTRMO ) *KO/PI ) +2 ) *NCYCLE 
DKZ 1= < KO-DELTAl) /NSECTl 
NSECT=NSECri+NliELTA+2 

C n=FIRiir  INTEGRAL  (BETWEEN  0.  AND  K0>  [ 

I1=<0. rO. > 

DO  200  Isl»NSECT 
IF<  t.LE.NSECTl  + l)  GO  TO  1.20 
KZ=KO  -DEL.  FA  1 + < I-NSECTl  -2 ) *DKZ2 
DKZ=DKZ2 
GO  TO  140 
120  KZ-(I-l)*DKZl 

IF(KZ.EO.O)  KZ=0.00001!f(K0 
DKZ=DK2t 
140  CSN^l, 

IF((J .EO.l) .OR. (I.Ea.NSECT,l+l).0R.(I.Ea.NSECTl+2).0R. (I.EQ.NSECT) 

X)  CIN=0.5 

KT-SORT  ( KO^KO-KZJKKZ ) 

PSIEXP-(SIN<KZ*B/2.  )/(KZ*B/2. ) ))K)K2*CIN*DKZ*CEXP(  (0.  f-1 . )*KZ*ZO) 

MMAX1=MMAX 

ROKT=RHO*KT 

C COMPUTATION  OF  FM(N)=1 ./( JN<X))K)(!2+YN(X))(t!((2)  AND  FN(N)=1  ./(DJN(X)!)t)K2I  i 

C DYN(X))K)(!2)  for  X=R0KF  and  N~0  to  MMAXl  f WHERE  MMAXl  IS  A NUMBER  AFTER 

C WHICH  THE  CONTRIBUTIONS  OF  FM(N>  AND  FN(N>  TO  THE  INFINITE  SUM 

C BECOME  NEGLIGIBLE.  MMAXl  IS  A FUNCTION  OF  THE  ARGUMENT  X AND  IS  ALWAYS 

C LESS  THAN  OR  EOUAL  TO  MMAX.  MMAXl » FM(N)  AND  FN(N)  ARE  CALCULATED 
C BY  SUBRDUFINE  FMFN( X .MMAXl .FMrFN) . 

CALL  FMFN(ROKF. MMAXl .FM.FN) 

KZKTRO=  ( KZ/  ( K F*KO*RHO ) ) )t!*2 
DO  200  M=l. MMAXl 
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k .1=  ( I . /KT**.?  > * ( F M < M ) +M1  Jl<)K2)KKZKTRn*FN  ( M ) ) 

200  .[l~IllRi*fSIE'XP)KI~.t<M) 

1 1==  ( 2 , JCKU/  < Pi)i;kH0 ) ) * ( n--F.i  a ) jkcexp  < ( o . » -i . > )KKO!«zo  > * ( p,i*pi/  ( 2 . jkko  ) > 

X;K ( GIN  ( K0)icB/2 . ) / ( K0*Ei/2 . ) ) ttl/ ( 2 . * ( 0 . f5772.l56649  FALDG ( RHa)>:S«RT 
8 ( KO/2 , ) ) ) -FALUM  < DELTA ) ) ) 

L L'OMPUTATION  OF  12  (BETWEEN  ZERO  AND  E FAMAX  f WHERE  ETAMAX  IG  A NUMBER 
(.'  AFTER  WHIEH  THE  INTEGRAND  BECOMES  MERY  SMALL) 
i:.'‘-0, 

i;TAMAX-14,/<Z0~B) 

C THE  INIEGRAri'ON  3 G CARRIED  OUT  BY  TRAPEZOIDAL  RULE*  AT  FIRST  THE  WHOLE 
C RANGE  OF  INTEGRATION  (O.rETAMAX)  IS  DEMIDED  INTO  TWO  SUfJINTERMALS  J 
C (O.vFTAI)  and  (EIAi.ETAMAX)  » WHERE  ETAl=ETAMAX/2. . THEN  THE  NUMERICAL 
C COMPillATlON  OF  THE  INTEGRAL  IS  PERFORMED  IN  THESE  SUBJNTERMALS  WITH  THE 

C NO,  OF  SAMPLES  IN  THL  FIRS’)  SUBTNTERMAL  TWO  TIMES  THA'I  IN  THE  SECOND  ONE. 

ErAl=7./(7.0-B) 

nSF  C 1 1 - < J F ! X < SORT  ( KO!«KU  ) f;TAl*!)t2  )*RH0/P1 ) +2)  ^XNCYCLF 

DErAl~ETA.I./NSE(:ll 

DlL’iA’2-2.>^liETAl 

.')!>■!:  iV!  2 - I F ;r  X ( ( £ lAMAX-ETA  t > /DE » A2 ) M 
NGEC'n-N-iEC  I I t NSEC  (2+2 
DO  300  I~l,N;i£i:i' 

I!  (,  .I,n  .NSFCri  + l)  HO  TO  220 
ETA»E  CAT  )■  ( r -MSEC  1 3 -2 ) !i<DE'I  A2 
DErA-DEfA2 
HO  10  'J'lO 

22«  L’IO-.(,l-l):t:nETAt 

JF ;E lA.LO.O. ) ETA-O.OOOl/B 
DETAi-liFfAt 
?40  (.:tN=-l, 

IF  < ( ’.£0. 1) .OR. (1 .EQ.NSECTIID.OR. (I.E0.NSECTI  I2) .OR. ( I .EO.NSECT) 

8)  C1N=0,5 

PSEX  = ( G 3 NH < L-.TA*B/2 . ) / ( ETA*B/2 . ) ) ’l:)C2)KEXP  ( -ETA* ZO ) *DE  I A*0 IN 
KT-SORT ( KO*KU+E  rA**2 ) 

ETKIRO=(EIA/<RHO*KO*M ) )**2 
MMAXI=MMAX 

CALL  FMFN ( RHO*KT  r MMAX 1 r FM , FN ) 

D!J  .300  N~  I,  MMAX  I 
rtl=’M-l 

R1-(I/(K r*KT ) )*<FM<M)-Ml*Ml*FrKIRO*FN(M) ) 

.300  l2-=12»Ft(M)*Kl*PGEX 
J2- l2*2.*K0/<ri*RH0) 

Y12-  (T1  ( (0.  ,1 . ;-*12)*A*B*YO/(2.*Pr*PI*RHO) 

(;  NORMALIZATION  OF  THE  PHASE  OF  Y12 

YNI 2- Y 1 2*CEXP ( < 0 . r 1 . ) * ( KO*SOR  3 ( ZO*ZO + ( RH0*PH1 0 ) **2 ) ) ) 
n COMINJIAtlON  OF  THE  ACTUAL  PHASE  'PHASEY'  AND  NORMALIZED  PHASE  'PHASNM' 

C OF  Y12. 

PHASE Y-A  TAN2  < A t MAH ( Y 1 2 ) , REAL  < r 1 2 ) ) *180 . /PI 
PHitG.NM--AtAN2(AIMAlj(YNI2)  ,REAL(YN12)  )*180./.») 

C COMPUTATION  OF  THE  MAGNirUDE  Of  THE  Y12  IN  lERMS  OF  ’.MHO.-  AND  <DB>. 
AMPY=CABG<Y12) 

AMPVDB-ALOG I 0 ( AMPY/ABS ( Y 1 1 ) ) *20 . 

102 

90 


j,'piMJi,f«lWl|^M'  '■' W-'riJiWPW!W!^|W 


BEST  AVAIIABLE  COPY 

RPHIK=KC)*RHO)|fpHIO 

zoK=k:o*zo 

PHI0D=PHI0)lcl80./PI 

RETURN 

ENn 

E 

C IH.IS  SUBROUTINE  IS  USED  TO  CALCULATE  THE  FUNCTIONS  CVF»CUF»CUlFf CVPFrCUPF 
C 

SUBROUTINE  FOCK(X) 

IMPLICIT  REAL<A-B»D--H,P'-Y)» COMPLEX  <C.Z) 

REAL  TN(.tO)  »TNPI(10) 

COMMON/CF/CVF » CUF  > CVIF » CUPF » CUPF 

COMMON/BA  rA.l/TN»TNPI>RHn,Cl»C2rF2»I0P»CC»RABNrDEG 

Fl=SaRT(X) 

l-.7,-XiK*1.5 

CUF:=0 . 

CUF---0. 

cyiF=o. 

(;upF=o, 

CUF|:=(). 

BO  20  N=l»tO 
ZTN=TN(N)^C1 
ZTNPI  = TNPI(N)!KC1 

. C3»CEXP  < CMPLX  < 0 . 0 » -X ) itiZTNP  I ) 

C4=CEXP  < CMPLX  < 0 . 0 > -X ) #ZTN ) 

CVF=CVF»C;J/ZTNPI 

CUF=CUF+C4 

CU1F=CU1F+C3 

CUPF=<1.0-CMPLX<0.0>2)KX)*ZTNPI))tcC3/ZTNPI+CUPF 
CUPF=  < 1 . 0-CMPLX  ( 0 . 0 f 2 . )KX/3 . ) *ZTN ) *C4iCUPF 
20  CONTINUE 

CVF=F2*FI*CUF/C2 
CUF=2.)KF2)((F3*C2)BCUF 
CV1F=-.2.*F2*F3#C2*CU1F 
CVPF=F2)t!CVPF/  < 2 . »F1  #C2 ) 

CUPF=3 . *F2*Fl*C2))eCUPF 

RETURN 

ENTRY  FOCKl 

XTHREE=X)K)|t3 

F1=SQR TCX) 

F3=X THREE 
Zl=F2)|(C2)|iSaRT(F3) 

Z2=CMPLX ( 0 , 0 » 1 . 0/60 . ) #XTHREE 

23=F2*X*)lc4 . 5/  ( C2Jt(64 . ) ' 

F4=F3*)t:2 

CVF=  1.0-21/4,  +7*22+7 . *Z3/8 . -4 . 141E-3!KF4 
CUF- 1 . 0-2 1 /2 . +25 . *Z2+5 . *23-3 . 701E-2*F4 
CV1F=1 . 0+Zl  /2 . -35 , *72-7 . *Z3+4 . fi55E-2*F4 
CVPF= . 375*Fl*T-2/CC+21 . *Z2/X+63 . *Z3/ ( 1 6 . *X  > -2 . 485E-2*F4/X 
CUPF= . 75*F 1 *F2/CC  fCMPLX  < 0 , 0 » 1 . 25*X**2 ) +22 . 5*Z3/X 
1 -2.221E-1*F4/X 
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FiETURN 

END 
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SlUBROUTilNE  PLANAR  < IJ f ZSUM ) 

L MPL,  I C ;t  T CUMPLEX  < C » F^ » Z ) » REAL  ( A-B » D-0 » P- Y ) 

REAL  ZO 

COMMON  p;i  »TZ:L»TZ2yTY;L  »TY2»R>THETHA 
G0HM0N/DATA2/A  y B r ZO » YO 
(SO  TO  (:l.0y20)  ylJ 
;iO  XM:L^=(  (TZ2-"rZ.t)/R)!t:*2 
XM2==--2,~3.)KXM1 

xrl==xm:vr 

x:i  M:-"XMlTXM2/R)i()K2 

HA==CEXP  ( CMPLX < 0 ♦ EO » -R ) ) *CMPLX ( XRL y -XIM ) / ( 240 ♦ *R!lcPI##2 ) 
FA(:;T0R==C0S  ( PI*TYI  /a  ) JKCOS  < PI»  < TY2-Y0 ) /A ) 
ZSUM=ZS)UMiF"ACTOR*F(A 
k'ETURN 

20  ZA1"CMPLX  < J. . /R)K!K2  y 1 . /R ) * < 2 ♦ -3 . *AC02 ) 

ZA2 -CEXP  < CMPLX  ( 0 ♦ EO  y ~R ) ) * ( AC(J2+ZA1 ) /R 
FiA:»(0,  y ■1  . ))KZA2/ (240* *P  1**2) 

F^-'iCTOR"T;OB  ( P.l*TZ:l/B ) *COS  ( pi*  ( TZ2--Z0  > /B ) 
ZSlJM"-Z!iUM'fF'ACTOR*HA 

rl-jurn 

l;.ND 


THJ:i  CUBROUriNE  .U)  UiSED  TO  GET  THE  "CYLINDRICAL*  SOLUTION 

SUBROUTINE  CYLIND ( IJ y ZSUM ) 

.IMPL  ICIT  COMPLEXCCyllyZ)  yREAL(A-ByD“(3yP-Y) 

REAL  ZOvKA 

REAL  TN<10)  y TNPKlO) 

COMMON  PI yTZI yTZ2yTYl yTY2yRyTHETHA 
COMMON/CF /COI-  y CUE  y COIF y COPF y CUPF 

COMMON /DA lAl /TN  y TNPI y RHO  y Cl y C2  y F2  y I OP  y CC  y RADN  y DEG 
C0MM0N/DATA2/A  y D y ZO  y YO 

CUMHON/DA  CA 4 /AC02  y SN2  y TN2  y R2  y ACONl y AC0N2  y COl y C 1 0 
ZGR  = ( 0 . y - 1 . ) *CEXP  < CMPLX  < 0 ♦ EO  y -R  > ) / ( 240 . *R*PI **2 ) 

IT-^O 

AN(il.F-=  ArAN2(ABS(  rZ2-TZl ) y ABS<  TY2--TY1  > )*DEG 
IF  (ANCl  EJ  T.89.y9)  GO  TO  10 

ri-ii  THA  =-PU-o‘?.yy/ieo. 

ZW  :JO:l/R  FCEXP  < CMPLX  ( 0 . EO  y - PI/4 . EO ) ) *SORT  < PI*R/2 . ) /RHO 
j 0 RH0G“RH0/AC02 

K '4-  f;* ABS  ( (1  * '■  ( 2 ♦ *kH0G**2 ) ) **  ( 1 , /3 . ) ) 

IF  <!.A.LI\0./>  GO  TO  20 
CALL  FOCKU'.A) 

00  )C  30 

20  CALL  FCCKKKA) 
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BEST  AVAIUBIE  COPY 


30  IF(IT.En.l)  GO  TO  40 

ZU-C01)K  ( 1 , /AC02)  /R*  ( CUF-CV1F)KSN2 ) 

40  GO  TO  (21»22f23»24) rlOP 

21  ZHBl-CMPLX  ( 1 , EO » -1 . /R ) JtiCVF 
ZHB2=CIJF/R!K*2 

ZH»3=C01  / ( SORT  ( 2 , EO ) )KRHO  > )((!(!AC0N2 
ZHB4=AC02)K*AC0N2)KC0PF+SN2)K  < CUPF/AC02*#  ( 1 . /3 , ) ) 

HB=ZGR!k  ( 2HB1-ZHB2+ZHB3#ZHB4  > 

ZTM1=C01*CUPF/  ( SORT  ( 2 , EO ) iKRHOG ) )|t*AC0N2 

HT=C01)KZGR/R)K<C0F+CMPLX(l.E0r-2,/R>)lcCUF+ZTMl> 

1F(IJ*LT.2)G0T0  25 

HZ=HB3KAC02f-HTit!SN2 

GOTO  500 

25  HPHI=HB»SN2+Hr*AC02 
GO  ’to  500 

22  HB--ZGR)t!CVF 
l-ir=t;01)KZGR)KCIJF/R 
.lF<lJ.LT.2)G0Ta  26 
HZ=HB*AC02+HT*SN2 
GOTO  500 

26  HPHI=HB*SN2+HTJKAC02 
GO  TO  500 

23  ZTM2=C0lJK<l,-3,)(!SN2)/R 
IF<IJ»LT.2)G0Ta  27 

HZ-ZGRiKCOF*  ( AC02+C01)t!  < 2.-3.  »AC02 ) /R  > 

GOTO  500 

27  HPHI=ZGR*  < CUF*  < SN2+2TM2 ) +ZW ) 

GOTO  500 

24  IF<IJ.LT.2)G0T0  28 

HZA=CVF*  ( AC02*C  1 0+  ( 2 . -3 . *AC02  > K<CMPl.X  ( 1 . /R2  r 1 . /R  > > 
HZB=C01*<CyPF/(SQRT(2.  )J|!RH0G))|t»AC0N2> 

HZC-CMPLX  ( - 11 . / 1 2 . *AC:02 » < -1 1 . /6 . -2 . /3 . *TN2+AC0N1  iHACOZ ) /R ) 
HZ=ZGR*  ( HZA+HZBiKHZC ) 

GOTOSOO 

28  HPA=CVF*CMPLX ( SN2+ (2.-3. *SN2  > /R2 » ( 2 . -3 . *SN2  > /R > 

HPa=C01)|!  < CUF-CVF ) /AC02/R 

HPC=C01  / ( SORT  ( 2 . ) )|!RH0G!|<AC02  ) )|t!|iAC0N2 
HPl D= ( CVPF/AC02** ( 1 . /3 . ) ) 

HP2D=CMPLX ( 4 . /3 . *SN2- 1 1 . /12 . *SN2*AC02  r 
$ .75*AC02-7./12.J|iSN2+AC0Nl»SN2*AC02) 

HP3D=C01/12 . /R»CUPF/AC02**  < 1 . /3 . ) 

HP4D=HPC*  < HP 1 D*HP2IHHP3D  > 

HPHI=Z0R*(HPA+HPQ+HP4D) 

500  ZGREEN=HPHI 

IF(IJ.EQ.2)ZGREEN=H2 

FAC  TOR=COS  < PIiKT  Y1  /A ) *003  < PI  * ( TY2-Y0 ) /A  > 

IF(IJ.EQ.2)FACT0R=C0S<PI*TZ1/B)»C0S(PI»(TZ2-Z0>/B> 

ZSUM=ZSUM+FACTOR*ZGREEN 

RETURN 

END 
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1.  INTRODUCTION 

, t 

This  paper  contains  two  results  for  the  mutual  admittance  Yj^2  between 
two  slots  on  the  surface  of  a large  conducting  cylinder  (including  the 
conducting  plane  as  a special  case).  The  first  and  the  main  result  is 
that  an  approximate,  closed-form  solution  of  Yj^2  derived.  This  solution 
may  be  considered  as  a simplified  version  of  the  asymptotic  solution  of  Y^2 
reported  in  [1],  as  the  two  surface  integrals  over  the  apertures  of  the 
slots  are  no  longer  needed  in  the  present  approximate  solution.  Our 
second  result  concerns  the  derivation  of  an  exact  solution  of  Yj^2»  which  is 
given  in  terms  of  an  inverse  Fourier  transform  and  an  infinite  summation  of 
cylindrical  modes.  This  solution  is  based  on  the  original  expression  for 
Yj^2  described  by  Stewart,  Golden,  and  Pridmore-Brown  [2],  (3),  and  is  more 
suitable  for  numerical  calculation  for  some  cases.  , 

This  work  is  undertaken  for  the  following  reasons.  The  determination 
of  Yj^2  problem  for  2^2  between  two  dipoles)  is  not  only  a 

classical  problem  in  electromagnetics  that  has  attracted  wide  attention 
[1]  - [10],  but  also  an  integral  part  in  the  design  of  modern  conformal 
arrays  [11]  - [15].  In  the  latter  application,  Yj^2  be  repeatedly 

calculated  for  a large  number  of  times.  Thus,  a simple  closed-form  | 

• 1 

solution  should  greatly  reduce  the  computation  effort  and,  furthermore, 
provide  a better  physical  insight  for  the  design  problem  as  the  "cause" 
and  "effect"  can  be  readily  identified  in  a closed-form  solution. 

The  organization  of  this  paper  is  as  follows.  In  Section  2,  we  first 

i 

define  Yj^2»  then  give  the  final  form  of  its  approximate  solution.  j 

i 

Discussions  and  numerical  results  are  presented  in  Section  3.  In  the  ’ 1 
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2.  APPROXIMATE  FORMULA  FOR  MUTUAL  AIMITTANCE 

Referring  to  Figure  1,  consider  two  slots  on  the  surface  of  an 

infinitely  long  conducting  cylinder  with  radius  R.  The  orientation  of 

the  slots  may  be  either  circumferential  (Figure  lb  where  a^  > b^,  n * 1,2), 

or  axial  (Figure  Ic  where  a < b ) . The  problem  is  to  determine  the 

n n 

mutual  admittance  between  these  two  slots  when  kR  is  large. 

First  let  us  define  mutual  admittance.  Throughout  this  work  we 
always  assume  that 

(i)  the  slots  are  thin,  and  (2.1a) 

(il)  their  length  is  roughly  a half-wavelength.  (2.1b) 

Then  the  aperture  field  in  each  slot  can  be  adequately  approximated  by  a 
simple  cosine  distribution,  which  is  the  so-called  "one-mode"  approximation. 
For  example,  if  slot  1 is  circumferential  (Figure  lb),  its  aperture  field 
under  the  "one-mode"  approximation  is  given  by 


E - V^e^  , H - Ij^h^ 

(2.2a) 

where 

-►  A / 4 7t  ^ A 

T '^  « h^  =.  X X 

(2.2b) 

1 

ct: 

n 

(2.2c) 

(Vi.ip 

are  respectively  the  modal  (voltage,  current)  of  slot  1. 

The 

mutual 

admittance  Yj^2  defined  by 

^21 

Y = Y ■»  — ^ 

12  21  V, 

(2.3) 

where  1^^  is  the  Induced  current  in  slot  2 when  slot  1 is  excited  by  a 
voltage  and  slot  2 is  short-circuited.  An  alternative  expression  for 
is 


A 
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where 

A2  ® aperture  of  slot  2 

*>  magnetic  field  when  slot  1 is  excited  with  voltage  V^,  and 
slot  2 is  covered  by  a perfect  conductor 
^2  * electric  field  when  slot  2 is  excited  with  voltage  V2>  and 
slot  1 is  covered  by  a perfect  conductor. 

Because  ■ ^21^2  ^2  * ^2^2*  ^ simple  matter  to  verify  that 

(2.3)  and  (2.4)  are  equivalent  [16]. 

There  Is  an  alternative  definition  of  mutual  admittance.  Instead 
of  (2.2),  a modal  voltage  (with  a bar)  may  be  defined  through  the 
aq>res8lon  for  the  aperture  field  of  slot  1 as  follows; 


* 1 B IT 

E ■ z 7*  V,  COB  •—  y 
b 1 *1 

or  equivalently 


(2.5a) 


(2.5b) 


Then  a different  mutual  admittance  Yj^2  defined  by  (2.4)  after  replacing 
^'^1*^2^  by  (Vj^,V2).  It  can  be  easily  shown  that 


‘12 


'V2 

‘V2, 


12 


(2.6) 


Two  remarks  are  in  order:  (1)  In  the  limiting  case  that  b^^  and  b2  0, 

1/2  - 

Yj^2  goes  to  zero  as  (bj^b2)  , whereas  Yj^2  approaches  a constant 

Independent  of  b,  and  b,.  (11)  For  the  special  case  a^  = a,  = A/2  and 

^ ^ ^119 


R -»■  ",  it  is  Y^2»  not  Yj^2>  that  is  identical  to  the  mutual  impedance  Z 


between  two  corresponding  dipoles  calculated  by  the  classical  Carter's 
method  [5],  [8],  [9],  (Hi)  When  the  slots  are  excited  by  waveguides 
(transmission  lines),  one  often  uses  (^22).  From  here  on,  we  will 


12 


concentrate  on  Yj^2  instead  of  Y^2* 


For  the  two  slots  in  Figure  1,  the  final  form  of  an  approximate 
solution  of  Yj^2  is  as  follows  (for  exp  +jut  time  convention) : 
Circumferential  slots 


g i/2 

^12  ' 2 ^^l*^l®2^2^  ®^^1  nos  0)  g, 

TT  ^ J.  ^ ^ 

(2.7a) 


Axial  slots 


8 


w 1/2 

^2  ' " "T  SV2‘'2^  sin  0)  C(b„  sin  0)  g . 

TI  * * X ^ ”Z 


(2.7b) 


The  various  factors  in  (2.7)  are  explained  below.  S and  C are  simple 
trigonometric  functions 


^ • 

1 - (kx/ir) 


(2.8) 


The  (simplified)  Green's  functions  g^  and  g^  are  given  by 

\ = G(s)[v(0  (sin^  0 + ^ cos  2oj  + ^ u(0  cos^  0 
+ ju'(C)(/2  kR  cos  0)“^^^  sin^  o] 

h “ ® ~ ^ ®°®  ^ sin^  0] 


(2.9a) 

(2.9b) 


where 


G(s)  = 


k Y-  -jks 
0 e 


2Trj  ks 


. = 


0 I2O1T 


C = (k  cos^  0/2R^)^^^  s 


(2.10) 

(2.11) 
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6 

s - /zj  + (R^q)^  (2.12) 

e » tan  ^ (zQ/R(j)Q)  . (2.13) 

The  Fock  functions  u and  v are  explained  in  the  Appendix.  In  the  limiting 
case  kR  ->■  » (slots  on  a planar  surface),  (2.9)  is  further  simplified  to 
become 

g ■ G(s)  [sin^  e + ^ (2  - 3 sin^  6)]  , 

^ " kR  -*»  » . (2.14) 

® ^ ^ * 

The  formula  in  (2.4)  is  an  approximate  solution,  valid  under  the  condition 

kR  » 1 and  ks  » 1 . (2.15) 

The  numerical  accuracy  of  the  formula  is  discussed  in  Section  3,  and  its 


derivation  in  Section  4. 


> ->/s  . 
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3.  NUMERICAL  RESULTS  AND  DISCUSSION 


For  the  two  slots  In  Figure  1,  the  final  form  of  the  approximate 
solution  of  Y^2  given  in  (2.7).  Generally  speaking,  its  accuracy  is 
good  only  if 

(i)  the  size  of  the  slots  is  small  in  terms  of  wavelength, 
and/or 

(ii)  the  separation  of  the  slots  is  large  in  terms  of  wave- 
length. 

In  this  section,  we  will  give  some  numerical  examples  to  illustrate  the 
quantitative  accuracy  of  (2.7). 

(A)  Circumferential  Slot  - (Figures  2 and  3) . The  size  of  each  slot 

is  0.5X  X 0.2A,  and  the  cylinder  radius  is  lA.  7^2  is  presented  in  (dB, 
normalized  phase)  format,  where  dB  = 20  logj^Q  ^1^12^  normalized 

phase  is  equal  to  Arg(Yj^2  expjks).  Three  solutions  of  Yj^2  given:  the 
UI  exact  modal  solution  calculated  from  (5.2),  (5,3)  and  (5.9);  the  UI 
asymptotic  solution  reported  in  (1);  and  the  approximate  solution  in  (2,7), 

We  note  that  all  the  three  solutions  are  in  an  excellent  agreement. 

(B)  Percentage  Error  vs.  Slot  Position  - (Figures  4 and  5).  In  these 
figures,  the  coordinates  of  each  point  determine  the  center-to-center 
distance,  in  t and  z directions  between  two  slots.  The  pairs  of  numbers  in 
the  parentheses  r the  percentage  error  in  magnitude  and  the  absolute 
error  in  phase  of  Yj^2  calculated  by  the  approximate  formula,  respectively. 
For  the  circumferential  slots  (Figure  A),  the  accuracy  is  generally  very 
good.  For  the  axial  slots  (Figure  5),  the  approximate  formula  gives  erratic 
results  (as  high  as  27  percent  error  in  magnitude)  when  the  two  slots  are 
very  closely  displaced  in  the  ^i-dlrection.  The  reason  for  this  Inaccuracy 

is  that  the  surface  field  due  to  a magnetic  dipole  varies  very  rapidly  as 


I : 


t 

j 

i 


f 


i 

t 


f 


a function  of  z when  the  observation  point  is  close  by. 
122 


8 

(C)  Accuracy  vs.  Cylinder  Radius  (Figure  6).  The  accuracy  of  the 
approximate  formula  is  not  sensitive  to  the  radius  of  the  cylinder. 

(D)  Planar  Slots  (Tables  1 and  2).  The  mutual  admittance  between 
two  identical  slots  of  dimension  (a  = 0.69X,  b = 0.3X)  on  an  infinite 
conducting  plane  is  calculated  as  a function  of  Zq  and  y^  (the  center-to- 
center  distance  between  two  slots  in  z and  y directions,  see  Figure  lb). 

is  given  in  (dB,  phase  in  degrees).  In  both  E-plane  and  H-plane  couplings, 
the  approximate  formula  is  accurate  when  the  separation  is  at  least  two  wave- 
lengths (2.6").  It  should  be  also  remarked  that  the  present  slots  (0.69X  x 
0.3X)  are  relatively  large.  The  accuracy  of  the  approximate  formula  is  better 
when  the  slots  are  smaller. 


i'  "Ji'  . .-ii 


'>-v;',^:>s.'^-.><»srf4.'«£t?3«  ■; 


A.  DERIVATION  OF  APPROXIMATE  FORMULA 


We  will  now  give  the  derivation  of  the  formula  in  (2.7a) [that  of 
(2.7b)  is  very  similar].  Consider  a circumferential  infinitesimal  dipole 
located  at  Q'  on  the  surface  of  a cylinder  (Figure  7)  which  is  described  by 
the  magnetic  current  density 


K « (j,  I 6(r  - R)  6(z) 


(A.l) 


At  an  observation  point  Q on  the  cylinder,  the  (})-component  of  the  H field, 


denoted  by  g^,  is  determined  in  Eq.  (2.16b)  of  [1],  which  reads  in  the 


present  notation, 


g^(t,a)  G(t)  |j(C)^sin^  “ ^ ^“] 

+ [^)  u(C)  [cos^  a jl  - |i)  + sin^  a] 

+ j (/2  kR/cos^  a) 

« [v'(5)  sin^  a + |tan^  a + ;^|  u'(5)  cos^  a|^ 


(A. 2) 


where  (t,a)  are  the  cylindrical  coordinates  of  Q with  respect  to  the  origin 
at  Q'  on  a developed  cylinder,  and 


? = (k  cos^  e/2R^)^^^  t 


(A. 3) 


The  formula  in  (A. 2)  is  mainly  based  on  a classical  work  of  Fock  [17],  and 
contains  a modification  that  introduces  a field  dependence  on  the  surface 
curvature  in  the  binormal  direction  of  the  surface  ray  (see  Section  6 of 
[1]).  This  formula  is  asymptotically  valid  for  kR  «>,  and  may  be  used  to 
calculate  the  fiold  at  any  point  on  the  cylindrical  surface. 
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Making  use  of  the  Green's  function  In  (4.2),  we  next  calculate  the  surface 

field  H,  due  to  slot  1 on  a cylinder  (Figure  8).  The  aperture  dlstrl- 

<P 

butlon  of  slot  1 Is  described  In  '(2.2a),  which  may  be  replaced  by  an 
equivalent  magnetic  current  density  (p.  108  of  [18]) 


K » $ 6 (r  - R)  cos  (iry/a^^)  . 


(4.4) 


Then,  H.  at  an  observation  point  Q Is  obtained  by  superposition,  namely, 


ab  '^1 


cos  — y g. (t,a)  dy  dz 

I ®1  M 


(4.5) 


The  expression  for  calculating  the  mutual  admittance  between  the  two 
slots  In  Figure  8 Is  given  In  (2.4).  Note  that  Is  described  much  as 
(2.2a)  and  in  (4.5).  Then  (2.4)  becomes 


Y » 

* 1 f> 


-2 


ja  b a h J'A 

A 1 2 2 ^1 


I dy  dz  I dy2  dz2 


cos  — y 
\ ®1 


cos  — y2 
\ ®2 


g.{t,a)  . 


(4.6) 


The  distance  t in  (4.6)  Is  given  by 


2 1/2 
t ■ ((s  cos  e + yj  - y)  + (s  sin  0 + Z2  - z)]  ' 


(4.7) 


If  s Is  large  relative  to  the  length  of  either  slot,  t may  be  approximated  by 


t 


s 1 + cos  0 


72 -y 


+ sin  0 


Z2  - z 


(4.8a) 

(4.8b) 


In  evaluating  the  magnitude  of  g^  In  (4.6),  we  use  the  approximation  in 
(4.8a),  whereas  in  evaluating  Its  progressive  phase  term,  we  use  (4.8b). 
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5.  EXACT  MODAL  SOLUTION 


12 


The  admittance  Y^2  defined  in  (2.3)  may  be  calculated  exactly  by 
using  cylindrical  modes,  as  has  been  done  by  Stewart,  Golden  and 
Pridmore-Brown  [2],  [3].  Extensive  numerical  results  of  Y^2  calculated 
from  the  SGP  solution  are  reported  in  113),  [14].  As  will  be  explained 
below,  the  SGP  solution  is  not  suitable  for  numerical  calculations  when 
the  slot  separation  Zq  (Figure  la)  is  large.  In  this  section,  we  will 
derive  an  alternative  modal  solution  of  Yj^2  which  does  not  have  this 
difficulty. 

Let  us  first  consider  the  circumferential  slots  shown  in  Figure  lb. 
For  the  case  that  a^^  » 82  “ a and  b^^  ■ b2  ■ b (Identical  slots) , the 

A 

mutual  admittance  Y^2  given  In  Eq.  (8)  of  [3]  , which  reads  in  the 
present  notation. 


r ® -j  (tni}>o+k  z.) 

dk  I i^(m,k  )G(m,k  )e 

z z z> 

-00  m=-“ 


(5.1a) 


where  2 

sin  (kb/2) 

♦(«,k  ) = 5- 

^ Stt  R (k  b/2) 
z 


2 


^sln  (m(ji  + tt/2)  sin  (m(>>  - tt/2)^ 

+ — 7-: -7^^: — ) (5.1b) 


(m<ji  + v/2) 

3 


(mill  - 7i/2) 
3 


K = (a/2R) 

o 


G(m,k^)  = Yq 


jk 

H^^^'(k  R) 
m ^11 

mk  \ 
z 

2k 

t 

b(«(1c^R) 

■^t 

V.. 

H^^^(k  R) 
m t 

jk 

(5.1c) 


The  multiplication  factor  2 in  the  definition  of  iji^^  in  [3]  is  a misprint 
and  should  be  removed. 
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k.  = 


/ 2 2 

v-k  - k , if  k > k 

2 — Z 


i -j  /k^  - k^  , if  k < k 

V Z ’ — 2 


Rewrite  Y^2  terms  of  its  real  and  imaginary  parts: 


Y^2=G+jB  . 


It  can  be  shown  that  G is  given  by 


f 00 

= I 

^0  m=0 


“ cos  m<ii. 


— cos  k 2^  il<(m,k  )R(m,k  ) dk 

C Z \J  z z z 

m 


where 


2 ^ 1 mk  . 

R(m,k  ) - ~i."  p * * —————  + I-—-——  , 

2 TTk^R  kj.  R)  \K(.k  R N^(k  R) 

L m t . m t 

M^(X)  = J^(X)  + Y^(x) 

m m ''  m ^ 


N^(X)  = J'^(X)  + Y'^(x) 

mm  ra 


2 , m = 0 
1 , m 0 


(5.3a) 


(5.3b) 


(5.3c) 


(5.3d) 


(5.3e) 


We  note  that  G contains  a finite  integral  and  can  be  evaluated  in  a straight- 
forward manner  by  standard  numerical  integration  techniques.  The  imaginary 
part  of  Y,,  is  given  by 


I 

C , m= 


«>  cos  mi{i. 


m=0  m 


• cos  k z • <|>(m,k  ) • W(m,k  ) dk  (5.43) 


where  the  integration  contour  is  shown  in  Figure  9 and 
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Make  the  change  of  variable  k = jn  In  (5.8).  Substitution  of  the  resultant 

2 


equation  and  (5.7)  into  (5.5)  gives 

« cos  in(l>„  C • 

B-  I 

ni=0  *^ni  ' 


R(in,k^)tj)(in,k^)  sin  dk^ 


R(m,jn)>l'(m,jn)e 


-nz. 


(5.9) 


dn 


Our  final  expression  for  given  in  (5.2),  with  its  real  part  G in  (5.3) 

and  its  imaginary  part  B in  (5.9).  Several  remarks  are  in  order:  (i)  Not 

only  G but  also  B is  determined  by  R(m,k  ),  which  is  much  simpler  than  W(m,k  ) 

2 2 

defined  in  (5.4b).  (ii)  B contains  only  a finite  integral,  (iil)  The 
infinite  Integral  in  B,  i.e.,  the  second  integral  in  (5,9a),  contains  an 
exponentially  decaying  factor  exp[-ZQ  ~ a)n]  in  its  integrand.  The  emergence 
of  the  evaluation  of  B is  faster  for  larger  Zq.  This  is  in  contrast  to  the 
original  expression  of  given  in  (5.1).  (iv)  There  is  no  nonintegrable 
singularity  in  (5.3)  or  (5.9). 

The  same  method  applies  to  the  derivation  of  an  alternative  expression 
of  Yj2  for  two  identical  axial  slots  (Figure  Ic  with  a^^  = a2  = a and 
bj^  = b2  = b).  We  give  below  only  the  final  result: 


‘12 


abYn  «>  cos  m<t>rt 
0_  r 0 

2 ^ c 

TTk  R in=0  m 


lo 


-nz. 


+ j 'l>(m,jn)e 
^0 


?>(m,k  )e 
2 


dn 


R) 

m t 


(5.10a) 


where 


^(m.k^)  = 


“sin  (mi{'  ) 

3 

(mi}'  ) 

. O 


cos  (k^b/2) 


(k^b/2)^  - (ti/2)^J 


(5.10b) 
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In  summary,  the  alternative  expression  of  the  exact  modal  solutions  is  given 
in  (5.2),  (5.3),  and  (5.9)  for  two  Identical  circumferential  slots,  and  in 
(5.10)  for  two  identical  axial  slots. 
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APPENDIX 


FOCK  FUNCTIONS 


In  this  appendix  we  define  and  list  some  useful  formulas  of  the  functions 
Wj^(t),  w^(t),  v(0,  u(C),  and  Vj^(0.  These  functions  arc  commoniy  known  as 
Fock  functions. 

(i)  Definition;  For  a complex  t and  a real 


w,  (t)  = — dz  exp  (tz  “ 4 z^j 
^ h \ j ; 


1 r / 1 ^ \ ^ 

w„(t)  = — dz  exp  [tz  - TT  z I = w. (t) 


v(0  - i i I at 


ua).eJ'"V«i  ( ^.-Wat 
I'l-li.  I/O  1 f W~(t)  _tr,. 


(A-1) 


(A-2) 


(A-3) 


(A-4) 


JZ  iv 


where  integration  contour  rj^(r2)  goes  from  ® to  0 along  the  line 

Arg  z » -2it/3  (+2ti/3)  and  from  0 to  « along  the  real  axis,  Because  of 

different  time  conventions,  ”^(^2)  above  is  equal  to  w,(Wj^)  defined  In  [17]. 
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« ' 

^ : * 


(H)  Residue  scries  representation;  For  real  positive  C, 


v(5)  . 5‘'3 

* 1 -j^t' 

(A-6) 

u(0  « e^’'^^2/n 

i 

n*l 

(A-7) 

v^(5)  - e'^"^^24' 

« -jCt’ 

1 e " 

< 

(A-8) 

v'(0  - j I (1  - j2Ct’)(t’)"^e  (A-9) 

n=l  " " 

u'(0  - e^’'^S/ir  J (^  " J 3 

n*l 

where  {t^}  and  {t^  } are  zeros  of  W2(t)  and  w^(t),  respectively,  and 
are  tabulated  In  Il7]  and  11]. 

(ill)  Small  argument  asymptotic  expansion;  For  real  positive  5 and  ? •*'  0, 


v(0 . 1 . ^ * |i  .3  * ^ ,-3'«5’«  - 4.141  . 10-35*  * . ., 


.V,  1 _ e)’^/V3/2  5i  ,3  525  _ _ 


(A-11) 


,701  X 10"^e^  + ... 


(A-12) 


,^(5)  •,.  1 * f e3'/‘53/2  . a 53  . ^ ^-J./4^,/2  , ^ , j„-2^6  ^ 

(A-13) 

,</£)  ,v  ^ e-:13»'/^rl/2  . 21  2 63^  -jn/4.7/2  , ,.,  ^ ,.-2.5  ^ 

5 ^20^^  1024  ® ^ - 2.485  X 10  5 + . . . 

(A-14) 

u'tt)  ^ I .-13»/4,1/2  . a 53  + ^ .-J./4j;/2  , , ^„-l  5 . 


221  X 10" + ... 

(A-15) 


(iv)  Niiiiicrlcal  evaluation;  For  5 > Cq.  the  residue  series  representation 
with  the  first  ten  terms  in  the  summation  may  be  used.  For  C £ F,^,  the  small 
argument  asymptotic  expansion  with  the  first  five  terms  may  be  used.  It  has 


been  indicated  in  [12] that  the  smoothest  crossover  is  obtained  If  = 0.6. 
In  the  pie.'ii'ut  study,  we  set  F,q  " 0.7,  where  the  difference  in  the  two 
representai  ion.s  is  less  than  O.U  In  magnitude  and  0.9°  in  phase  [1  ]. 
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ABSTRACT 

An  approximate  asymptotic  Green's  function  for  the  surface  magnetic 
field  due  to  a magnetic  dipole  on  a general  convex  conducting  surface  is 
developed.  Based  largely  on  the  classical  work  of  V.  A.  Fock  and  the 
current  GTD  recipes,  this  solution  is  presented  in  a form  that  admits  ray 
interpretations,  and  can  be  simply  evaluated.  We  apply  the  Green's  function 
to  calculate  the  mutual  admittance  between  two  slots  on  a cone.  The 
numerical  results  are  in  very  good  agreement  with  experiments. 


* This  work  was  supported  by  Naval  Air  Systems  Command  under  Contract 
N00019-77-C~0127. 
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1.  INTRODUCTION 


In  the  past  few  years  there  has  been  an  increasing  interest  in  the 
study  of  conformal  arrays,  where  the  radiating  elements  are  arranged  on 
a curved  conducting  surface.  A crucial  as  well  as  challenging  problem  in 
the  design  of  a conformal  array  is  determining  the  mutual  coupling  among 
elements  [1],  [2].  The  present  paper  addresses  itself  to  one  such  problem 
when 

(i)  the  conducting  surface  is  an  infinite  cone; 

(ii)  the  radiating  elements  are  slots  with  thin  width  and  a length 
of  about  a half  wavelength;  and 

(iii)  all  the  elements  are  distributed  in  a region  which  is  away 
from  the  cone  tip  and  whose  radii  of  curvature  are  large  in 
terms  of  wavelength. 

Because  of  assumption  (ii),  the  aperture  field  of  a slot  can  be  well- 
approximated  by  a simple  cosine  distribution,  i.e.,  the  so-called 
"one-mode  approximation."  Then  it  has  been  established,  e.g.,  p.  53  of 
[3]  or  p.  8 of  [4],  that  the  calculation  of  the  mutual  coupling  is 
reduced  to  tl.it  of  a dyadic  Green's  function  due  to  a magnetic  dipole  on 
the  same  curved  conducting  surface  where  the  array  is  at. 

The  Green's  function  of  a cone  can  be  calculated  in  the  following 
two  ways.  !'  ,i;  , r.ormal  modes  involving  spherical  Bessel  functions  and 
associated  Legendre  functions,  it  can  be  expressed  exactly  in  terms  of  a 
doubly  infinite  series  [5],  [6].  The  numerical  evaluation  of  such  a 
series,  however,  is  quite  tedious,  especially  at  high  frequencies.  Thus, 
up  to  now,  no  systematic  numerical  results  have  been  generated  from  the 
series.  An  alternative  way  to  calculate  the  Green's  function  is  to  employ 
the  surface  rays,  which  would  yield  a simple  approximate  solution  valid 
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for  high  frequencies.  The  general  concept  of  surface  rays  was  introduced 
by  J.  B.  Keller  more  than  twenty  years  ago  [7]  - [9];  however,  a uniformly 
valid  formula  for  fields  on  the  ray  has  not  been  developed  until  recently. 
Among  the  several  comparable  formulas  [10]  - [12],  we  chose  the  one 
reported  in  [12]  for  the  present  application.  The  reason  for  our  choice 
is  that,  at  least  for  the  case  of  a cylinder,  the  formula  in  [12]  gives  the 
most  accurate  numerical  results  [13]. 

The  organization  of  this  paper  is  as  follows.  The  formula  in  [12] 
for  the  Green's  function  applies  only  to  a cylinder.  Following  the 
GTU  recipe,  we  generalize  it  to  an  arbitrary  convex  surface  with  its  final 
solution  presented  in  Section  2,  and  its  derivation  in  Section  6.  In 
Sections  3 to  5,  the  Green's  function  is  specialized  to  a cone  and  is  used 
to  calculate  the  mutual  admittance  between  two  slots  on  a cone.  A 
conclusion  is  given  in  Section  7.  The  two  appendices  contain  (A)  formulas 
for  the  Fock  functions,  and  (B)  the  computer  listing  for  calculating  the 
mutual  admittance  on  a cone. 
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2.  SOLUTION  OF  THE  GREEN'S  FUNCTION 

Consider  a perfectly  conducting  convex  surface  E (Fig.  1),  whose 
radii  of  curvature  at  any  point  are  large  in  terms  of  wavelength.  At  a 
point  described  by  position  vector  r^^  on  £,  there  is  a tangential 
magnetic  dipole  source  described  by  a magnetic  current  density  (for 
exp  + jwt  time  convention) 

K(?)  = M6(?  - (2.1) 

when  M is  the  magnetic  dipole  moment  and  lies  in  the  tangent  plane  of  E. 

The  problem  is  to  determine  a high-frequency  asymptotic  solution  of  H at 
a general  point  Q2  described  by  position  vector  ^2  on  E.  In  other  words, 
the  dyadic  Green  s function  for  the  surface  magnetic  field  for  points  r^^ 
and  r2  is  to  be  found. 

Before  presenting  the  solution,  let  us  introduce  several  definitions 
and  parameters.  According  to  GTD  18],  [9],  the  dominant  high-frequency 
contribution  to  H(r2)  is  the  field  on  the  surface  ray  from  r^^  to  r2.  The 
surface  ray  is  a geodesic  of  E.  Some  of  its  geometrical  properties  are 
described  by  .'ig.  1) 

(i)  the  arc  length  s which  is  chosen  such  that  s = 0 at  the  source 

-V 

point  and  s = s at  the  observation  point 

(li)  thf  tangent,  normal,  and  blnormal,  denoted  by 

r where  n = 1,2;  and 
n 

(ill)  its  two  radii  of  curvature  R^(s),  and  R^(s)  of  E at  point  s 

in  the  directions  of  tangent,  and  binormal,  respectively. 

(Un  a general  convex  surface,  both  radii  are  nonnegative.) 

From  the  above  parameters,  we  may  calculate  the  following  quantities 


that  are  needed  for  the  solution  of  the  Green's  function: 
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(i)  The  large  parameter  in  our  asymptotic  expansion  of  the  Green's 
function  is 

U/3 


m(s)  = 


(2.2) 


which  is  a function  of  position  along  the  ray  from  r^^  to 

(ii)  A distance  parameter  from  r^  to  is  defined  by 


2m^(i) 


ds 


(2.3) 


For  the  special  case  when  R^.  is  not  a function  of  s (constant  ray  curvature), 

2 

C is  reduced  to  (ks/2m  ),  a well-known  parameter  introduced  first  by 
Fock  (14]. 

(ill)  The  ray  curvatures  at  the  source  and  observation  points  enter  in 
a parameter  defined  by 


X 


ks 


11/2 


2m(0)  ro(s)  £ 


(2.4) 


which  is  positive  real  for  a convex  surface,  and  is  reduced  to  unity  for 
the  special  case  of  a constant  ray  curvature. 

(iv)  Consider  a small  pencil  of  surface  rays  originating  from 
and  propagating  toward  r2  (Fig,  1),  The  angle  extended  by  the  pencil  at 
r^^  is  di|jp  and  tliat  at  r2  is  d<^2'  The  divergence  factor  DF  of  the  pencil 
is  defined  by 


DF  = 


(2.5) 


where  p is  the  caustic  distance  of  the  wavefront  at  r2  and  is  always 

> 

positive.  For  example,  if  T is  a sphere  and  rj  is  the  nortli  pole,  DF  at 
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point  T2  = (r,0,(l))  is 


which  varies  from  one  at  the  north  pole  (0  = 0)  to  infinity  at  the  south 
pole  (0  = Ti)  as  T2  moves  along  a great  circle. 

(v)  The  "mean"  radii  of  curvature  between  r^^  and  r2  are  defined  by 

= [R^(0)  R^(s)]^^^  (2.6i 

= (Rj^(O)  R^(s)]^/^  . (2.61 


Throughout  this  work,  we  always  assume  that  E is  a smooth  surface  with  a 
slowly  varying  curvature.  Then  (Rj_,Rj^)  represents  a sort  of  average 
value  of  radii  of  curvature  along  the  ray. 

Return  to  the  electromagnetic  problem  in  Fig.  1.  We  assume  that 
m(s)  is  large  and  is  slowly  varying  for  all  s in  the  range  0 < s < s. 

Then  an  approximate  asymptotic  solution  for  the  surface  magnetic  field  at 
r2  due  to  the  dipole  source  in  (2.1)  is  given  by 


where 


H(r2)  “ M • (b’bUjj  + t’tHj.)(DF) 

= G(s)|Jl  - Tv(C)  - t\(0  + j(/2kR^) 


(2. 


o-..n  ^-2/3 


,3 


[Tv'(C)  + (Rj./Rj,)  < u'(0] 


“t  = 


i\U)  + 1 


1-li) 


(C)  + j(.^kR|.)“^^^  T^u'(U 


G(s)  = 


k^Y 


-1 


Inj  ks  ’ 


Y = (12071) 


6 


K 


Ji 


The  Fock  functions  u and  v and  their  derivatives  u'  and  v'  are  described  in 
Appendix  A.  Several  remarks  about  the  solution  in  (2.7)  are  in  order. 

(i)  It  is  derived  in  an  approximate  manner  from  the  classical  work  of 
Fock  [14]  and  the  recipe  of  GTD  [8],  [9],  as  detailed  in  Section  6.  This 
solution  is  certainly  not  valid  when  the  curvature  of  the  surface  E is 


large  or  rapidly  varying,  (ii) 

For  the  special  case 

that  E is 

a planar 

1 

surface  (R^,  = Rj^  ->■  "),  (2.7)  recovers  the  known  exact 

solution. 

namely. 

1 

2 

•i 

= G(s) 

■ A ■ fe] 

(2.8a) 

f 

< 

Hj.  = G(s) 

#)  k ■ i) 

(2.8b) 

f 

DF  = 1 

(2,8c) 

(iii)  The  solution  is  valid  for  any  combination  of  r^^  and  V2>  In  the 


penumbra  region  (r^  is  close  to  and  C <<  1),  (2.7)  is  nearly  the  planar 
solution  in  (2,8).  In  the  deep  shadow  (C  <<  1),  the  residue  series 
representation  of  the  Foci'  functions  in  Appendix  A may  be  used,  and  (2.7) 
is  identified  as  the  creeping-wave  contribution,  (iv)  When  E is  a 
cylindrical  surface,  the  formula  (2.7)  has  been  used  to  calculate  the 
mutual  admittance  between  two  slots  on  E.  It  has  been  shown  [12],  [13] 
that  the  numerical  results  are  in  excellent  agreement  with  a known  exact 
solution  [15]  - [17].  (v)  Except  for  the  very  simple  surfaces  such  as  a 

cylinder,  cone,  or  sphere,  no  explicit  parametric  equations  can  be 
found  for  the  geodesics.  Thus,  for  a general  surface,  one  may  have  to 
rely  on  numerical  techniques  for  determining  the  geodesics  and  the 
divergent  factor,  as  has  been  done  in  [18], 
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3.  GREEN'S  FUNCTION  OF  A CONE 

Let  us  apply  the  formula  (2.7)  to  the  field  on  an  infinite  cone, 
described  by  the  equations  (Fig.  2a) 

X = r sin  6^  cos  ifi  , y = r sin  6^  sin  <j:  , z = r cos  0^  (3.1) 

where  6^  is  the  half  cone  angle  (0  ' < Tr/2).  Since  the  cone  is  a 

developable  surface,  the  rays  (geodesics)  on  a developed  cone  (Fig.  2b) 
are  straight  lines.  Due  to  the  source  at  r^  = main  con- 

tribution of  the  field  at  ^2  = (r2,6Q,<j>2)  comes  from  the  shortest  ray 
described  by 

r^  sin  = r2  sin  ^2  • (3.2) 

As  the  ray  propagates  away  from  the  source  point  r^,  it  reaches  the  highest 
altitude  at  M where  fi  = v/2.  After  M,  the  ray  travels  downward  away  from 
the  cone  tip.  The  various  parameters  defined  in  Section  2 can  be  simply 
calculated  from  the  cone  geometry  [6],  [19],  and  expressed  in  terms  of 
coordinates  and  (r2,<ti2)*  The  arclength  is 

s = (r^  + r2  - 2rjr2  cos  [(<)>j^  - sin  • (3.3) 

•*> 

The  angle  52^  at  r^^  is 

' sin  ^ " '*’1^ 

2 2 2 

We  choose  (i!jj  7i/2  if  r2  < s + rj^  , and  > 7r/2  if  otherwise.  The 

other  parameters  are 
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^2  = ^1+  (‘l>2  ~ '*’1^  ®o 


5 = Z!i!2_!!!L!o  5 Ai^2  % 

t sin  sin  ^2  ’ “ cos  cos  ^2 


C = ( 2 >^^2  sin  sin  Q^j  |<J.2  - <(>^1 


2/3  , 
cos  0, 


1/2/  2 \ ■^'”  ] 
T = (ks/5)  ' 2k  ^2^2)  (sin  sin  ri^  cot  6^) 


DF  = 1 


When  the  above  parameters  in  (3.3)  through  (3.9)  are  substituted  into  (2.7), 
we  obtain  an  approximate  solutioi.  for  the  surface  field  on  a cone  due  to 
a direct  surface  ray  contribution.  Let  us  consider  a special  observation 
point  r2  such  that 

ks  >>  1 , and  fi2  close  to  v/2  , (3.10) 

Then  the  two  components  of  the  field  in  (2.7)  are  reduced  to,  after  making 
used  of  the  residue  series  renreseatations  for  the  Fock  functions 
(Appendix  A)  and  keeping  only  the  leading  terms. 


k^(sin  Q.  sin  .T,  cot  f*  /r 

— — n — ir/b-Ti7r-  - ^[ti  ^ ^ H n n . 

1528  1k  rj^r2l  (kr.)  L.  ' U(3.11a) 


^ 0|(ks)'^^‘] 


(3.11b) 


which  agrees  with  the  rigorous  asymptotic  solution  given  in  Eqs.  (50)  and 
(53)  of  [6].  (In  making  the  comparison,  note  tlie  corresponding  notations 
used  in  (6)  and  here:  -i  -►  j,  6^  9q>  -*■  s,  r^  -*■  rj^,  r^  -»  r2, 

-»  ti/2  - and  |t^|.)  We  emphasize  that  the  result  in  (3.11)  or 

that  in  (6]  is  valid  only  under  the  conditions  in  (3.10).  For  an 
arbitrarily  located  observation  point,  (2.7)  should  be  used. 
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Two  final  remarks  about  the  formula  in  (2.7)  are  in  order,  (i)  For 
a given  source  and  observation  point,  there  are  infinitely  many  rays 
(geodesics)  passing  through  them.  The  contribution  from  each  ray  may  be 
calculated  from  (2.7),  and  the  final  field  solution  is  the  superposition 
of  all  ray  contributions.  In  most  practical  problems  (all  the  numerical 
computations  presented  in  this  paper),  only  the  ray  with  the  shortest 
arclength  gives  the  significant  contribution  to  the  field  solution, 
whereas  all  other  rays  may  be  ignored,  (ii)  Depending  on  the  polarization 
and  the  distances  of  the  source  and  observation  points  from  the  cone  tip, 
there  may  be  another  significant  contribution  to  the  field  from  the 
diffraction  at  the  tip.  In  such  a case,  the  total  field  at  any  point 
contains  two  dominant  contributions;  one  from  the  direct  ray  according 
to  formula  (2.7),  and  one  from  the  tip-diffracted  ray.  More  about  the 
latter  will  be  given  in  Section  4. 
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4.  MUTUAL  ADMITTANCE  ON  A CONE 

On  the  surface  of  a cone,  let  us  consider  two  arbitrarily  oriented 
slots.  Under  the  assumption  that  the  dimensions  of  the  slots  are 
relatively  small  compared  with  the  radii  of  curvature  of  the  cone  surface, 
the  shapes  of  slots  are  taken  to  be  rectangular  on  a developed  cone. 

Note  that,  depending  on  the  exact  manner  in  which  the  feeding  waveguide 
is  fitted  into  the  cone  surface,  the  shape  of  a slot  mapped  on  a developed 
cone  may  be  quite  irregular.  Our  assumption  of  rectangular  shapes 
represents  a good  approximation  for  practical  cases;  at  the  same  time, 
it  simplifies  the  subsequent  calculations. 

Referring  to  Fig.  3,  we  describe  the  dimensions  and  the  positions 
of  the  two  slots  by 


(a^,b  ) and  [c  , (n 
n n n 


1)  ‘*>0’“’^  ’ n = 1,2 


Thus,  the  radial  separation  of  the  two  slots  is  (c^  - c^)  and  the  angular 

separation  is  41^.  The  angle  measures  the  deviation  of  the  longitudinal 

direction  of  slot  n from  the  radial  direction  of  the  cone.  If  w = 0, 

n 

slot  n is  radial;  if  u = ii/2,  slot  n is  circumferential.  The  mutual 

n 

admittance  Yj^2  between  tlie  two  slots  is  defined  as  follows.  Throughout 
this  work  we  always  assume  that 

(i)  the  slots  are  thin,  and  (4.1a) 

(ii)  their  length  is  roughly  a half-wavelength.  (4.1b) 

Then  the  aperture  field  in  cacli  slot  can  be  adequately  approximated  by  a 
simple  cosine  distribution,  which  Is  the  so-called  "one-mode"  approximation 
The  aperture  field  of  slot  1 under  the  "one-mode"  approximation  is  given  by 


M.q 


H.  ,,hj 


(4.2a) 


I-  ■ .ail 


11 


where 


Xixei 


(4.2b) 


Here  are  the  local  rectangular  coordinates,  with  the  origin  at  the 

center  of  slot  1 and  y^-axis  parallel  to  the  longer  dimension  of  the  slot. 
(Vi,Ii)  are  respectively  the  modal  (voltage,  current)  of  slot  1.  The 
mutual  admittance  Y^2  defined  by 


Y = Y 
12  ^^21 


(4.3) 


where  is  the  Induced  current  in  slot  2 when  slot  1 is  excited  by  a 
voltage  and  slot  2 is  short-circuited.  An  alternative  expression  for 


= aperture  of  slot  2, 

» magnetic  field  when  slot  1 is  excited  with  voltage  V^,  and 
slot  2 is  covered  by  a perfect  conductor, 

E2  *=  electric  field  when  slot  2 is  excited  with  voltage  V2.  and 
slot  1 is  covered  by  a perfect  conductor. 

Because  = 121^2  ^2  ~ '^2^2*  ^ simple  matter  to  verify  that 

(4.3)  and  (4.4)  are  equivalent. 

At  high  frequencies,  Hj^  in  (4.4)  has  two  dominant  contributions: 
one  from  the  direct  rays  going  from  siot  1 to  slot  2,  and  the  other  from 
the  rays  diffracted  at  the  tip  of  the  cone,  namely, 
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Hi  % + Hi 


(4.5) 


Accordingly,  Y^2  also  has  two  parts 


^12  ^12  ■*■  ^12 


(4.6) 


Let  us  concentrate  on  Y°2  first.  Making  use  of  the  Green's  function  in  (2.7) 


and  the  aperture  distribution  in  (4.2),  Y°2  may  be  explicitly  written  as 


a,/2  b,/2  a,/2 


12 


(a, 


u,  / ay.  / £. 

f 1 f 2 

dYl 

dz  dy 

-b,/2  ^ J-a,/2  ^ 

J-a,/2  ^ J 

dz„ 


|cos  f- y^j  ||cos  f- 


(4.7a) 


where 


8(yj^.Z|!y2*Z2)  = cos  cos  sin  sin 


(4.7b) 


The  Green's  function  components  (•*jj»Hj.)  are  given  in  (2.7),  and  angles 
(oiyoi^)  are  shown  in  Fig.  3.  In  evaluating  the  integrals  in  (4.7a),  for 
two  given  points  (yj^,Zj^)  and  (y2,Z2)>  '*'a  must  calculate  some  geometrical 
parameters  appearing  in  and  H^.  Those  calculations  lead  to  the 
following  results 


r = 
n 


/~9~  9 

/y  +7." 


2 2’ 
c + y + z"  - 2c  / y~  + 7~  cos  (w  - to  , , ) 
nn  n n.^n  n n n+4 


1/2 


= (sin  Oq)  ^ sin  ^ 


/ 2 ^ 2 -1  . , , 
/ y + z r sin  (to  - to  , , ) 
/ ■'n  n n n n+4 


(4.8a) 

(4.8b) 


to  , , = tan  ^ (z  /y  ) 
n+4  ' n n 


<*<  .o  “ + (n/2)  - to  - ({i  sin  0„  + (n  - 1)  it>,,  sin  0. 

n+2  n n n 0 0 0 


(4.8c) 

(4.8d) 


where  n = 1 and  2.  We  evaluate  the  integrals  in  (4.7a)  numerically  with  the 


aid  of  a computer. 
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Next  let  us  consider  Y^2»  part  of  mutual  admittance  due  to  the 
rays  diffracted  at  the  cone  tip.  For  the  special  case  of  circumferential 
slots,  an  approximate  expression  of  given  in  [16],  which  reads 


yJ^2  “ T,  if  0)^  = 0)2  “ ’^''2 


(4.9a) 


where 


T = o 


(3^6^3262)  /tan 


sin  (kb,/2)  sin  (kb*/2) 


^ 30tt^Cj^C2  sin  6^ 


exp  - kc^  - kc2j 


2it 


(kb^/2)(kb2/2) 


(4.9b) 


Here  is  the  zeroth-order  tip  diffraction  coefficient  and  is  a function 
of  the  half  cone  angle  0^.  A numerical  table  of  for  several  typical 
values  of  0q  is  given  in  [16].  We  have  fitted  those  values  by  a simple 
expression,  viz., 


Oq  = A exp  jB  , (4.10) 

where 

A = 1.30570”^  - 1.755  + 2.7720^  - 1.4590^ 

B = 2.7195  + 1.46O80Q  - 1.12950^  + 0.65660^ 

Both  0 and  B are  in  radians. 

As  may  be  seen  from  Fig.  4,  the  numerical  values  of  o^  calculated  from 
(4.10)  are  in  excellent  agreement  with  those  tabulated  in  [16].  For 
the  special  case  of  axial  slots,  Y^2  according  to  [16]  is  approximately 

y52  ' 0 , if  o)j  = 0)2  = 0 . (4.11) 

In  the  present  paper,  we  are  Interested  in  the  general  case  that  the  two 
slots  have  arbitrary  orientations.  Before  a more  exact  solution  can  be 
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found,  we  will  use  the  following  formula 


T sin  sin 


(4.12) 


which  matches  the  two  extreme  cases  in  (4.9)  and  (4.11),  and  interpolates 
the  in-between  cases  by  regarding  each  slot  as  a thin  magnetic  dipole. 
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5.  NUMERICAL  RESULTS  OF  MUTUAL  ADMITTANCE 

The  final  solutions  for  (total  mutual  admittance)  and  Y^2 
(partial  mutual  admittance)  are  given  in  (4.6),  (4.7),  (4.12),  and  (4.9b). 
For  a given  geometry  of  the  slots  and  cone,  the  two  surface  integrals  in 
(4.7a)  are  evaluated  numerically  by  choosing  an  integration  grid 
roughly  equal  to  0.05A  x 0.05A.  Thus,  for  two  typical  0.5A  x 0.2A  slots, 
the  integrals  are  replaced  by  a summation  of  1600  terms.  Fortunately, 
the  integrand  is  simple  enough  that  each  Y^^^  calculation  takes  about  one 
second  on  the  CDC  Cyber  170  Series  Computer  Systems.  We  have  analyzed  a 
large  number  of  cases  for  Y^2‘  Typical  results  are  summarized  belo". 

Unless  specified  otherwise,  all  numerical  computations  are  based  on 
two  identical  slots  with 

slot  length  = 0.5A  , width  = 0.2A  , (5.1) 

The  other  parameters  are  the  half  cone  angle  0^,  the  slot  orientations 
(a)j^,a)2),  the  distances  from  slot  centers  to  the  cone  tip  (Cj^,C2),  and 
the  slot  angular  separation  ‘^'q  ‘(’2  ~ (fig-  3). 

(A)  "Equivalent"  cylinder.  It  has  been  conjectured  in  [16]  that, 
in  calculating  Y^2  (•^he  contribution  from  the  direct  rays)  approximately, 
the  cone  may  be  replaced  by  an  "equivalent"  cylinder  with  parameters 
(Fig.  5) 

^0  " ^2  ’ ‘*’0  ' *^0  ’ i ^‘^l 


This  conjecture  can  be  now  quantitatively  checked  out.  In  Table  I we 
compare  Y^2  ^ cone  with  6^  = 15°  or  30°  calculated  from  (4.7),  and 

Yj^2  ^ cylinder  calculated  by  a similar  GTD  solution  reported  in  [12]. 
All  values  of  Yj^2  (o’^  ^22^  listed  in  (DB  = 20  log^^Q 
degree)  format.  For  the  cone  with  the  smaller  angle  (6^  = 15°),  the 


165 


16 


"equivalent"  cylinder  method  gives  a good  approximation  for  the 

cases  listed  in  Table  I,  the  magnitude  error  of  Y^2  within  0.5  dB 
(6  percent)  and  phase  error  within  15“.  For  the  cone  with  the  larger  angle 
(6q  = 30°),  however,  the  "equivalent"  cylinder  method  is  not  very  accurate 
with  magnitude,  and  phase  errors  as  large  as  2.5  dB  (33  percent),  and  56°, 
respectively. 

(B)  Comparison  with  experiments.  A set  of  experimental  data  on  the 
mutual  coupling  between  two  X-band  open-ended  waveguides  (0.9"  x 0.4") 
on  a cone  was  reported  in  [16].  As  a function  of  frequency,  measurements 
were  done  on  the  coupling  coefficient  which  is  related  to  Yj^2  through 

the  relation 


-Y  y,- 

S U.  ..  - 

^ - '‘h 


-1  2 1/2 

Here  Y^  = (120it)  [k“  - (n/a)  ] is  the  admittance  of  the  mode  in 

the  feed  waveguide.  Y^^^  is  the  self-admittance  of  a slot  on  the  cone. 

In  the  present  calculations,  we  use.  Instead,  the  Yj^j^  on  an  "equivalent" 
cylinder  which  is  calculated  by  the  exact  modal  solution  described  in  (16] 
(for  example,  “ 0.8178  + jO.3886  at  8.5  GHz,  and  0.8591  + jO.3828 

at  9 GHz).  Since  Y^^  is  least  sensitive  to  the  geometry,  tlie  approximation 
of  a cone  by  a cylinder  should  not  introduce  any  significant  error  in  3^2* 
In  Figs.  6 and  7,  three  sets  of  data  are  presented:  (i)  the  experimental 


data;  (ii)  the  theoretical  results  from  the  present  analysis  in  wliich  tlie 
calculation  of  Y^  is  based  on  a cone,  e.g.,  Equation  (4.7);  (iii)  the 
theoretical  results  from  [16]  in  which  Y^  is  calculated  from  tlie  exact 
modal  solution  of  an  "equivalent"  cvlinder.  Several  observations  can  be 


made.  (a)  Both  theoretical  results  are  in  good  agreement  with  the 
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experimental  data  (with  the  present  result  being  slightly  better) . As 
explained  in  (A) , the  "equivalent"  cylinder  method  works  because  the  cone 
angles  (6^  10°)  are  small,  (b)  The  peaks  and  valleys  are  caused  by  the 

interference  between  Y^2  ^12’  comparable  magnitudes  due 

to  the  large  angular  separations  (60.8°  and  80°).  (c)  There  exists  a 

slight  shift  in  frequency  (Af/f  = 3 percent)  between  the  theoretical  and 
experimental  valleys  in  Fig.  6.  We  speculate  that  this  may  be  due  to  a 
slight  phase  inaccuracy  in  As  a final  remark,  it  has  been  found 

experimentally  (private  communication  from  G.  E.  Stewart  and  K.  E.  Golden 
of  Aerospace  Corporation)  that  the  Y22  contribution  is  sensitive  to  the 
exact  shape  of  the  cone  tip.  When  the  tip  is  not  extremely  sharp,  the 
peaks  and  valleys  in  Figs.  6 and  7 become  much  less  predominant. 

(C)  Mutual  admittances  of  circumferential  slots.  In  Figs.  8 to  10, 
Yj^2  circumferential  slots  are  displayed  as  functions  of 

angular  separation  <{>q  and  the  radial  separation  (c^  - C2).  We  note  that 
the  effect  of  Y22  can  modify  the  curves  of  Y22  in  several  different  ways. 
When  the  slots  are  at  the  same  latitude  (Fig.  8) , the  direct  coupling  is 
weak.  Thus,  Lrp  contribution  is  noticeable  even  at  a small  angular 
separation.  As  the  radial  separation  is  Increased  (Fig.  9),  the  tip 
contribution  is  almost  negligible  for  (J>q  < 65°.  When  the  two  slots  are 
widely  separai  in  the  radial  direction  with  one  slot  near  the  tip 
(Fig.  10),  the  tip  contribution  gets  stronger,  and  the  direct  contribution 
becomes  insensitive  to  iJiq.  Hence,  the  oscillation  on  the  Yj^2  curve  has 
a much  larger  period.  In  fact,  there  is  only  a half  "cycle"  in  the  range 
0 < i})q  < 90°,  and  Y^^  appears  to  be  shifted  from  Y22  by  a fixed  amount. 
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(D)  Effect  of  slot  orientation  on  mutual  admittance.  Consider  two 
slots  separated  by  1 A along  the  radial  direction.  The  magnitude  of  ¥^^2 
as  functions  of  the  slot  orientation  angles  and  is  plotted  in 
Fig.  11.  As  expected,  the  maximum  value  (-73  dB)  occurs  when  both  slots 
are  circumferential  = 90°).  This  value  is  above  14  dB  higher  than 

that  when  both  slots  are  radial  ~ *^2  ~ minimum  value  (-113  dB) 

of  ¥^2  occurs  when  the  top  slot  is  radial  and  the  bottom  one  is 
circumferential.  This  result  confirms  a common  belief  that  the  mutual 
coupling  between  two  orthogonal  slots  is  generallv  negligible. 
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6,  DERIVATION  OF  THE  GREEN'S  FUNCTION 


We  will  now  describe  briefly  the  derivation  for  the  Green's  function 
for  the  general  convex  surface  given  in  (2.7). 

Our  starting  point  is  the  corresponding  Green's  function  for  a cylinder. 
For  the  latter  problem,  several  versions  of  the  asymptotic  solutions 
[10]  - [12]  exist.  All  the  versions  are  of  similar  nature,  and  contain 
some  approximations  that  have  not  yet  been  fully  justified.  For  the 
cylinder  problem,  both  [11]  and  [12]  give  excellent  numerical  results 
(with  [12]  being  slightly  better  as  demonstrated  in  [13]).  We  quote 
below  the  Green's  function  for  a cylinder  reported  in  [12],  which  again 
can  be  written  in  the  form  of  (2.7)  with  DF  = 1 and 


where 


n 


G(s)\  1 


- ± 
ks 


• [v*(0  + (R 


(O  - U(C)  + j(»^  kRj.)"^^^ 

,/«b)  -’(s)i} 

v(C)  + u(0  + j(^  u'(C) 


(6.1a) 

(6.1b) 


2 2 
R^  = R/cos  0 , R^  “ R/sin  0 

C = 1^8(003*^  6/2kV)^''^ 

R = radius  of  cylinder 

0 = angle  between  the  ray  and  the  (Ji-direction. 


Since  0 is  a constant  along  a ray  (geodesic),  so  are  the  three  parameters 
^t’  (6.1).  The  formula  (6.1)  is  basically  derived 

from  Fock’s  classical  solution  for  vector  potentials  of  a sphere  [14], 
but  contains  a modification,  namely,  the  last  term  in  (6.1a) 


G(s)  j(^  kR|.)"^''^(Rj./R^)  u'(0 


(6.2) 


was  added  to  the  Fock's  solution  in  an  arbitrary  manner.  Note  that  this 
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additional  term  introduces  a very  small  contribution  for  a sphere,  or  for 
a cylinder  as  long  as  6 is  not  close  to  ir/2.  For  0 = Tr/2  on 
a cylinder,  in  (6.1a)  becomes  (see  Eq.  (2.10)  of  [12]) 


Here  is  the  corresponding  solution  on  an  infinite  plane  and  it 

decays  as  (ks)  ^ for  large  ks.  The  second  term  in  (6.3)  comes  from  the 

additional  term  in  (6.2),  and  is  the  dominant  contribution  at  large  ks 

-1/2 

because  it  decays  as  (ks)  . Recently,  J.  Boersma  (private  communication) 
has  shown  that  (6.3)  is  in  exact  agreement  with  a rigorous  asymptotic 
solution  for  the  cylinder.  Thus,  the  additional  term  in  (6.3)  is 
justified  fcr  6 = ir/2  where  its  contribution  is  most  significant. 

Now,  let  us  generalize  (6.1)  to  a general  convex  surface  sketched 
in  Fig.  1 following  the  GTD  recipe  [3],  (9),  [18]:  (1)  The  divergence 

factor  in  (2.5)  is  introduced  from  the  consideration  of  energy  conservation, 
(ii)  The  generalized  C in  (2.3)  and  x in  (2.4)  are  based  on  the  rigorous 
solutions  of  two-dimensional  canonical  problems.  The  only  remaining 
problem  is  the  generalization  of  R^  and  Rj^  in  (6.1).  We  note  that  the 
above  GTD  recipe  is  valid  only  if  and  Rj^  are  slowly  varying  along  a ray. 
Thus,  any  sort  of  average  values  of  R^  and  R^  should  give  approximately 
the  same  result.  We  choose  the  geometrical  mean  in  (2.6)  for  its 
simplicity  and  symmetry  between  the  source  and  observation  points.  From 
these  considerations,  we  obtain  the  generalized  solution  (2.7)  from  (6.11. 
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TABLE  I 


ON  CONE  ANDY^2 


^0 

z A 

0 

Cone 

6 = 15° 

0 

30° 

0 

- 91.35dB 

155® 

0 

0 

0 

-110.97 

116® 

0 

1 

- 73.35 

73® 

0 

4 

- 83.89 

77® 

0 

0 

1 

- 86.52 

- 74® 

30® 

4 

- 86.95 

24° 

60® 

1 

-104.64 

-35® 

0 

0 

4 

- 94.41 

-118® 

The  parameters  are  * W2  “ 90®,  R = 


•k 

ON  CYLINDER 


cylinder 

Cone 

e = 30° 

0 

- 91.47 

- 90.99 

153® 

160® 

-111.28 

-108.78 

101® 

151® 

- 73.64 

- 73.49 

73® 

74® 

- 84.30 

- 84.06 

75® 



75® 

- 86.59 

- 86.60 

- 76® 

-70® 

- 87.14 

- 86.83 

20° 

35® 

-104.17 

-42® 

- 94.64 

- 94.12 

-133® 

-77® 

2A,  a = 0.5A,  and  b = 0.2A. 


f 


I 


j 

V 

I 

I 


J.  Two  icctangular  slots  on  a developed  cone. 
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Mil 'll!  iiirai  1 1 Jiu-es  on  a cono  calculaiod  from  (JTl)  lit'sci  ibed 
Sci^lion  -t. 


APPENDIX  A 
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FOCK  FUNCTIONS 


In  this  appendix  we  define  and  list  some  useful  formulas  of  the  functions 
Wj^(t),  W2(t),  v(0>  u(5)»  and  These  functions  are  commonly  known  as 

Fock  functions. 

(i)  Definition;  For  a complex  t and  a real 


dz  exp  Itz  ~ -j  z 


W2(t) 


— dz  exp  Itz  - 2^]  = w- (t) 


v(C)  - i -i- 

^ A 


u(«  . i 

/n 


a. 


P WjCt) 

'l  ^ 


p W,(t) 
il  z 


Vj^(0 


gj3ir/A  3/2  _1. 


-1£t 

1 


(A-1) 


(A-2) 


(A-3) 


(A-A) 


(A-5) 


where  integration  contour  rj^(r2)  goes  from  “ to  0 along  the  line 

Arg  z = -2ti/3  (+2tt/3)  and  from  0 to  <»  along  the  real  axis.  Because  of 

different  time  conventions,  above  is  equal  to  defined  in  [14]. 

« 
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(11)  Residue  series  representation;  For  real  positive  i. 


v(0 

“ 1 

1 " 

(A-6) 

n-1  " 

u(f.) 

« -jet 

r n 

I e 

(A-7) 

n-1 

v,(0 

« -jet' 

r n 

1 e 

(A-8) 

1 

n**l 

v’(0 

. i r- 

1 (1  - j2et;)(t;)  S " 

(A-9) 

4 

n-1 

u'(0 

“ , 9 1 “J^^n 

n-1 

(A-10) 

where  {t  } and  (t'  ) are  zeros  of  w„(t)  and  w’(t),  respectively,  and  they 
n n * * 

are  tabulated  in  (12],  [14]  and  p.  478  of  (20l» 

(ill)  Small  argument  asymptotic  expansion;  For  real  positive  5 and  ^ 0, 


v(C) 


1 . £ _ J.iu  , ,0-356  „ 


512 


(A-11) 


u(0  'V 


1 . 4 * |i  s’ t - 3. 


,-2.6 


701  X 10  ‘•5“  + ... 

(A-12) 


1+25  ^ ^ g-1iT/4^9/2  ^ ^ 555  ^ 10"^5^  + ... 

(A-13) 


v^(0  X ^ ^ . S - 12  64 


v'  (O  ^ ^-j3it/4^1/2  71  ^2  ^ |3^  g-jit/4^7/2  _ 2.435  x + . . . 


1024 


(A-14) 


„'(5)  , I . 2.211  X lo-lfS  + ... 


(A-15) 


u 


(iv)  Wuinerlcal  evaluation;  For  C i Cq.  the  residue  series  representation 

with  the  first  ten  terms  in  the  summation  may  be  used.  For  C 1 5qi  the  small 

argument  asymptotic  expansion  with  the  first  five  terms  may  be  used.  It  has 

been  indicated  in  [11]  that  the  smoothest  crossover  is  obtained  if  Cq  ■ 0.6. 

Ill  the  present  study,  we  set  Cq  - 0.7,  where  the  difference  in  the  two 

reiiresentat  Ions  is  less  than  0.1%  in  magnitude  and  0.9*  in  phase. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c; 

c 

c: 

c 

c 

n 

(; 

c 

c 
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'-'  APPENDIX  B 

&MPUfER  LISTING  FOR  =CALCULATiNG- MUTUAL 
ADMITTANCE  ON vA  c6NE"BASED:'^dN  RAY  TECfflJIQll, 


PROGRAM  . ' 

i SCONE  <INPUT»TAPE<i,0lJTPUT»TAPE5=lNPUT) 


* 

)ic 

t 

t 

t 

t 

* 

* 

)}: 

t 

t 

t 

t 

» 

* 


MUTUAL  ADMITTANCE  OF  SLOTS  ON  A CONE 


by: 


C. 

P. 


W.  LEE 
L.  LAW 
CHANG 


date:  10/18/77 


UNIVERSITY  OF  ILLINOIS 


t 

* 

* 

f.  f 

■M 

!f: 

>K 

* 

* 

* 

* 

* 

* 

* 

* 

it: 

t 


C 
C 
C 
C 

C**fr** 

CH-J-** 

C4lt*H- 

C 

c 

c 

C5-ll*H- 

C 

c 

c 


PRECISION 

LANGUAGE 

MACHINE 


SINGLE 

FORTRAN 

cdc  cyber 
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SUBROUTINE  NEEDED  : FOCK. 
INf-UT  rOKMAT  : FREE. 
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LATEST  REV IS  IN 


10/10/77 


IMf-'LlCl  COMPLEX  <CyH>-Z)  ,REA1.  ( A-CyD-Gf  Kf  0--Y) 

REAL  ANGLE<70) yRl U20)yR2S(n0) 

UIMENGIitN  TITLE!  13 > 

REAL  TN<10)  yTNPKlO) 

REAL  W<7) 

COMMON  /DATAl/  TNy  TNPIyRHOyCI  yC2yF2y  lOF'yCCyRADNyDECi 
COMMON  /CF/  CVFyCUFyCVlFyCVPFyCUPs- 
DATA  TN/2. 33811 , 4. 087, 5y5. 52156.6. 70671 .7.99417. 
i 9 . 02265 .10. 04017 .11. 00052 . U . 9360? .12. 820 ’O/ 

DATA  rTU.E/13C  *)/ 

DATA  TNfT/ I ,010/9. 3. 24020. 4. 82010. 0.16331. 7. 37210. 
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BBT  AVAIWBIE  COPY 


J 8.48849*9,.53545»i6«52'7A<Srii.47566»i2r38479/ 

■ ' ' 

N « 

C 

PI=4i)|!ATAN(l.E6)  ,. 

Pi02=2i)KpI  $ SriRT2*SQRf(2.) 
iiEG=l86,/Pi  $ .AC0Nl=2i/3, 

RADN^PlViBO*. 

f-RINtik,"  INFHJT  A blNE  OF  MESSAGE  ABOUT  THE  CURRENT  JOB* 

RE(AD(a>9ii'-)  title 
.9Ti  FORMAtdiiAIO) 

PRINtiy"  HALF  CONE  ANGLE=*r 

REAfiJlCytHAtHE 

PRINT**’  A1*B1*A2*B2="* 

READ*vAl*BJ.*A2*B2 
PRINT**"  Wl.l*W22--=** 

REAri**Wll*W22 
PRINT**"  i OF  PHI="* 

REAB**NP 

PRINT**"  INPUT  PHKN)*  ONE  ON  EACH  LINE" 
no  9.13  INP».I*NP 
943  READ* * ANGLE <INP) 

PRINT**'  =1  OF  SETS  OF  C"* 

READ** NSC 

f-'RINT**"  INPUT  C*  ONE  SET  ON  EACH  LINE- 
DO  916  INSC«1*NSC 
9i6  READ**R11<INSC)*R22<INSC) 

PRINT**"  INPUT  THE  INTEGRATION  GRIDS* 

PRINT**'  WHICH  ARE  CORRESPONDED  TO  (A1*B1 *A2*B2) ' 

READ** IPO* TP1*IP4* IPS 

PRINT***  THANK  YOU  FOR  YOUR  ACCURATE  INPUT!!!* 

WRITE<6*S1S)TITLE 

WRITE<6*111) 

ill  FORMAT<'  ■*1AX*46(***)  /*1SX* **" *4AX* ***/*15X 

$***  MUTUAL  ADMITTANCE  OF  SLOTS  ON  A CONE  **/*15X****fA4X*"*" 
$/*15X*A6( **" )////) 

WRITE<6*1111) 

nil  F0RMAT</*20(  "**)/* '*'/***  FREQUENCY  } K---^0.6238D  01  a/WAVELENGT 

$H>*/**"/*20< "** ) ) 

C 

ccccccccccccccc 

c 

c 

T[IETIIA«THATHE*RAnN 

SN-S1N(TH[:TMA) 

ATN=-TAN(TIIE)HA) 

ACO~COS<TIIFTIIA) 

WRirE(6.j;:2)THATHE 

222  FORMAT  (/ *20 ("*")■  /*'*•/*•*  GEOMETRY  J HALF  CONE  ANGLL^^*  *F7.2* 

*"  <H[i:G>“/****/20(’**)) 

515  rORMAK "1*  *10X*0A10/10X*UA10///) 

Wl-WlUMvADN 

W2-^^U:i2*RAl.iN 
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BBT  AVAIWBiE  COPY 


39 


WRLtE(6,r  333)  AlVUl  > A2  »B2 

333  f6rMAT(/i-20<-"*'‘)  /r")K‘'/»"*  SLOT  HIHeNSION  <WAVEl;FN6fH>  r Al-“ 

s!.P7.375Xv  “Bl=iSF7.3y  " i * f5X»  ■A2='’’  yF,2.,3r5Xf  “B2=V;,F,7.»3/f  !')K"/r 

i20('*^-*)') 

WRITE  (6 » 444‘)'U:l  ,;L  j U22 

444  F0RMAT</»20C)k")/»')i:"/F"*  SLOT  ORIENTATION  J W1  = ”fF7.2» 

<DFG>"  f5XF."W2=  SF7i2»"  <DEG>"/r*i’/»26(  "sK” ) ) 
WRltE(6F555),IP.2.FlPl 

GGS  .F0RMAt<:/F'20<")i<*  i/f ’iVF'*  INTERGRATION  GRID  { ">I2»"  X " f;I2/ 

$F/“M/F20r":)i<'!)-) 

C1=CEXP  ( (iiMPLX  <0 . EO  F -PI/3  *.  X) 

C2=0EXp(ChipL-X(d.E0F  Pi/4. ).) 

-bCsC2i(<)l!3' 

ANSRal.30S7/TMETI-IA“i'.755+2.,772)(ttHETHA-l,459)KThttHA**2 
ANSL=a . 719S+i . 4<i<)8)i<'fHETHA-1 . 1295)KTHETMA*!t!2+0 . 6566ji:fHETHA;K!K3 
K-*i2'.*Pl 
H:aitH:*ai 

I<(?i2=X)KA2 

'KB2«H:)KB2 

F:2=saRf(Pi-) 

‘W.:cBt)-i3«KBn/iPi 

UrnTH2=KA2/lP2 

;UIDtH3=!KBI7.IP3 

.WIBTH.4=KA1/IP4 

,XLl--KB2/2,rWintHl/2» 

XL2=-KA2/2,-WintH.1/2. 

XL3--KB1 /2 ♦ rWi BtH3/2 . 

XU4“-KAl/2.  ^Wibt.H4/2. 
rib:  100.  -ii---trNsb 

■URITE(6f6A(&) 

666  FbRMAT(////FlOXF  "$$$$$  BATA  OUTPUT  $$$44"') 

KR1=K#RT1<II') 

KR2»K)|<R22(Tr) 

TERMi=KA:l*<<KRl+KHl/2,  ))K)i(2-(KRi-KBl/2.  )**2)/(2-..A'KRi) 

TERH2=KA2*  ( ( kR24KB2/2 . ) iitc2~  < KR2-KB2/2 ) *>|t2 ) / ( 2 , i|cKR2.) 
PARTi=KAl*KBl)t(kA2)(cKB2*SiiRT<ATN/(PI02)ltTERNl*TERM2) ) 

P ART2“i . / < ,30 . #P  I *)K4*KR  1 *kR2*SN ) 

PARt3=SIN  < KBl/2 . ) mSIN  < KB2/2 . ) / ( KBl*KB2/4 . ) 

bnbl=ANSR)t:CbS<ANSL) 

rji:iri2=ANSR#SlN(ANSL) 

ZYTIP=CMPLX < DDDl  f lUin2.) 

$ iKGEXP (CNPLX ( 0 , EO  f PI/4 . -KR1--KR2 ) ) *PART1)*!PART2)KPART3 

$.  *SIN(W1.).*SIN<W2) 

TMAG=CABS(ZYTIP) 

.IFltiiAG.EO.O.  >GOTO  37 
tDB=26.*ALOGiO(TMAG) 

TPHASE*ATAN2  < AIMAG ( Z YtiP ) f REAL ( Z YTIP ) ) *DEG 
37  CONTINUE 

DO  iOO  I^IfNP 

777  FORMAt<//F*  PHI='fF7,2>"  <DEG>  ) ■ f3Xf ‘C1=‘ fE12.7f “ <WAVELENGTH> 
4,>2Xf'C2^-“fE12,7f'  <WAUELEN6TH>  f'  fSXf  ’SLOT  DXST="  fE12f7f 
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iobb: 

20.06. 


t*"'  <WAVl::t.liNGTH>:V/). 
RirX=ANGL§(:lf>*RAi:iN 
IR  ABS I fHul.  .'()):)  PHi-  idl 

. ri6'  'i6uNi=i.rip4', 

■rid",26 

tm=Xl;3iWIti.frt3’«Nij 
'W5=iAfAN20tKZi»TK'yfri 
iRtiSM  bob » 2606, »;2666 
ws^w5-l:pidil- 
cnHfiWuE.' 

.Tl<Yi7.  l«‘f  K'Y  i )k*2+.f.KZl)i(  Y2 


TKR:l=bl.i8 < S«RT CKRl)K>i<2-l  TKYlZl-2 i *KRlYvS«RT (TKYiZl ) JkOOS OWl-WS).) ) 
TPil  11= 1 .♦./SMiKbS i' N,<  SORT  OTK Y 1 Zi;)7tK:i^  1*8 IN  ( W 1 “W5  > >! 

-.DO;  >36'  ;N3=1»;IP2 , 

fKY-2=XU2+UlbfH2*N3 
b,d'  !lO  ^n4“1'»  IPi 
.TRZ2=XL  l.H-W  I DTM 1 *N4 
U6i=ATAN2  aKZIl*  TKY2 ) 

.IF  < WA  ),30b-b » 406o.»  4000 
) 02' 

);CONfiNUE  . 

TlbY2Z2=:i‘KY2**2+TKZ2*>l<2 

fkR2~ADSl80Rt  < l<R2*!i(2+tKY2Z2-2»  *KR2*SQRt  (TKY2Z2 ) *008  ( W2-W6')  > > 
tPHiiJ-PHU  i t/SN*A8]:N  OSOR.f  ( TK.Y2Z2 ) /TKr<2*SIN  ),  > 

.kS«8QRT  < TKR 1 **2-fYKR2**2-2 , *TKRi*TKR2*C0S.<-  ( .TPHl  2-tTPHI  1) 
i'  *SN),) 

dM£GAli=!ASIN  < f KR2*8I N,(  < TPHI 2-TPHIl  i*SN)  /KS ). 
fF.<TKR2**2. Gf..TKRl**2H\T>**2)  .dMEGAl=PI-OMEGAl 
dHEGA2=0HEGA.i'K4rPHI2“tPHiib*8N 
- dMlStN=SIN(dMEGAlb 
dM2SlN=BiN<dHEGA2) 

OM 1 dn8?=cds  ( OMEGA  :io 
0M2C0S=CdS<dMEGA2) 

0>il'2S=AD8<UMl8:iN*0H2SIN) 
dM120--ADS,(  0MlCd8*dM2n()S ) 
iF-(  dM120 . L.f . 3 . E-i> ) UM12C=3 . E-A 
(X)Mi  28 . LT , 3 , E-6 ) drtl2S=3 . E-4 
UPPfi:R=£)dRT<.f'KRi*TKR2->*AfN 
KRt2~l)PPER/dMl2S 
KRn2=UPPER/tiM12d 
kA---ABS<fkRl.*0H.l8IN.*8N*AC:0**2/2.) 

:$  **(  i./3,  )*Ai:!8(tPHI2~'fPHIi) 

.TAIJ=  < < OrtlGIN*OM28lN/A:rN ) **2/ < 2 . *TKR1*TKR2 ) > 

$ **<:i./3,)*K8/l\A 

iF(i\'A.L.t,br/),  do  to  i 
CALL  F-dCKCKA) 

GO'  TO  2 
CALL  FddKHKA) 

Zd=  ( 0 . » 1 , ) *CEXP  ( CMPl.X  ( 0 . EO » -K8 ) ) / ( 240 . *K8*P]  **2  > 
80TAij=SdRT(TAi))  * j AlJltl=tAll**l  ,G 


190 


% 


a 


I 


HB=ZG:t:(.CMF-L)(.(  1 .E6»-i  ,/KS))l<SQf  AUJXCyF-CUF’/KSJitKtiS'CTAt'.ia 
$ HO.  H.  ))l<Cl./(S«Rta#KRT'2-y)>li*ACdi^ii|{Sat^U 

$ HO.  fi . )/.<SdRmj<RT2)*)(cAl:6Ni*KRt2/kRB2!i<tAU3.S)l<(;UIV> 

HT--<26*(0.7 i . )/KS,HSQTAU*GVFfCMPLX(i .Ed>-2./kd:))|<t^^^ 
!t:CUF+<_0..H.)/-<SQR;r2*kRt2')»*ACONl)l<TAUi5- 
$ :i(ClJPF) 

W3--0iiEGAi-TFHIl)|!nN+PI/2 . -Wl 
y^=6MEG^2- (TPHI'2-tPHio  )l<SN+PI/2 . -W2 
ifiREEN^CdS  ( W3 ) >KCf)S  < U4 ) jKHB+SiN  ( W3 ) iK.SIN  ( W4  ).)KHT 
ZSUHsZddM-ffcOS  ( Pl/Kdi)KtkYi ) *G0S  < P17kA2!KTKY2 ) JjtZGREEN. 

40:  Cdi^tlNljE  ... 

■io-  COktiNU'E 
20  GONTiNUE 
10  CQN.tiNUE 

.2;Yi2n«ZdUrt)t!WIDtHl*UIDfH2*WinTH3HtWiriTH4 
$ )/sdRT(kAi*kB.l*KA2»kB2^ 

XHAGi=CABS(Zyi2n). 

PHASE--*ATAN2  ( A IMAG  < Z.Y12D ) r REAL  ( Z Y12n  ) > SBEG 
-DBr26)|iALdGi6<XHAG). 

liiIST=k)((SdRt  < R1  ,i  d I ) *)i(2+R22 (.iiO  ii()*2-2JKRl  t (;l i > i(cR22  ( U h 
$))(cbs<pHr*sN).>, 

SDISTnDldt/K 

WRITE  < 6 , 7-7'/' ) ANGLE  <10  » R1 H 1 1.)  r,R22  ( H ) r SB  I ST 
ZZZZP='ZY12n)KGEkP.<  CHPLX  < P ♦ EO.f  BIST).) 

PHASEN=AtGN2<Ai‘HAG<ZZZ2P)  fREALdZZZP) ) 

PHASEN=PHASEN)i!DEG 

WRITE  U t B8Q ) >:M  AG » PHASE » DB  t PHASER 

SGG  FC)RMAt<2X»’Yi.2B=''/E13.4»'  <HH0>HF7.2i-  • <BEG7T"  f GXHBB™  ^/E12.5? 

1 “{  NORM  PHASE=  NF7.2/  ‘ <BEG>V) 

WRITE  < 6 y 999  > TMAG  y TPHASE  * TBB 

999  F0RMAT<2X»  ''Yl2T="  yE13.4y.“  <MH0>-yF7.2y  " <riEG>y  *■  y5Xy 'IiRu  SE12.6V) 
Zyi2i=ZY.12n+2YTIP 
)<MAG=CABSCZY12) 

PHASE=;ATAN2  < Ai MAG d.Yl 2 ) r REAL ( Z Y1 2 ) ) )KDEG 
DB=2bi*AL0Gi0<XMAG) 

Z'ZZ2P=ZY  1 ZiXCEXP  ( CMPLX  < 0 ♦ EO  y BIST  > > 

PNASEN=At AN2  < AIMAG ( ZZZZP ) r REAL  < ZZZZP ) ) 

PHASEN=PHASEN!((BEG 

WRITE  < 6 y B09  J XM AG  y PHASE  y DB  y PHASER 

889  f 0RMAT<2Xy  ‘Y12  ='  yG13i4y*  <MH0>*yF7.2y  ‘ <DEG>y  “ y 5X,y ' BB=5  VyE12<5y 
1 ‘i  NORM  PHASE=  •.yF7,.2^  ' <DEG>’/) 

iOO  CONTINUE 

WRITE  d > 191 9 ) BDBl y DDD2 

1919  F0RMAT<GXy////5Xy “REAL=*yGi0.2HIMAG=-yG10.2) 

Stop 

END 


lOsU.tJ 
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tNi^,  SUDROiJtlNE  is  USED  fO  CALCULATE  THE  RUNCTIdNS  CyEf CUE fGVlF»Cy.PF» CURE 

SUDRdUfjiNE  FCiCK^;)-  B 

iMRLid.it  'REALiA=^,B»D~H>P  -V*)  rCdHRLEX  (CVZ): 

Rlidh  tNXTOHtNpiiid). 

(yOMMON/CF^f^CVE » bUEf.G.U  J.E>  CURE  r CURE 

COHMON/DAtAl  /TN  r tNPIi  RHO  i clfC2r.F2  r i OR  f CC » RADN  r DEG 

ET-SQRti-X)' 

F3=X:it)|c;ii5 

XiUEsOi 

:dUE==6v 

cOlF=b'. 

CVRF=0*.  . 

CUPF=0,‘. 

iifj  20  N=l-»10 

ZTN=t,N(N')i)?CL 

ZtNP.i=tNPi’(i^))|!Cl 

C3teCEXR<CMPLX:(6.0j-Xy*ZfHPI) 

c4?=CEXF  (CHPLX  ( 0 . 0 » ~xt  iZTN ) 

cOF=icOF;fC3/Z.tNPi 

CUi'£>CUE+C4 

cuiit-cOiF+ca 

CORE-  H ♦0-CMPLX  (0.0  »,2*X ) JfiZTNPI ) )KC3/ZtNPI+CVPF 
CUPF=  <1  i 0-rCMPLX  ( 0 i 0 » 2;.  )KX/3  iO'iteZTN ) !|tC4+CUPF 
l20  -CONTiNilE  ,1, 

CVF=F.2*F1)|!GVF./C2 
:CUF=2V*F2)|<'E3*C2iKCIJF 

CyiF=2 .0(!E2!»!E3i|tC2*C,VlF  >*'■ 

CVPF=F2)t!C0PF/  i2 . >KEi  *02 ) 

CURF-3 . *F2!>:F1  *d2*CUPF 

RStliRN 

ENTRT  EOCKi 

XtHREE«Xj|i*3 

Pi«S(tRt<X) 

F3=XtHREE' 

Z 1 -FaiKCd^SCJRT  ( F3 ) 

Z2=CMPLX  ( 0 . 0 . O/AO . ) ^XTHREE 

Z3=F2iKy()i!*4 . S/  < d2il<64 , ) 

F4=F3**2 

CVE--!  . O-Zl/4  .TZJKZZt? . jkZS/S . -4  ♦ i41Er3*F4 
nUE=l . O-Zl/2 . +25 . )KZ2+5 . *'Z3-3'i  701E-^2*F4 
CV1E=1 . 0+Zl/2 . -35 . *Z2-7  *Z3+4  .•555E,-2*F4 
qVPE= . 375^:Fl>yE2/CC+,2rl , )|fZ2/X+63 . )(cZ3/  ( 16  < *X ) -2 . 485E-2*E4/X 
CUPE== . ZSikEl  S:r2/CC+CHRLX.<  0 . 0 > 1 , 25*X*!l<2 ) +22 , 5StZ3/X 
1 -2<221E^n<E4/X 

RETURN 
END 


